On Tue, Mar 16, 2021 at 9:51 AM Markus Neteler <neteler@osgeo.org> wrote:
Hi,
On Sun, Mar 14, 2021 at 3:32 PM javier Garcia Prieto
<fjgarciaprieto@hotmail.com> wrote:
>
> Thank you!
> You are really great! This is going to be so much useful to Mexican people, You are very kind.
Glad to help. Now it is merged, i.e. it is already available in the
current source code of GRASS GIS 7.8.git and will be part of GRASS GIS
7.8.6 and later.
BTW, I made a test which you indicated to me off-list, with the
udpated GRASS GIS 7.8.git version:
gdalinfo republica2017_r15m.bil
Driver: EHdr/ESRI .hdr Labelled
Files: republica2017_r15m.bil
republica2017_r15m.bil.ovr
republica2017_r15m.bil.aux.xml
republica2017_r15m.hdr
republica2017_r15m.stx
republica2017_r15m.prj
Size is 218763, 130943
Coordinate System is:
GEOGCRS["ITRF92",
DATUM["International Terrestrial Reference Frame 1992",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]],
ID["EPSG",6651]],
PRIMEM["Greenwich",0,
ANGLEUNIT["Degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["longitude",east,
ORDER[1],
ANGLEUNIT["Degree",0.0174532925199433]],
AXIS["latitude",north,
ORDER[2],
ANGLEUNIT["Degree",0.0174532925199433]]]
Data axis to CRS axis mapping: 1,2
Origin = (-117.124268882302502,32.718626154577500)
Pixel Size = (0.000138888889000,-0.000138888889000)
Corner Coordinates:
Upper Left (-117.1242689, 32.7186262) (117d 7'27.37"W, 32d43' 7.05"N)
Lower Left (-117.1242689, 14.5320984) (117d 7'27.37"W, 14d31'55.55"N)
Upper Right ( -86.7405189, 32.7186262) ( 86d44'25.87"W, 32d43' 7.05"N)
Lower Right ( -86.7405189, 14.5320984) ( 86d44'25.87"W, 14d31'55.55"N)
Center (-101.9323939, 23.6253623) (101d55'56.62"W, 23d37'31.30"N)
Band 1 Block=218763x1 Type=Int16, ColorInterp=Undefined
Min=-32.000 Max=5610.000
Minimum=-32.000, Maximum=5610.000, Mean=1041.610, StdDev=820.476
NoData Value=-32768
Overviews: 109382x65472, 54691x32736, 27346x16368, 13673x8184,
6837x4092, 3419x2046, 1710x1023, 855x512, 428x256, 214x128
Metadata:
SourceBandIndex=0
STATISTICS_SKIPFACTORX=1
STATISTICS_SKIPFACTORY=1
# generate location from dataset
grass78 -c republica2017_r15m.bil /mnt/data/grassdata/mexico
# verify projection
GRASS 7.8.6dev (mexico):~ > g.proj -w
GEOGCRS["ITRF92",
DATUM["International Terrestrial Reference Frame 1992",
ELLIPSOID["GRS 1980",6378137,298.257222101,
LENGTHUNIT["metre",1]],
ID["EPSG",6651]],
PRIMEM["Greenwich",0,
ANGLEUNIT["Degree",0.0174532925199433]],
CS[ellipsoidal,2],
AXIS["longitude",east,
ORDER[1],
ANGLEUNIT["Degree",0.0174532925199433]],
AXIS["latitude",north,
ORDER[2],
ANGLEUNIT["Degree",0.0174532925199433]]]
# import DEM
r.import republica2017_r15m.bil output=republica2017_r15m
# set color table
r.colors republica2017_r15m color=srtm_plus
The result looks good to me.
best,
Markus