def plot(): ps = PlotSphere(38, 127.5, figsize=(15,15)) # Korea centered ps.draw_cube_panel(2) #ps.draw_cube_panel(3) ps.draw_cube_panel(6) #latlon1 = xyp2latlon(0.5, 0.4, 2) #latlon2 = xyp2latlon(0.4, 0.5, 6) #latlon1 = xyp2latlon(at, 0.4, 2) #latlon2 = xyp2latlon(0.4, at, 6) #latlon1 = xyp2latlon(0.4, 0.3, 2) #latlon2 = xyp2latlon(-0.4, 0.2, 3) latlon1 = xyp2latlon(0.5, 0.4, 2) latlon2 = xyp2latlon(0.4, 0.5, 6) latlon3 = xyp2latlon(0.35, 0.5, 6) latlon4 = xyp2latlon(0.5, 0.35, 2) #latlon1 = xyp2latlon(0.5, 0.5, 2) #latlon2 = xyp2latlon(0.35, 0.5, 6) #latlon3 = xyp2latlon(at, at, 6) #latlon4 = xyp2latlon(0.5, 0.35, 2) #draw_greatcircle_line(ps.bmap, latlon1, latlon2, c='r') #draw_greatcircle_line(ps.bmap, latlon3, latlon4, c='b') latlons = [latlon1, latlon4, latlon3, latlon2] poly = draw_polygon(ps.bmap, latlons) poly.update(dict(fc='r')) draw_points(ps.bmap, latlons) ps.show(True)
def draw_cube_panel(self, panel, **kargs): kwds = dict(c='0.5', lw=1., ls='-') kwds.update(kargs) for side, (corner1,corner2) in panel_side2corners.items(): xy1, xy2 = panel_corner_xys[corner1], panel_corner_xys[corner2] latlon1 = xyp2latlon(xy1[0], xy1[1], panel) latlon2 = xyp2latlon(xy2[0], xy2[1], panel) draw_line(self.bmap, latlon1, latlon2, **kwds)