[GRASSLIST:6195] Problems with r.out.bin
Patricio Toledo Peña
patoledo at ing.uchile.cl
Wed Mar 16 18:32:48 EST 2005
Hi, i am a happy user of grass. I think i found a little bug, i
hope.
I have two rasters. A SRTM imported without problems with r.in.srtm
and SPOT image imported from a TIFF. Grass can show me the info
without problems with d.rast.
I would like to make an output to be used with GMT's grdimage. If i
use
[]% r.out.bin -h null=NaN in=srtm out=srtm.grd
and then grdimage, there is no problem!, so i tried with the spot,
but i couldnt render the image.
i decided to use old school output, i use
[]% r.stats -1x in=srtm
and the output looks like this
1 1 1911.8974609375
2 1 1021.7065429688
3 1 683.3014526367
4 1 547.7792358398
5 1 563.8315429688
6 1 697.7428588867
7 1 674.881652832
8 1 885.9045410156
9 1 1811.9166259766
10 1 1894.6561279297
11 1 2184.4465332031
...
and for the spot this
r.stats -1x spot
1 1 10667
2 1 8483
3 1 12844
4 1 8466
5 1 20080
6 1 13835
7 1 12741
8 1 10631
9 1 12743
10 1 6406
11 1 6368
...
so everything looks ok.
After many cups of coffee i managed to look at the bin output. I
tried
[]% r.out.bin -h null=NaN in=spot out=spot.grd
and
[]% grdreformat srtm.grd=1 srtm.nc
[]% grdreformat spot.grd=1 spot.nc
then i use ncdump to look at the grid.
[]% ncdump srtm.nc
looks like this
-----------------CUT--------------------------
netcdf srtm {
dimensions:
side = 2 ;
xysize = 240 ;
variables:
double x_range(side) ;
x_range:units = "Meters ,
\v@¤íÿ_Ð@d@@`@È\f@è\f at X\v\f at l@Ôo@Ë\t@¤íÿ¿Ï@Ë\t@µØ" ;
double y_range(side) ;
y_range:units =
"Meters¿o at X\v\f@`íÿ¿ì\bíÿ¿o@µØ\\îÿ¿lf at 1ú@íÿ¿ @È\f at pëÿ¿" ;
double z_range(side) ;
z_range:units = "elevationM8n@£\bíÿ¿Wû@\b`l at lf@
´îÿ¿z£@Ôo at X\v\f@" ;
double spacing(side) ;
int dimension(side) ;
float z(xysize) ;
z:scale_factor = 1. ;
z:add_offset = 0. ;
z:node_offset = 1 ;
// global attributes:
:title = "dem_srtm_90mx@\r\blf at lf@Âc@`x@\n",
"@@@lf at l@" ;
:source = "r.out.bin -h input=dem_srtm_90m
output=srtm.grd0 used for NULL" ;
data:
x_range = 310000, 370000 ;
y_range = 6235000, 6335000 ;
z_range = 0, 2804.87915039062 ;
spacing = 5000, 5000 ;
dimension = 12, 20 ;
z = 1911.897, 1021.707, 683.3015, 547.7792, 563.8315, 697.7429,
674.8817,
885.9045, 1811.917, 1894.656, 2184.447, 0, 1524.52, 1867.924,
663.9888,
602.0908, 519.0523, 765.1163, 594.4846, 706.0734, 1212.255,
1479.04,
2702.787, 0, 1720.611, 832.5412, 738.3867, 559.8695, 496.3092,
648.4741,
564.0004, 685.2112, 1330.027, 1833.741, 2804.879, 0, 764.7325,
1113.144,
884.9847, 515.4147, 488.7629, 506.5552, 530.0245, 599.58,
1238.697,
1361.697, 2373.276, 0, 1638.143, 1045.054, 740.9617, 486.9698,
483.3979,
488.5664, 747.3654, 615.978, 1140.735, 971.7466, 1293.763, 0,
683.3162,
817.0809, 625.1544, 475.649, 476.8131, 487.9839, 501.0298,
639.7838,
1023.9, 820.8611, 896.476, 0, 596.2768, 1182.775, 597.7131,
467.4986,
477.4478, 488.2751, 515.9766, 540.3224, 664.6289, 763.2197,
1317.51, 0,
208.3247, 852.842, 495.4633, 469.7721, 473.3067, 485.7872,
517.0618,
564.0345, 605.4836, 688.7662, 1530.948, 0, 188.037, 269.8225,
844.0717,
466.8229, 455.4855, 482.079, 523.5785, 556.554, 587.5854,
782.3171,
1458.525, 0, 641.4139, 439.9863, 402.5107, 472.7809, 452.5125,
499.4549,
550.4744, 594.2014, 630.0211, 852.4738, 2028.455, 0, 216.0057,
361.0359,
548.9712, 403.6217, 452.0291, 716.2377, 574.1116, 621.963,
672.7051,
776.1525, 826.4809, 0, 569.1906, 673.9302, 368.0906, 417.1798,
468.4817,
518.9277, 574.379, 624.7796, 677.2015, 728.1828, 1031.293, 0,
276.125,
311.6169, 350.5661, 391.2746, 792.7806, 511.8949, 557.0769,
778.4391,
668.5333, 702.291, 1478.079, 0, 254.2038, 279.732, 311.3417,
356.5961,
768.6725, 485.8281, 512.5512, 786.2317, 806.8865, 746.6079,
906.1571, 0,
602.7701, 583.3986, 331.9892, 372.4999, 410.9167, 444.7338,
460.9563,
458.8165, 1277.886, 877.2421, 1959.047, 0, 812.3603, 722.5879,
809.0806,
354.405, 373.225, 395.8081, 410.2335, 405.3456, 663.4382,
1451.57,
2679.614, 0, 808.3639, 452.7131, 364.8825, 362.1678, 526.1887,
376.366,
378.2451, 418.3407, 616.1348, 1688.328, 2589.296, 0, 1906.097,
1647.922,
810.8583, 578.948, 841.3588, 395.6625, 838.9558, 425.5698,
953.824,
1469.772, 1608.909, 0, 537.739, 1550.32, 1003.811, 1071.668,
1459.342,
884.1784, 465.0422, 584.5891, 1117.134, 2319.113, 2493.716, 0,
866.7081,
1447.663, 1515.252, 1140.653, 1452.902, 798.4017, 486.7825,
536.4229,
727.7171, 1633.92, 1886.89, 0 ;
}
---------------------------------------------
[]% ncdump spot.nc
like this
-----------------------CUT---------------------------
netcdf spot {
dimensions:
side = 2 ;
xysize = 240 ;
variables:
double x_range(side) ;
x_range:units = "Meters ,
\v@¤íÿ_Ð@d@@`@È\f@è\f at X\v\f at l@Ôo@Ë\t@¤íÿ¿Ï@Ë\t@µØ" ;
double y_range(side) ;
y_range:units =
"Meters¿o at X\v\f@`íÿ¿ì\bíÿ¿o@µØ\\îÿ¿lf at 1ú@íÿ¿ @È\f at pëÿ¿" ;
double z_range(side) ;
z_range:units = "elevationM8n@£\bíÿ¿Wû@\b`l at lf@
´îÿ¿z£@Ôo at X\v\f@" ;
double spacing(side) ;
int dimension(side) ;
float z(xysize) ;
z:scale_factor = 1. ;
z:add_offset = 0. ;
z:node_offset = 1 ;
// global attributes:
:title = "imag_spot_0`x@\r\blf at lf@Âc@`x@\n",
"@@@lf at l@" ;
:source = "r.out.bin -h input=imag_spot_0
output=spot.grd0 usedfor NULL" ;
data:
x_range = 310000, 370000 ;
y_range = 6235000, 6335000 ;
z_range = 8.79034526670958e-42, 1.47562183307065e+36 ;
spacing = 5000, 5000 ;
dimension = 12, 20 ;
z = 5.52817e-19, 4.953313e-19, 2.075827e-06, 6.003838e-14,
6.93769e-24,
8.923469e-42, 2.395605e-16, 2.536351e-06, 5.854575e-09,
2.580565e-06,
2.185477e-21, 8.923469e-42, 5.948305e-24, 8.041159e-09,
6.814512e-09,
1.988653e-16, 0.0009047714, 1.493924e-41, 3.990042e-09,
7.512333e-14,
1.293629e-16, 26386.59, 1.802539e-21, 1.19783e-41, 5.121541e-24,
1.424025e-16, 6.119387e-09, 1.924446e-11, 6.094157e-14,
8.839391e-42,
2.695473e-11, 0.236614, 3.568534e-09, 0.00077581, 7683.11,
1.329552e-41,
4.965599e-24, 0.1427695, 5.223165e-19, 0.0003209909,
1.911767e-11,
1.33768e-41, 6.88098e-24, 2.042624e+09, 9.468347e+13,
1.151268e+09,
0.0006491734, 8.980922e-42, 0.0009161414, nan, 5.846576e+33,
28.90538,
1.736398e-16, 8.931876e-42, 2.308851e-16, 8.589542e-24,
1.23291e-16,
2741403, 5.824033e-24, 1.647086e-41, 6.988471e-24, 1.675731e-21,
6.835183e-24, 1.634152e+31, 4.740647e-14, 8.934679e-42,
4.141374e-24,
5.969695e-19, 1.688843e-21, 1.682113e-21, 5.382549e-24,
8.790345e-42,
4.067404e-19, 5.458837e-24, 9.093348e-07, 1.763027e-21,
1.501376e-16,
1.04635e-41, 3.178451e-14, 1.319227e-16, 4.444136e-14,
2.596219e-06,
5.255992e-19, 8.839391e-42, 6.984338e-24, 1.957237e-11,
1.381501e-16,
5.302054e-24, 5.225773e-24, 4.589813e-41, 6.051259e-24,
5.508448e-24,
51.29015, 1.41027e-11, 3.365455e-24, 8.833786e-42, 3.440471e-24,
5.921985e-24, 4.238196e-19, 1.684212e-11, 6.677936e-19,
1.933512e-41,
11271.27, 4.991422e-24, 5.430862e-24, 1.475622e+36,
5.017299e-24,
8.929074e-42, 5.172394e-24, 5.094835e-24, 3.31417e-24,
1.582978e-21,
7.376226e-14, 3.266707e-41, 4.293486e-24, 5.096956e-24,
7.143683e-24,
nan, 4.965612e-24, 1.638398e-41, 5.172394e-24, 5.094835e-24,
3.31417e-24,
1.582978e-21, 7.376226e-14, 3.266707e-41, 4.293486e-24,
5.096956e-24,
7.143683e-24, nan, 4.965612e-24, 1.638398e-41, 5.172394e-24,
5.094835e-24, 3.31417e-24, 1.582978e-21, 7.376226e-14,
3.266707e-41,
4.293486e-24, 5.096956e-24, 7.143683e-24, nan, 4.965612e-24,
1.638398e-41, 5.172394e-24, 5.094835e-24, 3.31417e-24,
1.582978e-21,
7.376226e-14, 3.266707e-41, 4.293486e-24, 5.096956e-24,
7.143683e-24,
nan, 4.965612e-24, 1.638398e-41, 5.172394e-24, 5.094835e-24,
3.31417e-24,
1.582978e-21, 7.376226e-14, 3.266707e-41, 4.293486e-24,
5.096956e-24,
7.143683e-24, nan, 4.965612e-24, 1.638398e-41, 5.172394e-24,
5.094835e-24, 3.31417e-24, 1.582978e-21, 7.376226e-14,
3.266707e-41,
4.293486e-24, 5.096956e-24, 7.143683e-24, nan, 4.965612e-24,
1.638398e-41, 5.172394e-24, 5.094835e-24, 3.31417e-24,
1.582978e-21,
7.376226e-14, 3.266707e-41, 4.293486e-24, 5.096956e-24,
7.143683e-24,
nan, 4.965612e-24, 1.638398e-41, 5.172394e-24, 5.094835e-24,
3.31417e-24,
1.582978e-21, 7.376226e-14, 3.266707e-41, 4.293486e-24,
5.096956e-24,
7.143683e-24, nan, 4.965612e-24, 1.638398e-41, 5.172394e-24,
5.094835e-24, 3.31417e-24, 1.582978e-21, 7.376226e-14,
3.266707e-41,
4.293486e-24, 5.096956e-24, 7.143683e-24, nan, 4.965612e-24,
1.638398e-41, 5.172394e-24, 5.094835e-24, 3.31417e-24,
1.582978e-21,
7.376226e-14, 3.266707e-41, 4.293486e-24, 5.096956e-24,
7.143683e-24,
nan, 4.965612e-24, 1.638398e-41, 5.172394e-24, 5.094835e-24,
3.31417e-24,
1.582978e-21, 7.376226e-14, 3.266707e-41, 4.293486e-24,
5.096956e-24,
7.143683e-24, nan, 4.965612e-24, 1.638398e-41 ;
}
-------------------------------------
so, the data is spreaded between
8.79034526670958e-42, 1.47562183307065e+36
and i think that is the problem. I think i could use oldschool
output and then xyz2grd but is very interesting these possible bug.
As an additional info i tried grass53 and grass60 (from debian
unstable) and is the same problem.
Any ideas?
Thanks a lot.
--
Patricio Toledo Peña
Departamento de Geofísica
Universidad de Chile
More information about the grass-user
mailing list