def save(self, outpath, resolution=75., oversample=2., size=None,
             width=None, height=None):
        gmt = self.gmt
        self.draw_labels()
        self.draw_axes()
        self.draw_cities()
        size_factor = 9.
        if self.show_topo and self.show_topo_scale:
            self._draw_topo_scale()

        if self.stations:
            lats = [s.lat for s in self.stations]
            lons = [s.lon for s in self.stations]

            self.gmt.psxy(
                in_columns=(lons, lats),
                S='t10p',
                G='black',
                *self.jxyr)

            for i_station, s in enumerate(self.stations):
                if self.station_label_mapping:
                    label = self.station_label_mapping[i_station]
                else:
                    label = "%s"%s.station
                self.add_label(s.lat, s.lon, label)
        if self.event:
            e = self.event
            if e.moment_tensor:
                if e.magnitude == None:
                    e.magnitude = e.moment_tensor.magnitude
                size_cm = math.sqrt(math.sqrt(mag2mom(e.magnitude) / 100e17)) * size_factor
                m = e.moment_tensor
                m6 = m.m6_up_south_east()
                data = (e.lon, e.lat, 10, m.strike1, m.dip1, m.rake1, 1, 0, 0, 'M %1.1f' % e.magnitude)
                if True:
                    self.gmt.psmeca(
                        S='%s%g' % ('a', size_cm*2.0),
                        G='red',
                        #W='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        #L='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        in_rows=[data],
                        *self.jxyr)
        gmt = self._gmt
        if outpath.endswith('.eps'):
            tmppath = gmt.tempfilename() + '.eps'
        elif outpath.endswith('.ps'):
            tmppath = gmt.tempfilename() + '.ps'
        else:
            tmppath = gmt.tempfilename() + '.pdf'

        gmt.save(tmppath)

        if any(outpath.endswith(x) for x in ('.eps', '.ps', '.pdf')):
            shutil.copy(tmppath, outpath)
        else:
            convert_graph(tmppath, outpath,
                          resolution=resolution, oversample=oversample,
                          size=size, width=width, height=height)
Пример #2
0
    def save(
        self,
        outpath,
        resolution=75.0,
        oversample=2.0,
        size=None,
        width=None,
        height=None,
    ):
        self.draw_labels()
        self.draw_axes()
        self.draw_cities()
        size_factor = 9.0
        if self.show_topo and self.show_topo_scale:
            self._draw_topo_scale()

        if self.stations:
            lats = [s.lat for s in self.stations]
            lons = [s.lon for s in self.stations]

            self.gmt.psxy(in_columns=(lons, lats), S="t10p", G="black", *self.jxyr)

            for i_station, s in enumerate(self.stations):
                label = self.station_label_mapping.get(i_station, str(s.station))
                self.add_label(s.lat, s.lon, label)

        if self.event:
            e = self.event
            if e.moment_tensor:
                if e.magnitude is None:
                    e.magnitude = e.moment_tensor.magnitude
                size_cm = (
                    math.sqrt(math.sqrt(mag2mom(e.magnitude) / 100e17)) * size_factor
                )
                m = e.moment_tensor
                data = (
                    e.lon,
                    e.lat,
                    10,
                    m.strike1,
                    m.dip1,
                    m.rake1,
                    1,
                    0,
                    0,
                    "M %1.1f" % e.magnitude,
                )
                if True:
                    self.gmt.psmeca(
                        S="%s%g" % ("a", size_cm * 2.0),
                        G="red",
                        # W='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        # L='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        in_rows=[data],
                        *self.jxyr
                    )
        gmt = self._gmt
        if outpath.endswith(".eps"):
            tmppath = gmt.tempfilename() + ".eps"
        elif outpath.endswith(".ps"):
            tmppath = gmt.tempfilename() + ".ps"
        else:
            tmppath = gmt.tempfilename() + ".pdf"

        gmt.save(tmppath)

        if any(outpath.endswith(x) for x in (".eps", ".ps", ".pdf")):
            shutil.copy(tmppath, outpath)
        else:
            convert_graph(
                tmppath,
                outpath,
                resolution=resolution,
                oversample=oversample,
                size=size,
                width=width,
                height=height,
            )
Пример #3
0
def make_map(lat=None,
             lon=None,
             radius=None,
             outfn=None,
             stations=None,
             events=None,
             stations_label_mapping=None,
             map_parameters=None,
             show_topo=True):
    if map_parameters:
        lat = map_parameters.lat
        lon = map_parameters.lon
        radius = map_parameters.radius
        outfn = map_parameters.outfn
        stations = map_parameters.stations
        station_label_mapping = map_parameters.station_label_mapping
        events = map_parameters.events
        station_colors = map_parameters.station_colors
        show_topo = map_parameters.show_topo
        color_wet = map_parameters.color_wet
        color_dry = map_parameters.color_dry
        show_grid = map_parameters.show_grid
        height = map_parameters.height
        width = map_parameters.width

    _map = JMap(width=width,
                height=height,
                lat=lat,
                lon=lon,
                radius=radius,
                topo_resolution_max=200,
                topo_resolution_min=40.,
                show_topo=show_topo,
                show_grid=show_grid,
                illuminate=True,
                illuminate_factor_land=0.5,
                illuminate_factor_ocean=0.25,
                color_wet=color_wet,
                color_dry=color_dry)

    #_map.draw_cities()
    if stations:
        for i_station, s in enumerate(stations):
            lats = [s.lat]
            lons = [s.lon]
            if station_colors:
                color = "%s" % station_colors[s.nsl()]
            else:
                color = "black"

            _map.gmt.psxy(in_columns=(lons, lats),
                          S='t20p',
                          G=color,
                          *_map.jxyr)

            if station_label_mapping:
                label = station_label_mapping[i_station]
            else:
                label = "%s" % s.station
            _map.add_label(s.lat, s.lon, label)
    if events:
        shifts = [(-0.6, -0.6), (-0.6, 0.)]
        for i_e, e in enumerate(events):
            if e.moment_tensor:
                _map.gmt.psxy(in_columns=([e.lon], [e.lat]),
                              S='a25p',
                              G='red',
                              *_map.jxyr)

                break

                size_cm = math.sqrt(
                    math.sqrt(mag2mom(e.moment_tensor.magnitude) /
                              10e7)) * size_factor
                m = e.moment_tensor
                #mc = uniso(m)
                #mc = mc / e.moment_tensor.scalar_moment() * mag2mom(5.0)
                #m6 = to6(c)
                m6 = m.m6_up_south_east()
                #data = (e.lon), (e.lat) + (10,) + m6 + (1, 0, 0)
                #data = lonlat_to_en1_km(e.lon, e.lat) + (10,) + m6 + (1, 0, 0)
                #data = [e.lon, e.lat, 5, m6, 1,0,0]
                eshift, nshift = shifts[i_e]
                #data = (e.lon, e.lat, 10, m.strike1, m.dip1, m.rake1, 1, e.lon+eshift, e.lat+nshift, 'Test site')
                data = (e.lon, e.lat, 10, 1, 1, 1, 0, 0, 0, 1, e.lon, e.lat,
                        'Test site')
                if True:
                    _map.gmt.psmeca(
                        S='%s%g' % ('m', size_cm * 2.0),
                        #G=gmtpy.color(colors[e.cluster]),
                        #G=colors[i_e],
                        G='red',
                        C='3p,0/0/0',
                        #W='thinnest,%i/%i/%i' % (255, 255, 255),
                        #L='thinnest,%i/%i/%i' % (255, 255, 255),
                        in_rows=[data],
                        *_map.jxyr)
    _map.save(outpath=outfn)
def make_map(lat=None, lon=None, radius=None, outfn=None, stations=None, events=None, stations_label_mapping=None, map_parameters=None, show_topo=True):
    if map_parameters:
        lat=map_parameters.lat
        lon=map_parameters.lon
        radius=map_parameters.radius
        outfn = map_parameters.outfn
        stations = map_parameters.stations
        station_label_mapping = map_parameters.station_label_mapping
        events = map_parameters.events
        station_colors = map_parameters.station_colors
        show_topo = map_parameters.show_topo
        color_wet = map_parameters.color_wet
        color_dry = map_parameters.color_dry
        show_grid = map_parameters.show_grid
        height = map_parameters.height
        width = map_parameters.width

    _map = JMap(width=width, height=height, lat=lat, lon=lon,
        radius=radius,
        topo_resolution_max=200,
        topo_resolution_min=40.,
        show_topo=show_topo,
        show_grid=show_grid,
        illuminate=True,
        illuminate_factor_land=0.5,
        illuminate_factor_ocean=0.25,
        color_wet=color_wet,
        color_dry=color_dry)

    #_map.draw_cities()
    if stations:
        for i_station, s in enumerate(stations):
            lats = [s.lat]
            lons = [s.lon]
            if station_colors:
                color = "%s" % station_colors[s.nsl()]
            else:
                color = "black"

            _map.gmt.psxy(
                in_columns=(lons, lats),
                S='t20p',
                G=color,
                *_map.jxyr)

            if station_label_mapping:
                label = station_label_mapping[i_station]
            else:
                label = "%s"%s.station
            _map.add_label(s.lat, s.lon, label)
    if events:
        shifts = [(-0.6, -0.6), (-0.6, 0.)]
        for i_e, e in enumerate(events):
            if e.moment_tensor:
                _map.gmt.psxy(
                    in_columns=([e.lon], [e.lat]),
                    S='a25p',
                    G='red',
                    *_map.jxyr)

                break

                size_cm = math.sqrt(math.sqrt(mag2mom(e.moment_tensor.magnitude) / 10e7)) * size_factor
                m = e.moment_tensor
                #mc = uniso(m)
                #mc = mc / e.moment_tensor.scalar_moment() * mag2mom(5.0)
                #m6 = to6(c)
                m6 = m.m6_up_south_east()
                #data = (e.lon), (e.lat) + (10,) + m6 + (1, 0, 0)
                #data = lonlat_to_en1_km(e.lon, e.lat) + (10,) + m6 + (1, 0, 0)
                #data = [e.lon, e.lat, 5, m6, 1,0,0]
                eshift, nshift = shifts[i_e]
                #data = (e.lon, e.lat, 10, m.strike1, m.dip1, m.rake1, 1, e.lon+eshift, e.lat+nshift, 'Test site')
                data = (e.lon, e.lat, 10, 1,1,1,0,0,0, 1, e.lon, e.lat, 'Test site')
                if True:
                    _map.gmt.psmeca(
                        S='%s%g' % ('m', size_cm*2.0),
                        #G=gmtpy.color(colors[e.cluster]),
                        #G=colors[i_e],
                        G='red',
                        C='3p,0/0/0',
                        #W='thinnest,%i/%i/%i' % (255, 255, 255),
                        #L='thinnest,%i/%i/%i' % (255, 255, 255),
                        in_rows=[data],
                        *_map.jxyr)
    _map.save(outpath=outfn)
Пример #5
0
    def save(self,
             outpath,
             resolution=75.,
             oversample=2.,
             size=None,
             width=None,
             height=None):
        gmt = self.gmt
        self.draw_labels()
        self.draw_axes()
        self.draw_cities()
        size_factor = 9.
        if self.show_topo and self.show_topo_scale:
            self._draw_topo_scale()

        if self.stations:
            lats = [s.lat for s in self.stations]
            lons = [s.lon for s in self.stations]

            self.gmt.psxy(in_columns=(lons, lats),
                          S='t10p',
                          G='black',
                          *self.jxyr)

            for i_station, s in enumerate(self.stations):
                label = self.station_label_mapping.get(i_station,
                                                       str(s.station))
                self.add_label(s.lat, s.lon, label)

        if self.event:
            e = self.event
            if e.moment_tensor:
                if e.magnitude == None:
                    e.magnitude = e.moment_tensor.magnitude
                size_cm = math.sqrt(math.sqrt(
                    mag2mom(e.magnitude) / 100e17)) * size_factor
                m = e.moment_tensor
                m6 = m.m6_up_south_east()
                data = (e.lon, e.lat, 10, m.strike1, m.dip1, m.rake1, 1, 0, 0,
                        'M %1.1f' % e.magnitude)
                if True:
                    self.gmt.psmeca(
                        S='%s%g' % ('a', size_cm * 2.0),
                        G='red',
                        #W='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        #L='thinnest,%i/%i/%i' % darken(gmtpy.color_tup(colors[e.cluster])),
                        in_rows=[data],
                        *self.jxyr)
        gmt = self._gmt
        if outpath.endswith('.eps'):
            tmppath = gmt.tempfilename() + '.eps'
        elif outpath.endswith('.ps'):
            tmppath = gmt.tempfilename() + '.ps'
        else:
            tmppath = gmt.tempfilename() + '.pdf'

        gmt.save(tmppath)

        if any(outpath.endswith(x) for x in ('.eps', '.ps', '.pdf')):
            shutil.copy(tmppath, outpath)
        else:
            convert_graph(tmppath,
                          outpath,
                          resolution=resolution,
                          oversample=oversample,
                          size=size,
                          width=width,
                          height=height)