[GRASS-dev] Strange import of .asc

Hi all.
Importing an .asc file (with -o) results in a raster which is not aligend with the
original. When vectorizing it, the vector is aligned with the original, not with the
imported one.
Original file and GRASS Location, with imported raster and resulting vector, here:
https://www.dropbox.com/s/w60ywwddl3aj03h/cavallini.zip
Anyone confirms? If so I'll open a ticket.
All the best.
--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

Hi Paolo,

Importing an .asc file (with -o) results in a raster which is not aligend

with the

original.
Anyone confirms?

gdalinfo 23B07.asc
Driver: AAIGrid/Arc/Info ASCII Grid
Files: 23B07.asc
Size is 1640, 1240
Coordinate System is `'
Origin = (1568780.000000000000000,4896020.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
  AREA_OR_POINT=Point
Corner Coordinates:
Upper Left ( 1568780.000, 4896020.000)
Lower Left ( 1568780.000, 4894780.000)
Upper Right ( 1570420.000, 4896020.000)
Lower Right ( 1570420.000, 4894780.000)
Center ( 1569600.000, 4895400.000)
Band 1 Block=1640x1 Type=Float32, ColorInterp=Undefined
  NoData Value=-9999

and

r.info map=testimport@user1

+----------------------------------------------------------------------------+
| Layer: testimport@user1 Date: Mon Jan 14 15:45:08 2013
|
| Mapset: user1 Login of Creator: myricaria
|
| Location: lidarprova
|
| DataBase: C:\dl\gl\cavallini
|
| Title: ( testimport )
|
| Timestamp: none
|

|----------------------------------------------------------------------------|
|
|
| Type of Map: raster Number of Categories: 255
|
| Data Type: FCELL
|
| Rows: 1240
|
| Columns: 1640
|
| Total Cells: 2033600
|
| Projection: Transverse Mercator
|
| N: 4896020 S: 4894780 Res: 1
|
| E: 1570420 W: 1568780 Res: 1
|
| Range of data: min = 76.581 max = 214.843
|
|
|
| Data Description:
|
| generated by r.in.gdal
|
|
|
| Comments:
|
| r.in.gdal -o input="C:\dl\gl\cavallini\23B07.asc" output="testimport"
|
|
|

+----------------------------------------------------------------------------+

the extent seems to be the same?

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/Strange-import-of-asc-tp5027341p5027346.html
Sent from the Grass - Dev mailing list archive at Nabble.com.

Il 14/01/2013 15:51, Helmut Kudrnovsky ha scritto:

Hi Paolo,

the extent seems to be the same?

yes - have you tired displaying it under the vector?
additionally, in qgis-dev I get errors querying the raster:
168 bytes expected but 3200000 byte were read from qgis.d.rast
thanks a lot for your help.
--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

yes - have you tired displaying it under the vector?

yes, the newly imported raster and also your raster delivered with the
location are aligned with the vector.

I've also tested a raster link with r.external, this linked is also aligned
to the vector.

HTH

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/Strange-import-of-asc-tp5027341p5027409.html
Sent from the Grass - Dev mailing list archive at Nabble.com.

Il 14/01/2013 20:17, Helmut Kudrnovsky ha scritto:

yes - have you tired displaying it under the vector?

yes, the newly imported raster and also your raster delivered with the
location are aligned with the vector.

I've also tested a raster link with r.external, this linked is also aligned
to the vector.

I get this strange artifact (both on GRASS and on QGIS GUIs), see attached (in red
the vector derived from the imported raster, on the background the imported raster;
the original one is perfectly aligned with the vector). I suppose as a consequence of
this, querying values gives wrong results.
Thanks a lot.
--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

(attachments)

shot2013-01-15 10:24:33.png

Il 14/01/2013 20:17, Helmut Kudrnovsky ha scritto:

yes - have you tired displaying it under the vector?

yes, the newly imported raster and also your raster delivered with the
location are aligned with the vector.

I've also tested a raster link with r.external, this linked is also aligned
to the vector.

at other zoom levels, I get different misaligment, see attached

--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

(attachments)

shot2.png

I get this strange artifact (both on GRASS and on QGIS GUIs)

here on my side it's aligned, see attached screenshot (your vector and an
imported raster of your asc-file).

screenshot.png
<http://osgeo-org.1560.n6.nabble.com/file/n5027529/screenshot.png&gt;

I suppose as a consequence of
this, querying values gives wrong results.

probably.

which grass version are you using? here on my side wingrass6.4.3RC2.

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.n6.nabble.com/Strange-import-of-asc-tp5027341p5027529.html
Sent from the Grass - Dev mailing list archive at Nabble.com.

Hi Paolo

I have imported your layers in grass 6.4 and the vector and both the rasters (original and the asci after import) align perfectly well. I tried the same in grass 7.0, with the same result.

Cheers,

Paulo

ยทยทยท

On 01/15/2013 10:43 AM, Paolo Cavallini wrote:

Il 14/01/2013 20:17, Helmut Kudrnovsky ha scritto:

yes - have you tired displaying it under the vector? 

yes, the newly imported raster and also your raster delivered with the
location are aligned with the vector.

I've also tested a raster link with r.external, this linked is also aligned
to the vector.

at other zoom levels, I get different misaligment, see attached

_______________________________________________
grass-dev mailing list
[grass-dev@lists.osgeo.org](mailto:grass-dev@lists.osgeo.org)
[http://lists.osgeo.org/mailman/listinfo/grass-dev](http://lists.osgeo.org/mailman/listinfo/grass-dev)

Il 15/01/2013 14:42, Paulo van Breugel ha scritto:

Hi Paolo

I have imported your layers in grass 6.4 and the vector and both the rasters
(original and the asci after import) align perfectly well. I tried the same in grass
7.0, with the same result.

Thanks a lot to both of you for checking. I have this issue both with 6.4.2 on Debian
and (I suppose) 6.4.3RC2 from osgeo4w.
Really puzzling. Ideas on wht to test?
All the best.
--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

On Tue, Jan 15, 2013 at 10:27 AM, Paolo Cavallini <cavallini@faunalia.it> wrote:

Il 14/01/2013 20:17, Helmut Kudrnovsky ha scritto:

yes - have you tired displaying it under the vector?

yes, the newly imported raster and also your raster delivered with the
location are aligned with the vector.

I've also tested a raster link with r.external, this linked is also aligned
to the vector.

I get this strange artifact (both on GRASS and on QGIS GUIs)

GRASS GUI? Could you try with pure GRASS d.mon, d.rast, d.vect? The
problem could be the same as in QGIS, lost window precision.

Radim

, see attached (in red
the vector derived from the imported raster, on the background the imported raster;
the original one is perfectly aligned with the vector). I suppose as a consequence of
this, querying values gives wrong results.
Thanks a lot.
--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario

_______________________________________________
grass-dev mailing list
grass-dev@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-dev

Il 18/01/2013 11:55, Radim Blazek ha scritto:

GRASS GUI? Could you try with pure GRASS d.mon, d.rast, d.vect? The
problem could be the same as in QGIS, lost window precision.

It was in fact a QGIS bug, beautifully fixed by Radim.
Thanks, and sorry for the noise.

--
Paolo Cavallini - Faunalia
www.faunalia.eu
Full contact details at www.faunalia.eu/pc
Nuovi corsi QGIS e PostGIS: http://www.faunalia.it/calendario