Bug #667

-JG/g projection (General Perspective)

Added by Chenjian about 5 years ago. Updated over 4 years ago.

Status:In ProgressStart date:2015-01-27
Priority:NormalDue date:
Assignee:Paul% Done:

100%

Category:-
Target version:Candidate for next bugfix release
Affected version:5.2-svn Platform:

Description

Hi,
When I use -JG projection, some projection centers cause weird behaviors of latitudinal gridlines (See attached script and images), e.g.

lon0=188.144925731
lat0=80.7167845195
#lon0=210.651901558 
#lat0=78.8472134393 
#lon0=190.718850211
#lat0=82.2203558311
psbasemap -Rg -JG$lon0/$lat0/1900/0/0/0/170/180/6i -Ba30f10g10/a5f5g5 -P > test.ps

JG.sh Magnifier (202 Bytes) Chenjian, 2015-01-27 11:33

test1.png (101 KB) Chenjian, 2015-01-27 11:33

test2.png (98.2 KB) Chenjian, 2015-01-27 11:33

test3.png (97.1 KB) Chenjian, 2015-01-27 11:33

test5.png (92.8 KB) Chenjian, 2015-02-10 17:20

History

#1 Updated by Paul about 5 years ago

  • Status changed from New to Resolved
  • Assignee set to Paul
  • Target version set to Candidate for next bugfix release
  • % Done changed from 0 to 100

Thanks for reporting this bug. I believe I have fixed it in r14037 as all tests currently run. I added one of your examples as a new test script (tests/genper/polar.sh) in case the problem resurface later...

#2 Updated by Chenjian about 5 years ago

Thank you for your time! Unfortunately, this one:

lon0=280.969616634
lat0=58.7193421224
psbasemap -Rg -JG$lon0/$lat0/2513/0/0/0/200/180/6i -Ba30f30g30/a5f5g5NWSE -P > test.ps

still gives those lines. My projection centers and altitudes change with different point picks in a loop. Most results are good though.

#3 Updated by Paul about 5 years ago

  • Status changed from Resolved to In Progress

I see. Seems our algorithms do not check for points outside this funny octagon map boundary. Will have to consider various strategies here.

#4 Updated by Paul about 5 years ago

  • Status changed from In Progress to Resolved

I've implemented a new function to deal with gridline clipping. In r14056. Also added your other examples which all work for me now.

#5 Updated by Chenjian about 5 years ago

Sorry. It still fails sometimes. Please see this case (test5.png attached):

lon0=225.989251068
lat0=80.3192229138
psbasemap -Rg -JG$lon0/$lat0/6166.78/0/0/54.6281322919/68.2/68.2/6i -Ba30f30g30/a5f5g5NW -P > test.ps

#6 Updated by Chenjian almost 5 years ago

Please excuse me for noise. It still fails at times. Please see the case and its reasult test5.png in the last post.

#7 Updated by Paul almost 5 years ago

  • Status changed from Resolved to In Progress

We are aware of that. I changed the status to in progress - nothing new to report at this time.

#8 Updated by Paul over 4 years ago

No progress except I have turned your test5.png case into a test script so that we are reminded that it is still failing...

Also available in: Atom PDF