d.out.file: world (wld) file support issue

classic Classic list List threaded Threaded
9 messages Options
Reply | Threaded
Open this post in threaded view
|

d.out.file: world (wld) file support issue

Markus Neteler
Hi,

I wanted to add world file support to d.out.file but it seems that
d.info -g doesn't work correctly (-g just added in SVN) in case that
the "res=" parameter of d.out.file is used. The generated world file
remains identical and thus the resulting map shifted.

What's wrong here? Patch attached to generate the world file.

thanks
Markus

_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev

d_out_file.diff (1K) Download Attachment
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

Dylan Beaudette-2
On Friday 25 January 2008, Markus Neteler wrote:

> Hi,
>
> I wanted to add world file support to d.out.file but it seems that
> d.info -g doesn't work correctly (-g just added in SVN) in case that
> the "res=" parameter of d.out.file is used. The generated world file
> remains identical and thus the resulting map shifted.
>
> What's wrong here? Patch attached to generate the world file.
>
> thanks
> Markus

I am not sure what is wrong, but... would it make sense to use the facillities
for world file creation that GDAL provides ?  via r.out.gdal ?

maybe I am not seeing the subtle differences..

cheers,

--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

Markus Neteler
On Jan 25, 2008 11:21 PM, Dylan Beaudette <[hidden email]> wrote:

>
> On Friday 25 January 2008, Markus Neteler wrote:
> > Hi,
> >
> > I wanted to add world file support to d.out.file but it seems that
> > d.info -g doesn't work correctly (-g just added in SVN) in case that
> > the "res=" parameter of d.out.file is used. The generated world file
> > remains identical and thus the resulting map shifted.
> >
> > What's wrong here? Patch attached to generate the world file.
> >
> > thanks
> > Markus
>
> I am not sure what is wrong, but... would it make sense to use the facillities
> for world file creation that GDAL provides ?  via r.out.gdal ?
>
> maybe I am not seeing the subtle differences..

I want to write out the *monitor* content with d.out.file (r.out.gdal cannot
do that).

Markus
_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

Dylan Beaudette-2
On Friday 25 January 2008, Markus Neteler wrote:

> On Jan 25, 2008 11:21 PM, Dylan Beaudette <[hidden email]> wrote:
> > On Friday 25 January 2008, Markus Neteler wrote:
> > > Hi,
> > >
> > > I wanted to add world file support to d.out.file but it seems that
> > > d.info -g doesn't work correctly (-g just added in SVN) in case that
> > > the "res=" parameter of d.out.file is used. The generated world file
> > > remains identical and thus the resulting map shifted.
> > >
> > > What's wrong here? Patch attached to generate the world file.
> > >
> > > thanks
> > > Markus
> >
> > I am not sure what is wrong, but... would it make sense to use the
> > facillities for world file creation that GDAL provides ?  via r.out.gdal
> > ?
> >
> > maybe I am not seeing the subtle differences..
>
> I want to write out the *monitor* content with d.out.file (r.out.gdal
> cannot do that).
>
> Markus

right.. (whacks head against table...)

--
Dylan Beaudette
Soil Resource Laboratory
http://casoilresource.lawr.ucdavis.edu/
University of California at Davis
530.754.7341
_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

hamish-2
In reply to this post by Markus Neteler
Markus Neteler wrote:
> I wanted to add world file support to d.out.file but it seems that
> d.info -g doesn't work correctly (-g just added in SVN) in case that
> the "res=" parameter of d.out.file is used. The generated world file
> remains identical and thus the resulting map shifted.
>
> What's wrong here? Patch attached to generate the world file.

Hi,
I haven't looked at the patch yet, but some comments.

first of all why not just use g.region for this? when are display NSEW
and g.region NSEW different? I don't understand the point.

display pixel resolution is resampled. It has nothing to do with region
resolution. so to caclulate it you would need to do
(east-west)/$GRASS_WIDTH   not  (east-west)/ewres

have a look at d.out.gpsdrive. This does the same, it creates world
file(like) output of another flavour, dealing with things like double
and quadruple sizing. [For that I found it nice to quadruple size then
use GIMP's cubic resampler to mix back to target size. Much nicer
output than GRASS's nearest neighbor resampling.]


A solution I have used in the past is to use d.out.file, then
r.in.gdal(XY) + r.region. Then r.out.gdal could make the geotiff of the
image.

ISTR something about using r.mapcalc's r#map, g#map, b#map operators
but don't quite remember how that fit into the puzzle. That was done in
relation to allowing multiple raster overlays in ps.map AFAICR. (??)


other Qs:
how to deal with white bands if region aspect ratio does not match
display window aspect ratio? (see previous changes to d.info needed for
d.out.gpsdrive)

need to check that changes to d.info work correctly in sub-frames.



Hamish



      ____________________________________________________________________________________
Never miss a thing.  Make Yahoo your home page.
http://www.yahoo.com/r/hs

_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

hamish-2
> Markus Neteler wrote:
> > I wanted to add world file support to d.out.file but it seems that
> > d.info -g doesn't work correctly (-g just added in SVN)

Hamish wrote:
> first of all why not just use g.region for this? when are display
> NSEW and g.region NSEW different? I don't understand the point.


oh, NSEW bounds *including* the white out-of-region bands on the side
of the display if the aspect ratios don't match.

anyway, the problem already solved (??) in d.out.gpsdrive in what I
think is a cleaner way. ('d.info -b' + g.region -g,-e)

Looking at it I now remember d.out.gpsdrive tiles are of a fixed size,
I was doubling/quadrupling the display window size & cubic rescaling
with GIMP in a local modification version.


Hamish



      ____________________________________________________________________________________
Never miss a thing.  Make Yahoo your home page.
http://www.yahoo.com/r/hs

_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

Markus Neteler
On Jan 26, 2008 1:02 AM, Hamish <[hidden email]> wrote:

> > Markus Neteler wrote:
> > > I wanted to add world file support to d.out.file but it seems that
> > > d.info -g doesn't work correctly (-g just added in SVN)
>
> Hamish wrote:
> > first of all why not just use g.region for this? when are display
> > NSEW and g.region NSEW different? I don't understand the point.
>
>
> oh, NSEW bounds *including* the white out-of-region bands on the side
> of the display if the aspect ratios don't match.

Exactly. That's why I wanted to get the coordinates from the monitor
itself.

> anyway, the problem already solved (??) in d.out.gpsdrive in what I
> think is a cleaner way. ('d.info -b' + g.region -g,-e)

Thanks, I wasn't aware of d.out.gpsdrive...

> Looking at it I now remember d.out.gpsdrive tiles are of a fixed size,
> I was doubling/quadrupling the display window size & cubic rescaling
> with GIMP in a local modification version.

mhh mhh mhh, mumble. I'll think about it.

Markus
_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

hamish-2
> > Hamish wrote:
> > anyway, the problem already solved (??) in d.out.gpsdrive in what I
> > think is a cleaner way. ('d.info -b' + g.region -g,-e)

Markus:
> Thanks, I wasn't aware of d.out.gpsdrive...
>
> > Looking at it I now remember d.out.gpsdrive tiles are of a fixed
> > size, I was doubling/quadrupling the display window size & cubic
> > rescaling with GIMP in a local modification version.
>
> mhh mhh mhh, mumble. I'll think about it.

the gimp thing was just a local hack to get nicer output when creating
the default GPSDrive backdrop set (Blue Marble).

May I ask what are you trying to include in the GeoTIFF(?)?
Decorations? vector overlays? grid lines? '+' fiducial marks? (ps.map
can make those but d.grid still lacks a flag for that)


Hamish



      ____________________________________________________________________________________
Be a better friend, newshound, and
know-it-all with Yahoo! Mobile.  Try it now.  http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ 


_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev
Reply | Threaded
Open this post in threaded view
|

Re: d.out.file: world (wld) file support issue

hamish-2
Hamish wrote:
> '+' fiducial marks? (ps.map can make those but d.grid still lacks a
> flag for that)

new fiducial mark symbol added in 6.3svn. see attached preivew + detail
and
 http://grass.gdf-hannover.de/wiki/IconSymbols#Existing_symbols

size = 12,15 looks nice in an xmon.
use 45 degree symbol rotation if you wish to align ticks with grid
lines


It would be possible use D_symbol() in d.grid to draw them. (propose
new -c cross & -f fiducial flags for d.grid) Because the way d.grid is
written some new code would be needed to caclulate grid line
interesctions. 45 deg rotation should be applied to align ticks with
grid lines. for geogrids re-use text label rotation angle code to align
ticks with geo grid lines. probably need to calculate that per grid
interection if the map scale is such that the geo-lines are
significantly curved. Maybe for -c just use simple cross icon code from
d.graph, not full D_symbol(cross1).


enjoy,
Hamish



      ____________________________________________________________________________________
Looking for last minute shopping deals?  
Find them fast with Yahoo! Search.  http://tools.search.yahoo.com/newsearch/category.php?category=shopping
_______________________________________________
grass-dev mailing list
[hidden email]
http://lists.osgeo.org/mailman/listinfo/grass-dev

symb-fiducial.png (1K) Download Attachment
fiducial_2k.png (19K) Download Attachment