Switch example shpfile for cropped region

spireglobal · May 14, 2020 · 0fc4755 · 0fc4755
1 parent f211fa1
commit 0fc4755
Show file tree

Hide file tree

Showing 10 changed files with 9 additions and 9 deletions.
diff --git a/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_aviation_data_in_area.py
@@ -13,7 +13,7 @@
 
 # Load the shapefile area
 driver = ogr.GetDriverByName('ESRI Shapefile')
-shpfile = driver.Open('shpfile/france.shp')
+shpfile = driver.Open('shpfile/italy.shp')
 AREA = shpfile.GetLayer()
 
 # Check if point is inside of shapefile area
@@ -55,8 +55,8 @@ def parse_data(filepath):
     # Retrieve flight level values
     flightlevels = df.index.get_level_values('lv_AMSL0')
     # Filter to a specific flight level:
-    # FL100 = 10000 feet = 3048 meters
-    df = df.loc[(flightlevels == 3048)]
+    # FL450 = 45000 feet = 13716 meters
+    df = df.loc[(flightlevels == 13716)]
     # Get longitude values from index
     lons = df.index.get_level_values('lon_0')
     # Map longitude range from (0 to 360) into (-180 to 180)
@@ -98,10 +98,10 @@ def plot_data(data):
         c=color,
         s=10,
         cmap='Spectral_r',
-        edgecolors='gray',
+        # edgecolors='gray',
         linewidths=0.1
     )
-    plt.title('Clear-Air Turbulence % at FL100')
+    plt.title('Clear-Air Turbulence % at FL450')
     plt.colorbar()
     plt.show()
 

diff --git a/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_precip_data_in_area.py
@@ -13,7 +13,7 @@
 
 # Load the shapefile area
 driver = ogr.GetDriverByName('ESRI Shapefile')
-shpfile = driver.Open('shpfile/france.shp')
+shpfile = driver.Open('shpfile/italy.shp')
 AREA = shpfile.GetLayer()
 
 # Check if point is inside of shapefile area

diff --git a/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py b/examples/working_with_grib_data/cropped_area/get_upper_air_data_in_area.py
@@ -13,7 +13,7 @@
 
 # Load the shapefile area
 driver = ogr.GetDriverByName('ESRI Shapefile')
-shpfile = driver.Open('shpfile/france.shp')
+shpfile = driver.Open('shpfile/italy.shp')
 AREA = shpfile.GetLayer()
 
 # Check if point is inside of shapefile area
@@ -100,7 +100,7 @@ def plot_data(data):
         c=color,
         s=10,
         cmap='coolwarm',
-        edgecolors='gray',
+        # edgecolors='gray',
         linewidths=0.1
     )
     plt.title('Temperature (C) at 20 hPa')

diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.prj b/examples/working_with_grib_data/cropped_area/shpfile/france.prj
diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.shp b/examples/working_with_grib_data/cropped_area/shpfile/france.shp
diff --git a/examples/working_with_grib_data/cropped_area/shpfile/france.shx b/examples/working_with_grib_data/cropped_area/shpfile/france.shx
diff --git a/...grib_data/cropped_area/shpfile/france.dbf → ..._grib_data/cropped_area/shpfile/italy.dbf b/...grib_data/cropped_area/shpfile/france.dbf → ..._grib_data/cropped_area/shpfile/italy.dbf
diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.prj b/examples/working_with_grib_data/cropped_area/shpfile/italy.prj
@@ -0,0 +1 @@
+GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]
diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.shp b/examples/working_with_grib_data/cropped_area/shpfile/italy.shp
diff --git a/examples/working_with_grib_data/cropped_area/shpfile/italy.shx b/examples/working_with_grib_data/cropped_area/shpfile/italy.shx
Original file line number	Diff line number	Diff line change
		@@ -0,0 +1 @@
		GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]