Using: Raster Spatial Allocator Version 4.2, 05/2014 Compute grid landuse percentage information... GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Getting modeling grid domain variables... Grid name: 2017test5 Rows=42 Cols=60 Total gird cells=2520 xmin=-450000.000000 ymin=-1062000.000000 xcell=15000 ycell=15000 Grid domain projection proj4 = +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 proj type = 1 Domain BELD4 netCDF file: /home/u220220948411/festc/scenarios/china5/share_data/beld4_2017test5_2011_5.nc U.S. county shapefile: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/country_xiugai_lcc6.shp Region country shapefile: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/province_xiugai_lcc6.shp 8-digit HUC shapefile: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/china_basins_lcc7.shp region elevation image: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img region slope image: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img Output text file name: /home/u220220948411/festc/scenarios/china5/share_data/EPICSites_Info14.csv Output text file name: /home/u220220948411/festc/scenarios/china5/share_data/EPICSites_Crop14.csv Reading grid landuse data to select EPIC sites... Reading BELD4 netCDF file: /home/u220220948411/festc/scenarios/china5/share_data/beld4_2017test5_2011_5.nc Obtaining all dimension IDs in input NetCDF file... BELD4 file has: 8 dims, 16 variables, 32 global attributes, 0 unlimited variable ID BELD4 file dimension: dimName = Time dimSize = 1 BELD4 file dimension: dimName = DateStrLen dimSize = 19 BELD4 file dimension: dimName = west_east dimSize = 60 BELD4 file dimension: dimName = south_north dimSize = 42 BELD4 file dimension: dimName = land_cat dimSize = 40 BELD4 file dimension: dimName = crop_cat dimSize = 42 BELD4 file dimension: dimName = tree_cat dimSize = 194 BELD4 file dimension: dimName = nameStrLen dimSize = 50 Getting IOAPI file grid information... Rows=42 Cols=60 Total gird cells=2520 xmin=-450000.000000 ymin=-1062000.000000 xmax=0.000000 ymax=-nan xcell=15000.00000000 ycell=15000.00000000 proj4=+proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30.000 +lat_2=60.000 +lat_0=40.000 +lon_0=100.800 Variable name = LANDUSEF var_type = 5 ndims=4 Dimensions: 1 40 42 60 Total size = 100800 Obtained variable: LANDUSEF Reading BELD4 netCDF file: /home/u220220948411/festc/scenarios/china5/share_data/beld4_2017test5_2011_5.nc Obtaining all dimension IDs in input NetCDF file... BELD4 file has: 8 dims, 16 variables, 32 global attributes, 0 unlimited variable ID BELD4 file dimension: dimName = Time dimSize = 1 BELD4 file dimension: dimName = DateStrLen dimSize = 19 BELD4 file dimension: dimName = west_east dimSize = 60 BELD4 file dimension: dimName = south_north dimSize = 42 BELD4 file dimension: dimName = land_cat dimSize = 40 BELD4 file dimension: dimName = crop_cat dimSize = 42 BELD4 file dimension: dimName = tree_cat dimSize = 194 BELD4 file dimension: dimName = nameStrLen dimSize = 50 Getting IOAPI file grid information... Rows=42 Cols=60 Total gird cells=2520 xmin=-450000.000000 ymin=-1062000.000000 xmax=0.000000 ymax=-nan xcell=15000.00000000 ycell=15000.00000000 proj4=+proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30.000 +lat_2=60.000 +lat_0=40.000 +lon_0=100.800 Variable name = CROPF var_type = 5 ndims=4 Dimensions: 1 42 42 60 Total size = 105840 Obtained variable: CROPF Minimum crop percent used for EPIC site selection: 0.071944 Get dimention size of land_cat from file: /home/u220220948411/festc/scenarios/china5/share_data/beld4_2017test5_2011_5.nc Dimension of land_cat is: 40 Selected 1690 EPIC potential sites. Final sites will be selectled based on each crop fraction > 0.071944 Extracting polygon shapefile attributes for EPIC sites: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/country_xiugai_lcc6.shp... GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Polygon Shapefile projection: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs Grid domain projection: +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 No +a= in GRID_PROJ - +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs and set it to: -9999.000000 EPIC site x and y will be projected into the polygon shapefile to get attributes Completed extracting attributes from shapefile:/home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/country_xiugai_lcc6.shp Extracting polygon shapefile attributes for EPIC sites: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/province_xiugai_lcc6.shp... GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Polygon Shapefile projection: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs Grid domain projection: +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 No +a= in GRID_PROJ - +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs and set it to: -9999.000000 EPIC site x and y will be projected into the polygon shapefile to get attributes Completed extracting attributes from shapefile:/home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/province_xiugai_lcc6.shp Extracting polygon shapefile attributes for EPIC sites: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/china_basins_lcc7.shp... GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Polygon Shapefile projection: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs Grid domain projection: +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 No +a= in GRID_PROJ - +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs and set it to: -9999.000000 EPIC site x and y will be projected into the polygon shapefile to get attributes Completed extracting attributes from shapefile:/home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/china_basins_lcc7.shp Obtaining point values from image file: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Obtaining image information from /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img... Driver: HFA/Erdas Imagine Images (.img) Size is 6303x4480x1 Origin = (-3370785.005894,1926820.206051) Pixel Size = (1199.028,1199.028) Satellite image extent: minXY(-3370785.006,-3444825.837) maxXY(4186689.327,1926820.206) Proj4 projection string from image is: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs Grid domain projection: +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 Image projection: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs No +a= in GRID_PROJ - +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs and set it to: -9999.000000 Image file has different projection and reprojecting it... Project an image data file... Projected image data is stored in file: tmp1159602976_img.bil /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/gdalwarp -s_srs "+proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs " -t_srs "+proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8" -te -450000.000000 -1062000.000000 450000.000000 -432000.000000 -tr 15000.0Creating output file that is 60P x 42L. Processing input file /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img. Using internal nodata values (e.g. -32768) for image /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img. Copying nodata values from source /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img to destination tmp1159602976_img.bil. 0...10...20...30...40...50...60...70...80...90...100 - done. Creating output file that is 60P x 42L. Processing input file /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img. Using internal nodata values (e.g. -9999) for image /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img. Copying nodata values from source /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img to destination tmp1159602976_img.bil. 0...10...20...30...40...50...60...70...80...90...100 - done. 000 15000.0000 /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/10.img tmp1159602976_img.bil Completed projecting the image data. Obtaining image information from tmp1159602976_img.bil... Driver: GTiff/GeoTIFF Size is 60x42x1 Origin = (-450000.000000,-432000.000000) Pixel Size = (15000.000,15000.000) Satellite image extent: minXY(-450000.000,-1062000.000) maxXY(450000.000,-432000.000) Proj4 projection string from image is: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs Rows=42 Cols=60 Total gird cells=2520 xmin=-450000.000000 ymin=-1062000.000000 xmax=450000.000000 ymax=-432000.000000 xcell=15000.00000000 ycell=15000.00000000 proj4=+proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs Image Type = 3 Raster file exists and delete it: tmp1159602976_img.bil Completed value extraction from the image. Obtaining point values from image file: /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img GDAL bin directory: /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/ Obtaining image information from /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img... Driver: HFA/Erdas Imagine Images (.img) Size is 6144x4367x1 Origin = (-3370820.403228,1926842.576550) Pixel Size = (1230.023,1230.023) Satellite image extent: minXY(-3370820.403,-3444669.630) maxXY(4186443.393,1926842.577) Proj4 projection string from image is: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs Grid domain projection: +proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 Image projection: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs No +a= in GRID_PROJ - +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs and set it to: -9999.000000 Image file has different projection and reprojecting it... Project an image data file... Projected image data is stored in file: tmp1159602976_img.bil /home/u220220948411/festc-i/sa_052014/src/libs/gdal-1.11.0/local/bin/gdalwarp -s_srs "+proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs " -t_srs "+proj=lcc +a=6370000.0 +b=6370000.0 +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8" -te -450000.000000 -1062000.000000 450000.000000 -432000.000000 -tr 15000.0000 15000.0000 /home/u220220948411/festc-i/festc1.4/epic/common_data/china_gisfile5/slope_test17.img tmp1159602976_img.bil Completed projecting the image data. Obtaining image information from tmp1159602976_img.bil... Driver: GTiff/GeoTIFF Size is 60x42x1 Origin = (-450000.000000,-432000.000000) Pixel Size = (15000.000,15000.000) Satellite image extent: minXY(-450000.000,-1062000.000) maxXY(450000.000,-432000.000) Proj4 projection string from image is: +proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs Rows=42 Cols=60 Total gird cells=2520 xmin=-450000.000000 ymin=-1062000.000000 xmax=450000.000000 ymax=-432000.000000 xcell=15000.00000000 ycell=15000.00000000 proj4=+proj=lcc +lat_1=30 +lat_2=60 +lat_0=40 +lon_0=100.8 +x_0=0 +y_0=0 +a=6370000 +b=6370000 +units=m +no_defs Image Type = 3 Raster file exists and delete it: tmp1159602976_img.bil Completed value extraction from the image. Build geolocation point trees. Filling HUC8 numbers for sites in USA without 8-digit HUC numbers... Filling REG10,STFIPS,CNTYFIPS for sites in USA without the codes... Filling COUNTRY codes with water/agua/d'eau... Write EPIC site crop file: /home/u220220948411/festc/scenarios/china5/share_data/EPICSites_Crop14.csv... Write EPIC site info file: /home/u220220948411/festc/scenarios/china5/share_data/EPICSites_Info14.csv... Completed in computing EPIC site information. ==== Finished crop site info generation.