issue_comments: 503052359

This data as json

html_url issue_url id node_id user created_at updated_at author_association body reactions performed_via_github_app issue
https://github.com/pydata/xarray/issues/501#issuecomment-503052359 https://api.github.com/repos/pydata/xarray/issues/501 503052359 MDEyOklzc3VlQ29tbWVudDUwMzA1MjM1OQ== 37119533 2019-06-18T10:49:52Z 2019-06-18T10:49:52Z NONE

I am doing a project where we need to segment a geospatial data. The data set comprises the Fire Weather Index (FWI) -- a fire hazard classification system -- for the whole world from the 1980's to 2018. We want to segment data according to the NUTS (or GADM) in order to be able to predict and investigate possible patterns over time.

We have to our disposal:

  • the FWI dataset (resolution: 0.7 degrees ~ 80x80 Km)
  • the NUTS (and GADM) shape files

Our methodology was the following:

  1. import shape files with geopandas
  2. import the FWI dataset with xarray
  3. rasterize the shapefiles (adapting @shoyer code)
  4. get FWI data segmented by NUTS/GADM

A couple of issues were identified (as shown in the picture below):

As can be seen there seems to exist an offset between the segmented data and the actual map areas. Moreover, there are some islands that are not even captured by the segmentation.

The only modification to the original code was the introduction of the argument all_touched=True to the features.rasterize method.

Does anyone ever got a similar error?

 
{
    "total_count": 0,
    "+1": 0,
    "-1": 0,
    "laugh": 0,
    "hooray": 0,
    "confused": 0,
    "heart": 0,
    "rocket": 0,
    "eyes": 0
}
   98074194