re. r43579,
r.out.arc: flag added to refer to cell center in header
why is it necessary or a good thing to let the user export data in the
incorrect format? I presume that Arc's import tool is very specific in
it's definition of the header, and that we are already shipping the
correct way.
ie what's the justification for violating the read-sloppy, write-exact
rule-of-thumb?
curious,
Hamish
On Sun, Oct 17, 2010 at 8:43 AM, Hamish <hamish_b@yahoo.com> wrote:
re. r43579,
r.out.arc: flag added to refer to cell center in headerwhy is it necessary or a good thing to let the user export data in the
incorrect format?
Why incorrect?
I presume that Arc's import tool is very specific in
it's definition of the header, and that we are already shipping the
correct way.
There are two notations.
ie what's the justification for violating the read-sloppy, write-exact
rule-of-thumb?
Please explain what you mean... I don't see any violation yet nor I
would support it. If you look at the change:
153 if (G_projection() != PROJECTION_LL) { /* Is Projection
!= LL (3) */
154 if (!flag.ccenter->answer) {
155 G_format_easting(region.west, buf, region.proj);
X 156 fprintf(fp, "xllcorner %s\n", buf);
157 G_format_northing(region.south, buf, region.proj);
X 158 fprintf(fp, "yllcorner %s\n", buf);
159 }
160 else {
161 G_format_easting(region.west + cellsize / 2.,
buf, region.proj);
X 162 fprintf(fp, "xllcenter %s\n", buf);
163 G_format_northing(region.south + cellsize / 2.,
buf, region.proj);
X 164 fprintf(fp, "yllcenter %s\n", buf);
165 }
X--> See that the reference changes, too.
Some software packages do not understand the xllcorner/yllcorner notation
but only the xllcenter/yllcenter notation.
Markus