[GRASS-dev] [GRASS GIS] #2757: r.import: ERROR: Input raster map is outside current region

[GRASS GIS]

#2757: r.import: ERROR: Input raster map is outside current region
-------------------------+-------------------------------------------------
  Reporter:  neteler     |      Owner:  grass-dev@…
      Type:  defect      |     Status:  closed
  Priority:  normal      |  Milestone:  7.4.0
 Component:  Python      |    Version:  svn-releasebranch74
Resolution:  fixed       |   Keywords:  r.import, r.in.gdal, r.external,
       CPU:              |  r.region
  Unspecified            |   Platform:  All
-------------------------+-------------------------------------------------
Changes (by neteler):

 * status:  reopened => closed
 * version:  svn-trunk => svn-releasebranch74
 * resolution:   => fixed
 * milestone:  7.0.6 => 7.4.0


Comment:

 Checking this ticket with GRASS GIS 7.4:

 {{{
 # Fedora box
 GRASS 7.4.1svn (nc_basic_spm_grass7):/scratch > g.version -ger
 version=7.4.1svn
 date=2018
 revision=r72180
 build_date=2018-01-28
 build_platform=x86_64-pc-linux-gnu
 build_off_t_size=8
 libgis_revision=70829
 libgis_date="2017-04-04 09:43:02 +0200 (Tue, 04 Apr 2017) "
 proj4=4.9.3
 gdal=2.1.3
 geos=3.6.1
 sqlite=3.20.1
 }}}


 Test in Long-Lat:

 {{{
 GRASS 7.4.1svn (latlong_wgs84):/scratch > r.import input=bio1.bil
 output=bioclim1
 90 degree north is exceeded by 6.82121e-13 cells
 Importing raster map <bioclim1>...
  100%
 GRASS 7.4.1svn (latlong_wgs84):/scratch > g.region raster=bioclim1 -p
 projection: 3 (Latitude-Longitude)
 zone:       0
 datum:      wgs84
 ellipsoid:  wgs84
 north:      90N
 south:      60S
 west:       180W
 east:       180E
 nsres:      0:02:30
 ewres:      0:02:30
 rows:       3600
 cols:       8640
 cells:      31104000

 GRASS 7.4.1svn (latlong_wgs84):/scratch > g.region raster=bioclim1 -g
 projection=3
 zone=0
 n=90
 s=-60
 w=-180
 e=180
 nsres=0.0416666666666667
 ewres=0.0416666666666667
 rows=3600
 cols=8640
 cells=31104000
 }}}

 Test in NC sample dataset location:

 {{{
 GRASS 7.4.1svn (nc_basic_spm_grass7):/scratch > r.import input=bio1.bil
 output=bioclim1 resample=bilinear extent=region resolution=region
 WARNING: Projection of dataset does not appear to match current location.

          Location PROJ_INFO is:
          name: Lambert Conformal Conic
          proj: lcc
          datum: nad83
          a: 6378137.0
          es: 0.006694380022900787
          lat_1: 36.16666666666666
          lat_2: 34.33333333333334
          lat_0: 33.75
          lon_0: -79
          x_0: 609601.22
          y_0: 0
          no_defs: defined

          Dataset PROJ_INFO is:
          name: WGS 84
          datum: wgs84
          ellps: wgs84
          proj: ll
          no_defs: defined

          ERROR: proj

          In case of no significant differences in the projection
          definitions, use the -o flag to ignore them and use current
          location definition.
          Consider generating a new location from the input dataset using
          the 'location' parameter.
 90 degree north is exceeded by 6.82121e-13 cells
 Importing raster map <bioclim1>...
  100%
 Estimated target resolution for input band <bioclim1>: 4021.78765331
 Using current region resolution for input band <bioclim1>: nsres=500.0,
 ewres=500.0
 Reprojecting <bioclim1>...

 GRASS 7.4.1svn (nc_basic_spm_grass7):/scratch > g.region -p
 raster=bioclim1
 projection: 99 (Lambert Conformal Conic)
 zone:       0
 datum:      nad83
 ellipsoid:  a=6378137 es=0.006694380022900787
 north:      320000
 south:      10000
 west:       120000
 east:       935000
 nsres:      500
 ewres:      500
 rows:       620
 cols:       1630
 cells:      1010600
 }}}

 (the "ERROR: proj" is a bit disturbing as a message, since using r.import)

 It looks ok to me - thanks for the hard work esp. by Markus Metz on
 improving the global data import and treatment of global data with
 coordinate precision issues!

 Closing, this improvement is part of the recently published GRASS GIS
 7.4.0.

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/2757#comment:21>
GRASS GIS <https://grass.osgeo.org>