U xbv!@sddlZddlZddlmZddlmZddl m Z dddddgdde ddfdd Z dddd d dd gfd dZddddgdddfddZddZdS)N)pyplotcrs)tickerTc  Cst|||d\} } }} |r.td|d||jdd|dddd|d|d|d |d d |r|jd d |dd|jdd |ddtdd| } tjd| dd} tjd| dd}|j| |j ||j | |d|j | |d|S)aEAuto draw lat/lon label/tick based on coverage from geo_box Parameters: ax : cartopy axes. geo_box : 4-tuple of float, (W, N, E, S) in degree lalo_step : float, major tick interval for X/Y axes lalo_loc : list of 4 bool, positions where the labels are drawn as in [left, right, top, bottom] default: [1,0,0,1] lalo_max_num : int, maximum number of major ticks for X/Y axes lalo_offset : float, distance in points between tick and label. Set to negative value to move the ticklabel inside the plot. projection : cartopy.crs object ... Example: geo_box = (128.0, 37.0, 138.0, 30.0) m.draw_lalo_label(geo_box) ) lalo_step lalo_max_numzplot lat/lon label in step of z and location of bothinTrr) which direction labelsizeleftrighttopbottom labelleft labelrightlabeltop labelbottomxmajor)axisrpady.f) number_formatr) auto_lalo_sequenceprint tick_paramsmaxctickerZLongitudeFormatterZLatitudeFormatterxaxisset_major_formatteryaxis set_xticks set_yticks)axgeo_boxrlalo_locr lalo_offset projection font_size print_msglatslonsdigitZ dec_digitZ lon_formatterZ lat_formatterr5G/home/exouser/operations/rsmas_insar/sources/MintPy/mintpy/utils/map.pydraw_lalo_labels6   r7r r c st|d|d|d|dg|st}ttt|t|ddkrld7t|}fdd|D}fd d|D}||t |}ttt|t |dd dd d}t |dd dd d}t t ||d | t ||d |f}t t ||d | t ||d |f} t t|}t t| } |tt||dk||dk}| tt| |dk| |dk} || |fS) aAuto calculate lat/lon label sequence based on input geo_box Parameters: geo_box : 4-tuple of float, defining UL_lon, UL_lat, LR_lon, LR_lat coordinate lalo_step : float lalo_max_num : int, rough major tick number along the longer axis step_candidate : list of int, candidate list for the significant number of step Returns: lats/lons : np.array of float, sequence of lat/lon auto calculated from input geo_box lalo_step : float, lat/lon label step Example: geo_box = (128.0, 37.0, 138.0, 30.0) lats, lons, step = m.auto_lalo_sequence(geo_box) rr r r)89csg|]}|dqS) r5.0i)r4r5r6 _sz&auto_lalo_sequence..csg|]}|dqS)r r5r=)r max_lalo_distr5r6r@`sr<g$@)r$ round_to_1npintfloorlog10strroundindexminceilhstackarangesortuniquewhere logical_and) r,rr Zstep_candidateZlalo_step_candidatedistanceZ lat_majorZ lon_majorr2r3r5)r4r rAr6r!Fs0 $  &&&&r!degreesg?皙?g? kc Cstjdd}|r|nt}|d|d|d|d}|d|d|d|d} |drx|d|d} nF|dr||d|d|d|dd} |d|dd kr|St| |d} | d krt| d d } |dr| || d | d} t || } n|dr(| } || d }|| d }|j ||g| | g|d |j ||g| | d| g|d |j ||g| | d| g|d d}| d krd}| d9} |d|d|}|j || d | || dd|dd||d|S)a8draw a simple map scale from x1,y to x2,y in map projection coordinates, label it with actual distance ref_link: http://matplotlib.1069221.n5.nabble.com/basemap-scalebar-td14133.html Parameters: ax : matplotlib.pyplot.axes object geo_box : tuple of 4 float in (x0, y0, x1, y1) for (W, N, E, S) in degrees / meters unit : str, coordinate unit - degrees or meters loc : list of 3 float in (length, y, x) of scale bar location: length = ratio of the total width y / x = axes fraction of the scale bar center labelpad : float Returns: ax : matplotlib.pyplot.axes object Example: from mintpy.utils import plot as pp pp.draw_scale_bar(ax, geo_box) WGS84)ellpsrrr r meterdegg@@Zg@)colorrTmkmgMbP?z.0f center)rrsvahafontsizer]) pyprojGeodpltgca startswithinvrBrCrintfwdabsplottext)r+r,unitloclabelpadr0r]geodlon_clat_cZ scene_widthZ length_meterZlon_c2Z length_displon0lon1Z txt_offsetr5r5r6 draw_scalebarvsT            rycCs$tttt|}t|| S)z1Return the most significant digit of input number)rDrCrErFrnrH)rr4r5r5r6rBsrB)rfnumpyrC matplotlibrrhcartopyrccrsZ cartopy.mplrr% PlateCarreer7r!ryrBr5r5r5r6 s    10L