::
:: rename csv2shp to csv2shp.bat
::
  echo on
  set pp=./20200930/
  set fn=ArcRCC_regions_20200930
  set sql="SELECT CAST(FID_cnt AS character(1)), CAST (name AS character(26)) FROM %fn%"

  if exist %pp%%fn%*.shp del %pp%%fn%*.shp %pp%%fn%*.shx %pp%%fn%*.dbf %pp%%fn%*.prj %pp%%fn%*.tif

  ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:4326 %pp%%fn%_geo.shp %pp%%fn%.csv -sql %sql%
  ogr2ogr -s_srs EPSG:4326 -t_srs EPSG:3995 %pp%%fn%_stere.shp %pp%%fn%.csv -sql %sql%
  gdal_rasterize -a FID_cnt -te -4700000 -3500000 3500000 4700000 -tr 10000 10000 -a_nodata 99 -co "COMPRESS=LZW" %pp%%fn%_stere.shp %pp%%fn%_stere.tif