disdrodb 0.0.21__py3-none-any.whl → 0.1.0__py3-none-any.whl
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- disdrodb/__init__.py +132 -15
- disdrodb/_config.py +4 -2
- disdrodb/_version.py +9 -4
- disdrodb/api/checks.py +264 -237
- disdrodb/api/configs.py +4 -8
- disdrodb/api/create_directories.py +235 -290
- disdrodb/api/info.py +217 -26
- disdrodb/api/io.py +295 -269
- disdrodb/api/path.py +597 -173
- disdrodb/api/search.py +486 -0
- disdrodb/{metadata/scripts → cli}/disdrodb_check_metadata_archive.py +12 -7
- disdrodb/{utils/pandas.py → cli/disdrodb_data_archive_directory.py} +9 -18
- disdrodb/cli/disdrodb_download_archive.py +86 -0
- disdrodb/cli/disdrodb_download_metadata_archive.py +53 -0
- disdrodb/cli/disdrodb_download_station.py +84 -0
- disdrodb/{api/scripts → cli}/disdrodb_initialize_station.py +22 -10
- disdrodb/cli/disdrodb_metadata_archive_directory.py +32 -0
- disdrodb/{data_transfer/scripts/disdrodb_download_station.py → cli/disdrodb_open_data_archive.py} +22 -22
- disdrodb/cli/disdrodb_open_logs_directory.py +69 -0
- disdrodb/{data_transfer/scripts/disdrodb_upload_station.py → cli/disdrodb_open_metadata_archive.py} +22 -24
- disdrodb/cli/disdrodb_open_metadata_directory.py +71 -0
- disdrodb/cli/disdrodb_open_product_directory.py +74 -0
- disdrodb/cli/disdrodb_open_readers_directory.py +32 -0
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0.py +38 -31
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0_station.py +32 -30
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0a.py +30 -21
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0a_station.py +24 -33
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0b.py +30 -21
- disdrodb/{l0/scripts → cli}/disdrodb_run_l0b_station.py +25 -34
- disdrodb/cli/disdrodb_run_l0c.py +130 -0
- disdrodb/cli/disdrodb_run_l0c_station.py +129 -0
- disdrodb/cli/disdrodb_run_l1.py +122 -0
- disdrodb/cli/disdrodb_run_l1_station.py +121 -0
- disdrodb/cli/disdrodb_run_l2e.py +122 -0
- disdrodb/cli/disdrodb_run_l2e_station.py +122 -0
- disdrodb/cli/disdrodb_run_l2m.py +122 -0
- disdrodb/cli/disdrodb_run_l2m_station.py +122 -0
- disdrodb/cli/disdrodb_upload_archive.py +105 -0
- disdrodb/cli/disdrodb_upload_station.py +98 -0
- disdrodb/configs.py +90 -25
- disdrodb/data_transfer/__init__.py +22 -0
- disdrodb/data_transfer/download_data.py +87 -90
- disdrodb/data_transfer/upload_data.py +64 -37
- disdrodb/data_transfer/zenodo.py +15 -18
- disdrodb/docs.py +1 -1
- disdrodb/issue/__init__.py +17 -4
- disdrodb/issue/checks.py +10 -23
- disdrodb/issue/reader.py +9 -12
- disdrodb/issue/writer.py +14 -17
- disdrodb/l0/__init__.py +17 -26
- disdrodb/l0/check_configs.py +35 -23
- disdrodb/l0/check_standards.py +32 -42
- disdrodb/l0/configs/{Thies_LPM → LPM}/bins_diameter.yml +44 -44
- disdrodb/l0/configs/{Thies_LPM → LPM}/bins_velocity.yml +40 -40
- disdrodb/l0/configs/LPM/l0a_encodings.yml +80 -0
- disdrodb/l0/configs/{Thies_LPM → LPM}/l0b_cf_attrs.yml +62 -59
- disdrodb/l0/configs/{Thies_LPM → LPM}/l0b_encodings.yml +9 -9
- disdrodb/l0/configs/{Thies_LPM → LPM}/raw_data_format.yml +245 -245
- disdrodb/l0/configs/{OTT_Parsivel → PARSIVEL}/bins_diameter.yml +66 -66
- disdrodb/l0/configs/{OTT_Parsivel → PARSIVEL}/bins_velocity.yml +64 -64
- disdrodb/l0/configs/PARSIVEL/l0a_encodings.yml +32 -0
- disdrodb/l0/configs/{OTT_Parsivel → PARSIVEL}/l0b_cf_attrs.yml +22 -20
- disdrodb/l0/configs/{OTT_Parsivel → PARSIVEL}/l0b_encodings.yml +17 -17
- disdrodb/l0/configs/{OTT_Parsivel → PARSIVEL}/raw_data_format.yml +77 -77
- disdrodb/l0/configs/{OTT_Parsivel2 → PARSIVEL2}/bins_diameter.yml +64 -64
- disdrodb/l0/configs/{OTT_Parsivel2 → PARSIVEL2}/bins_velocity.yml +64 -64
- disdrodb/l0/configs/PARSIVEL2/l0a_encodings.yml +39 -0
- disdrodb/l0/configs/{OTT_Parsivel2 → PARSIVEL2}/l0b_cf_attrs.yml +24 -22
- disdrodb/l0/configs/{OTT_Parsivel2 → PARSIVEL2}/l0b_encodings.yml +20 -20
- disdrodb/l0/configs/{OTT_Parsivel2 → PARSIVEL2}/raw_data_format.yml +98 -98
- disdrodb/l0/configs/{RD_80 → RD80}/bins_diameter.yml +40 -40
- disdrodb/l0/configs/RD80/l0a_encodings.yml +16 -0
- disdrodb/l0/configs/{RD_80 → RD80}/l0b_cf_attrs.yml +3 -3
- disdrodb/l0/configs/RD80/l0b_encodings.yml +135 -0
- disdrodb/l0/configs/{RD_80 → RD80}/raw_data_format.yml +48 -48
- disdrodb/l0/l0_reader.py +216 -340
- disdrodb/l0/l0a_processing.py +237 -208
- disdrodb/l0/l0b_nc_processing.py +227 -80
- disdrodb/l0/l0b_processing.py +93 -173
- disdrodb/l0/l0c_processing.py +627 -0
- disdrodb/l0/readers/{ARM → LPM/ARM}/ARM_LPM.py +36 -58
- disdrodb/l0/readers/LPM/AUSTRALIA/MELBOURNE_2007_LPM.py +226 -0
- disdrodb/l0/readers/LPM/BRAZIL/CHUVA_LPM.py +185 -0
- disdrodb/l0/readers/LPM/BRAZIL/GOAMAZON_LPM.py +183 -0
- disdrodb/l0/readers/LPM/ITALY/GID_LPM.py +179 -0
- disdrodb/l0/readers/{UK → LPM/UK}/DIVEN.py +14 -35
- disdrodb/l0/readers/PARSIVEL/AUSTRALIA/MELBOURNE_2007_PARSIVEL.py +157 -0
- disdrodb/l0/readers/PARSIVEL/CHINA/CHONGQING.py +113 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/ARCTIC_2021.py +40 -57
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/COMMON_2011.py +37 -54
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/DAVOS_2009_2011.py +34 -51
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/EPFL_2009.py +34 -51
- disdrodb/l0/readers/{EPFL/PARADISO_2014.py → PARSIVEL/EPFL/EPFL_ROOF_2008.py} +38 -50
- disdrodb/l0/readers/PARSIVEL/EPFL/EPFL_ROOF_2010.py +105 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/EPFL_ROOF_2011.py +34 -51
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/EPFL_ROOF_2012.py +33 -51
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/GENEPI_2007.py +25 -44
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/GRAND_ST_BERNARD_2007.py +25 -44
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/GRAND_ST_BERNARD_2007_2.py +25 -44
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/HPICONET_2010.py +34 -51
- disdrodb/l0/readers/{EPFL/EPFL_ROOF_2010.py → PARSIVEL/EPFL/HYMEX_LTE_SOP2.py} +37 -50
- disdrodb/l0/readers/PARSIVEL/EPFL/HYMEX_LTE_SOP3.py +111 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/HYMEX_LTE_SOP4.py +36 -54
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/LOCARNO_2018.py +34 -52
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/LOCARNO_2019.py +38 -56
- disdrodb/l0/readers/PARSIVEL/EPFL/PARADISO_2014.py +105 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/PARSIVEL_2007.py +27 -45
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/PLATO_2019.py +24 -44
- disdrodb/l0/readers/PARSIVEL/EPFL/RACLETS_2019.py +140 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/RACLETS_2019_WJF.py +41 -59
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/RIETHOLZBACH_2011.py +34 -51
- disdrodb/l0/readers/PARSIVEL/EPFL/SAMOYLOV_2017.py +117 -0
- disdrodb/l0/readers/PARSIVEL/EPFL/SAMOYLOV_2019.py +137 -0
- disdrodb/l0/readers/{EPFL → PARSIVEL/EPFL}/UNIL_2022.py +42 -55
- disdrodb/l0/readers/PARSIVEL/GPM/IFLOODS.py +104 -0
- disdrodb/l0/readers/{GPM → PARSIVEL/GPM}/LPVEX.py +29 -48
- disdrodb/l0/readers/PARSIVEL/GPM/MC3E.py +184 -0
- disdrodb/l0/readers/PARSIVEL/NCAR/CCOPE_2015.py +113 -0
- disdrodb/l0/readers/{NCAR/VORTEX_SE_2016_P1.py → PARSIVEL/NCAR/OWLES_MIPS.py} +46 -72
- disdrodb/l0/readers/PARSIVEL/NCAR/PECAN_MOBILE.py +125 -0
- disdrodb/l0/readers/{NCAR/OWLES_MIPS.py → PARSIVEL/NCAR/PLOWS_MIPS.py} +45 -64
- disdrodb/l0/readers/PARSIVEL/NCAR/VORTEX2_2009.py +114 -0
- disdrodb/l0/readers/PARSIVEL/NCAR/VORTEX2_2010.py +176 -0
- disdrodb/l0/readers/PARSIVEL/NCAR/VORTEX2_2010_UF.py +183 -0
- disdrodb/l0/readers/{ARM/ARM_LD.py → PARSIVEL2/ARM/ARM_PARSIVEL2.py} +27 -50
- disdrodb/l0/readers/PARSIVEL2/BRAZIL/CHUVA_PARSIVEL2.py +163 -0
- disdrodb/l0/readers/PARSIVEL2/BRAZIL/GOAMAZON_PARSIVEL2.py +163 -0
- disdrodb/l0/readers/{DENMARK → PARSIVEL2/DENMARK}/EROSION_nc.py +14 -35
- disdrodb/l0/readers/PARSIVEL2/FRANCE/SIRTA_PARSIVEL2.py +119 -0
- disdrodb/l0/readers/PARSIVEL2/GPM/GCPEX.py +104 -0
- disdrodb/l0/readers/PARSIVEL2/GPM/NSSTC.py +176 -0
- disdrodb/l0/readers/PARSIVEL2/ITALY/GID_PARSIVEL2.py +32 -0
- disdrodb/l0/readers/PARSIVEL2/MEXICO/OH_IIUNAM_nc.py +56 -0
- disdrodb/l0/readers/PARSIVEL2/NCAR/PECAN_FP3.py +120 -0
- disdrodb/l0/readers/{NCAR → PARSIVEL2/NCAR}/PECAN_MIPS.py +45 -64
- disdrodb/l0/readers/PARSIVEL2/NCAR/RELAMPAGO_PARSIVEL2.py +181 -0
- disdrodb/l0/readers/PARSIVEL2/NCAR/SNOWIE_PJ.py +160 -0
- disdrodb/l0/readers/PARSIVEL2/NCAR/SNOWIE_SB.py +160 -0
- disdrodb/l0/readers/{NCAR/PLOWS_MIPS.py → PARSIVEL2/NCAR/VORTEX_SE_2016_P1.py} +49 -66
- disdrodb/l0/readers/PARSIVEL2/NCAR/VORTEX_SE_2016_P2.py +118 -0
- disdrodb/l0/readers/PARSIVEL2/NCAR/VORTEX_SE_2016_PIPS.py +152 -0
- disdrodb/l0/readers/PARSIVEL2/NETHERLANDS/DELFT.py +166 -0
- disdrodb/l0/readers/{NCAR/RELAMPAGO_RD80.py → RD80/BRAZIL/CHUVA_RD80.py} +36 -60
- disdrodb/l0/readers/{BRAZIL → RD80/BRAZIL}/GOAMAZON_RD80.py +36 -55
- disdrodb/l0/readers/{NCAR → RD80/NCAR}/CINDY_2011_RD80.py +35 -54
- disdrodb/l0/readers/{BRAZIL/CHUVA_RD80.py → RD80/NCAR/RELAMPAGO_RD80.py} +40 -54
- disdrodb/l0/readers/template_reader_raw_netcdf_data.py +62 -0
- disdrodb/l0/readers/{reader_template.py → template_reader_raw_text_data.py} +20 -44
- disdrodb/l0/routines.py +885 -581
- disdrodb/l0/standards.py +72 -236
- disdrodb/l0/template_tools.py +104 -109
- disdrodb/l1/__init__.py +17 -0
- disdrodb/l1/beard_model.py +716 -0
- disdrodb/l1/encoding_attrs.py +620 -0
- disdrodb/l1/fall_velocity.py +260 -0
- disdrodb/l1/filters.py +192 -0
- disdrodb/l1/processing.py +200 -0
- disdrodb/l1/resampling.py +236 -0
- disdrodb/l1/routines.py +357 -0
- disdrodb/l1_env/__init__.py +17 -0
- disdrodb/l1_env/routines.py +38 -0
- disdrodb/l2/__init__.py +17 -0
- disdrodb/l2/empirical_dsd.py +1735 -0
- disdrodb/l2/event.py +388 -0
- disdrodb/l2/processing.py +519 -0
- disdrodb/l2/processing_options.py +213 -0
- disdrodb/l2/routines.py +868 -0
- disdrodb/metadata/__init__.py +9 -2
- disdrodb/metadata/checks.py +165 -118
- disdrodb/metadata/download.py +81 -0
- disdrodb/metadata/geolocation.py +146 -0
- disdrodb/metadata/info.py +20 -13
- disdrodb/metadata/manipulation.py +1 -1
- disdrodb/metadata/reader.py +59 -8
- disdrodb/metadata/search.py +77 -144
- disdrodb/metadata/standards.py +7 -8
- disdrodb/metadata/writer.py +8 -14
- disdrodb/psd/__init__.py +38 -0
- disdrodb/psd/fitting.py +2146 -0
- disdrodb/psd/models.py +774 -0
- disdrodb/routines.py +1176 -0
- disdrodb/scattering/__init__.py +28 -0
- disdrodb/scattering/axis_ratio.py +344 -0
- disdrodb/scattering/routines.py +456 -0
- disdrodb/utils/__init__.py +17 -0
- disdrodb/utils/attrs.py +208 -0
- disdrodb/utils/cli.py +269 -0
- disdrodb/utils/compression.py +60 -42
- disdrodb/utils/dask.py +62 -0
- disdrodb/utils/decorators.py +110 -0
- disdrodb/utils/directories.py +107 -46
- disdrodb/utils/encoding.py +127 -0
- disdrodb/utils/list.py +29 -0
- disdrodb/utils/logger.py +168 -46
- disdrodb/utils/time.py +657 -0
- disdrodb/utils/warnings.py +30 -0
- disdrodb/utils/writer.py +57 -0
- disdrodb/utils/xarray.py +138 -47
- disdrodb/utils/yaml.py +0 -1
- disdrodb/viz/__init__.py +17 -0
- disdrodb/viz/plots.py +17 -0
- disdrodb-0.1.0.dist-info/METADATA +321 -0
- disdrodb-0.1.0.dist-info/RECORD +216 -0
- {disdrodb-0.0.21.dist-info → disdrodb-0.1.0.dist-info}/WHEEL +1 -1
- disdrodb-0.1.0.dist-info/entry_points.txt +30 -0
- disdrodb/data_transfer/scripts/disdrodb_download_archive.py +0 -53
- disdrodb/data_transfer/scripts/disdrodb_upload_archive.py +0 -57
- disdrodb/l0/configs/OTT_Parsivel/l0a_encodings.yml +0 -32
- disdrodb/l0/configs/OTT_Parsivel2/l0a_encodings.yml +0 -39
- disdrodb/l0/configs/RD_80/l0a_encodings.yml +0 -16
- disdrodb/l0/configs/RD_80/l0b_encodings.yml +0 -135
- disdrodb/l0/configs/Thies_LPM/l0a_encodings.yml +0 -80
- disdrodb/l0/io.py +0 -257
- disdrodb/l0/l0_processing.py +0 -1091
- disdrodb/l0/readers/AUSTRALIA/MELBOURNE_2007_OTT.py +0 -178
- disdrodb/l0/readers/AUSTRALIA/MELBOURNE_2007_THIES.py +0 -247
- disdrodb/l0/readers/BRAZIL/CHUVA_LPM.py +0 -204
- disdrodb/l0/readers/BRAZIL/CHUVA_OTT.py +0 -183
- disdrodb/l0/readers/BRAZIL/GOAMAZON_LPM.py +0 -204
- disdrodb/l0/readers/BRAZIL/GOAMAZON_OTT.py +0 -183
- disdrodb/l0/readers/CHINA/CHONGQING.py +0 -131
- disdrodb/l0/readers/EPFL/EPFL_ROOF_2008.py +0 -128
- disdrodb/l0/readers/EPFL/HYMEX_LTE_SOP2.py +0 -127
- disdrodb/l0/readers/EPFL/HYMEX_LTE_SOP3.py +0 -129
- disdrodb/l0/readers/EPFL/RACLETS_2019.py +0 -158
- disdrodb/l0/readers/EPFL/SAMOYLOV_2017.py +0 -136
- disdrodb/l0/readers/EPFL/SAMOYLOV_2019.py +0 -158
- disdrodb/l0/readers/FRANCE/SIRTA_OTT2.py +0 -138
- disdrodb/l0/readers/GPM/GCPEX.py +0 -123
- disdrodb/l0/readers/GPM/IFLOODS.py +0 -123
- disdrodb/l0/readers/GPM/MC3E.py +0 -123
- disdrodb/l0/readers/GPM/NSSTC.py +0 -164
- disdrodb/l0/readers/ITALY/GID.py +0 -199
- disdrodb/l0/readers/MEXICO/OH_IIUNAM_nc.py +0 -92
- disdrodb/l0/readers/NCAR/CCOPE_2015.py +0 -133
- disdrodb/l0/readers/NCAR/PECAN_FP3.py +0 -137
- disdrodb/l0/readers/NCAR/PECAN_MOBILE.py +0 -144
- disdrodb/l0/readers/NCAR/RELAMPAGO_OTT.py +0 -195
- disdrodb/l0/readers/NCAR/SNOWIE_PJ.py +0 -172
- disdrodb/l0/readers/NCAR/SNOWIE_SB.py +0 -179
- disdrodb/l0/readers/NCAR/VORTEX2_2009.py +0 -133
- disdrodb/l0/readers/NCAR/VORTEX2_2010.py +0 -188
- disdrodb/l0/readers/NCAR/VORTEX2_2010_UF.py +0 -191
- disdrodb/l0/readers/NCAR/VORTEX_SE_2016_P2.py +0 -135
- disdrodb/l0/readers/NCAR/VORTEX_SE_2016_PIPS.py +0 -170
- disdrodb/l0/readers/NETHERLANDS/DELFT.py +0 -187
- disdrodb/l0/readers/SPAIN/SBEGUERIA.py +0 -179
- disdrodb/l0/scripts/disdrodb_run_l0b_concat.py +0 -93
- disdrodb/l0/scripts/disdrodb_run_l0b_concat_station.py +0 -85
- disdrodb/utils/netcdf.py +0 -452
- disdrodb/utils/scripts.py +0 -102
- disdrodb-0.0.21.dist-info/AUTHORS.md +0 -18
- disdrodb-0.0.21.dist-info/METADATA +0 -186
- disdrodb-0.0.21.dist-info/RECORD +0 -168
- disdrodb-0.0.21.dist-info/entry_points.txt +0 -15
- /disdrodb/l0/configs/{RD_80 → RD80}/bins_velocity.yml +0 -0
- /disdrodb/l0/manuals/{Thies_LPM.pdf → LPM.pdf} +0 -0
- /disdrodb/l0/manuals/{ODM_470.pdf → ODM470.pdf} +0 -0
- /disdrodb/l0/manuals/{OTT_Parsivel.pdf → PARSIVEL.pdf} +0 -0
- /disdrodb/l0/manuals/{OTT_Parsivel2.pdf → PARSIVEL2.pdf} +0 -0
- /disdrodb/l0/manuals/{PWS_100.pdf → PWS100.pdf} +0 -0
- /disdrodb/l0/manuals/{RD_80.pdf → RD80.pdf} +0 -0
- {disdrodb-0.0.21.dist-info → disdrodb-0.1.0.dist-info/licenses}/LICENSE +0 -0
- {disdrodb-0.0.21.dist-info → disdrodb-0.1.0.dist-info}/top_level.txt +0 -0
|
@@ -0,0 +1,140 @@
|
|
|
1
|
+
#!/usr/bin/env python3
|
|
2
|
+
# -----------------------------------------------------------------------------.
|
|
3
|
+
# Copyright (c) 2021-2023 DISDRODB developers
|
|
4
|
+
#
|
|
5
|
+
# This program is free software: you can redistribute it and/or modify
|
|
6
|
+
# it under the terms of the GNU General Public License as published by
|
|
7
|
+
# the Free Software Foundation, either version 3 of the License, or
|
|
8
|
+
# (at your option) any later version.
|
|
9
|
+
#
|
|
10
|
+
# This program is distributed in the hope that it will be useful,
|
|
11
|
+
# but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
12
|
+
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
13
|
+
# GNU General Public License for more details.
|
|
14
|
+
#
|
|
15
|
+
# You should have received a copy of the GNU General Public License
|
|
16
|
+
# along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
17
|
+
# -----------------------------------------------------------------------------.
|
|
18
|
+
"""Reader for RACLETS 2019 campaign."""
|
|
19
|
+
import pandas as pd
|
|
20
|
+
|
|
21
|
+
from disdrodb.l0.l0_reader import is_documented_by, reader_generic_docstring
|
|
22
|
+
from disdrodb.l0.l0a_processing import read_raw_text_file
|
|
23
|
+
|
|
24
|
+
|
|
25
|
+
@is_documented_by(reader_generic_docstring)
|
|
26
|
+
def reader(
|
|
27
|
+
filepath,
|
|
28
|
+
logger=None,
|
|
29
|
+
):
|
|
30
|
+
"""Reader."""
|
|
31
|
+
##------------------------------------------------------------------------.
|
|
32
|
+
#### Define column names
|
|
33
|
+
column_names = ["TO_BE_SPLITTED"]
|
|
34
|
+
|
|
35
|
+
##------------------------------------------------------------------------.
|
|
36
|
+
#### Define reader options
|
|
37
|
+
reader_kwargs = {}
|
|
38
|
+
# - Define delimiter
|
|
39
|
+
reader_kwargs["delimiter"] = ";" # Used to not split anything !
|
|
40
|
+
# - Define encoding
|
|
41
|
+
reader_kwargs["encoding"] = "ISO-8859-1"
|
|
42
|
+
# Skip first row as columns names
|
|
43
|
+
reader_kwargs["header"] = None
|
|
44
|
+
# - Avoid first column to become df index !!!
|
|
45
|
+
reader_kwargs["index_col"] = False
|
|
46
|
+
# - Define behaviour when encountering bad lines
|
|
47
|
+
reader_kwargs["on_bad_lines"] = "skip"
|
|
48
|
+
# - Define reader engine
|
|
49
|
+
# - C engine is faster
|
|
50
|
+
# - Python engine is more feature-complete
|
|
51
|
+
reader_kwargs["engine"] = "python"
|
|
52
|
+
# - Define on-the-fly decompression of on-disk data
|
|
53
|
+
# - Available: gzip, bz2, zip
|
|
54
|
+
reader_kwargs["compression"] = "gzip"
|
|
55
|
+
# - Strings to recognize as NA/NaN and replace with standard NA flags
|
|
56
|
+
# - Already included: '#N/A', '#N/A N/A', '#NA', '-1.#IND', '-1.#QNAN',
|
|
57
|
+
# '-NaN', '-nan', '1.#IND', '1.#QNAN', '<NA>', 'N/A',
|
|
58
|
+
# 'NA', 'NULL', 'NaN', 'n/a', 'nan', 'null'
|
|
59
|
+
reader_kwargs["na_values"] = [
|
|
60
|
+
"na",
|
|
61
|
+
"",
|
|
62
|
+
"error",
|
|
63
|
+
]
|
|
64
|
+
|
|
65
|
+
##------------------------------------------------------------------------.
|
|
66
|
+
#### Read the data
|
|
67
|
+
df = read_raw_text_file(
|
|
68
|
+
filepath=filepath,
|
|
69
|
+
column_names=column_names,
|
|
70
|
+
reader_kwargs=reader_kwargs,
|
|
71
|
+
logger=logger,
|
|
72
|
+
)
|
|
73
|
+
|
|
74
|
+
##------------------------------------------------------------------------.
|
|
75
|
+
#### Adapt the dataframe to adhere to DISDRODB L0 standards
|
|
76
|
+
# Drop row that contains errors
|
|
77
|
+
df = df[~df["TO_BE_SPLITTED"].str.contains("Error in data reading! 0")]
|
|
78
|
+
|
|
79
|
+
# Check if file empty
|
|
80
|
+
if len(df.index) == 0:
|
|
81
|
+
raise ValueError("Error in all rows. The file has been skipped.")
|
|
82
|
+
|
|
83
|
+
# Split the column
|
|
84
|
+
df = df["TO_BE_SPLITTED"].str.split(",", expand=True, n=1111)
|
|
85
|
+
|
|
86
|
+
# Define auxiliary columns
|
|
87
|
+
column_names = [
|
|
88
|
+
"id",
|
|
89
|
+
"latitude",
|
|
90
|
+
"longitude",
|
|
91
|
+
"time",
|
|
92
|
+
"datalogger_temperature",
|
|
93
|
+
"datalogger_voltage",
|
|
94
|
+
"rainfall_rate_32bit",
|
|
95
|
+
"rainfall_accumulated_32bit",
|
|
96
|
+
"weather_code_synop_4680",
|
|
97
|
+
"weather_code_synop_4677",
|
|
98
|
+
"reflectivity_32bit",
|
|
99
|
+
"mor_visibility",
|
|
100
|
+
"laser_amplitude",
|
|
101
|
+
"number_particles",
|
|
102
|
+
"sensor_temperature",
|
|
103
|
+
"sensor_heating_current",
|
|
104
|
+
"sensor_battery_voltage",
|
|
105
|
+
"sensor_status",
|
|
106
|
+
"rainfall_amount_absolute_32bit",
|
|
107
|
+
"error_code",
|
|
108
|
+
]
|
|
109
|
+
df_variables = df.iloc[:, 0:20]
|
|
110
|
+
df_variables.columns = column_names
|
|
111
|
+
|
|
112
|
+
# Define raw fields
|
|
113
|
+
df_raw_drop_concentration = df.iloc[:, 20:52].agg(",".join, axis=1)
|
|
114
|
+
df_raw_drop_average_velocity = df.iloc[:, 53:85].agg(",".join, axis=1)
|
|
115
|
+
df_raw_drop_number = df.iloc[:, 86:1110].agg(",".join, axis=1)
|
|
116
|
+
|
|
117
|
+
# Combine together
|
|
118
|
+
df = df_variables
|
|
119
|
+
df["raw_drop_concentration"] = df_raw_drop_concentration
|
|
120
|
+
df["raw_drop_average_velocity"] = df_raw_drop_average_velocity
|
|
121
|
+
df["raw_drop_number"] = df_raw_drop_number
|
|
122
|
+
|
|
123
|
+
# Drop invalid rows
|
|
124
|
+
df = df.loc[df["id"].astype(str).str.len() < 10]
|
|
125
|
+
|
|
126
|
+
# Drop columns not agreeing with DISDRODB L0 standards
|
|
127
|
+
columns_to_drop = [
|
|
128
|
+
"datalogger_temperature",
|
|
129
|
+
"datalogger_voltage",
|
|
130
|
+
"id",
|
|
131
|
+
"latitude",
|
|
132
|
+
"longitude",
|
|
133
|
+
]
|
|
134
|
+
df = df.drop(columns=columns_to_drop)
|
|
135
|
+
|
|
136
|
+
# Convert time column to datetime
|
|
137
|
+
df["time"] = pd.to_datetime(df["time"], format="%d-%m-%Y %H:%M:%S", errors="coerce")
|
|
138
|
+
|
|
139
|
+
# Return the dataframe adhering to DISDRODB L0 standards
|
|
140
|
+
return df
|
|
@@ -15,23 +15,21 @@
|
|
|
15
15
|
# You should have received a copy of the GNU General Public License
|
|
16
16
|
# along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
17
17
|
# -----------------------------------------------------------------------------.
|
|
18
|
-
|
|
18
|
+
"""Reader for RACLETS WJF 2019 station."""
|
|
19
|
+
import pandas as pd
|
|
20
|
+
|
|
19
21
|
from disdrodb.l0.l0_reader import is_documented_by, reader_generic_docstring
|
|
22
|
+
from disdrodb.l0.l0a_processing import read_raw_text_file
|
|
20
23
|
|
|
21
24
|
|
|
22
25
|
@is_documented_by(reader_generic_docstring)
|
|
23
26
|
def reader(
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
station_name,
|
|
27
|
-
# Processing options
|
|
28
|
-
force=False,
|
|
29
|
-
verbose=False,
|
|
30
|
-
parallel=False,
|
|
31
|
-
debugging_mode=False,
|
|
27
|
+
filepath,
|
|
28
|
+
logger=None,
|
|
32
29
|
):
|
|
30
|
+
"""Reader."""
|
|
33
31
|
##------------------------------------------------------------------------.
|
|
34
|
-
####
|
|
32
|
+
#### Define column names
|
|
35
33
|
column_names = [
|
|
36
34
|
"id",
|
|
37
35
|
"latitude",
|
|
@@ -59,7 +57,7 @@ def reader(
|
|
|
59
57
|
]
|
|
60
58
|
|
|
61
59
|
##------------------------------------------------------------------------.
|
|
62
|
-
####
|
|
60
|
+
#### Define reader options
|
|
63
61
|
reader_kwargs = {}
|
|
64
62
|
# - Define delimiter
|
|
65
63
|
reader_kwargs["delimiter"] = ";"
|
|
@@ -77,59 +75,43 @@ def reader(
|
|
|
77
75
|
# - Available: gzip, bz2, zip
|
|
78
76
|
reader_kwargs["compression"] = "gzip"
|
|
79
77
|
# - Strings to recognize as NA/NaN and replace with standard NA flags
|
|
80
|
-
# - Already included:
|
|
81
|
-
#
|
|
82
|
-
#
|
|
78
|
+
# - Already included: '#N/A', '#N/A N/A', '#NA', '-1.#IND', '-1.#QNAN',
|
|
79
|
+
# '-NaN', '-nan', '1.#IND', '1.#QNAN', '<NA>', 'N/A',
|
|
80
|
+
# 'NA', 'NULL', 'NaN', 'n/a', 'nan', 'null'
|
|
83
81
|
reader_kwargs["na_values"] = ["na", "", "error"]
|
|
84
82
|
|
|
85
83
|
##------------------------------------------------------------------------.
|
|
86
|
-
####
|
|
87
|
-
|
|
88
|
-
|
|
89
|
-
|
|
90
|
-
|
|
91
|
-
|
|
92
|
-
|
|
93
|
-
|
|
94
|
-
# - Split th 'TO_BE_SPLITTED' column
|
|
95
|
-
df_splitted = df["TO_BE_SPLITTED"].str.split(",", expand=True, n=1)
|
|
96
|
-
df_splitted.columns = ["datalogger_voltage", "rainfall_rate_32bit"]
|
|
97
|
-
df["rainfall_rate_32bit"] = df_splitted["rainfall_rate_32bit"]
|
|
84
|
+
#### Read the data
|
|
85
|
+
df = read_raw_text_file(
|
|
86
|
+
filepath=filepath,
|
|
87
|
+
column_names=column_names,
|
|
88
|
+
reader_kwargs=reader_kwargs,
|
|
89
|
+
logger=logger,
|
|
90
|
+
)
|
|
98
91
|
|
|
99
|
-
|
|
100
|
-
|
|
101
|
-
|
|
102
|
-
|
|
103
|
-
"longitude",
|
|
104
|
-
"datalogger_temperature",
|
|
105
|
-
"datalogger_error",
|
|
106
|
-
"TO_BE_SPLITTED",
|
|
107
|
-
]
|
|
108
|
-
df = df.drop(columns=columns_to_drop)
|
|
92
|
+
##------------------------------------------------------------------------.
|
|
93
|
+
#### Adapt the dataframe to adhere to DISDRODB L0 standards
|
|
94
|
+
# Drop invalid rows
|
|
95
|
+
df = df.loc[df["id"].astype(str).str.len() < 10]
|
|
109
96
|
|
|
110
|
-
|
|
111
|
-
|
|
97
|
+
# Split th 'TO_BE_SPLITTED' column
|
|
98
|
+
df_splitted = df["TO_BE_SPLITTED"].str.split(",", expand=True, n=1)
|
|
99
|
+
df_splitted.columns = ["datalogger_voltage", "rainfall_rate_32bit"]
|
|
100
|
+
df["rainfall_rate_32bit"] = df_splitted["rainfall_rate_32bit"]
|
|
112
101
|
|
|
113
|
-
|
|
102
|
+
# Drop columns not agreeing with DISDRODB L0 standards
|
|
103
|
+
columns_to_drop = [
|
|
104
|
+
"id",
|
|
105
|
+
"latitude",
|
|
106
|
+
"longitude",
|
|
107
|
+
"datalogger_temperature",
|
|
108
|
+
"datalogger_error",
|
|
109
|
+
"TO_BE_SPLITTED",
|
|
110
|
+
]
|
|
111
|
+
df = df.drop(columns=columns_to_drop)
|
|
114
112
|
|
|
115
|
-
|
|
116
|
-
|
|
117
|
-
glob_patterns = "*.dat.gz*"
|
|
113
|
+
# Convert time column to datetime
|
|
114
|
+
df["time"] = pd.to_datetime(df["time"], format="%d-%m-%Y %H:%M:%S", errors="coerce")
|
|
118
115
|
|
|
119
|
-
|
|
120
|
-
|
|
121
|
-
run_l0a(
|
|
122
|
-
raw_dir=raw_dir,
|
|
123
|
-
processed_dir=processed_dir,
|
|
124
|
-
station_name=station_name,
|
|
125
|
-
# Custom arguments of the reader for L0A processing
|
|
126
|
-
glob_patterns=glob_patterns,
|
|
127
|
-
column_names=column_names,
|
|
128
|
-
reader_kwargs=reader_kwargs,
|
|
129
|
-
df_sanitizer_fun=df_sanitizer_fun,
|
|
130
|
-
# Processing options
|
|
131
|
-
force=force,
|
|
132
|
-
verbose=verbose,
|
|
133
|
-
parallel=parallel,
|
|
134
|
-
debugging_mode=debugging_mode,
|
|
135
|
-
)
|
|
116
|
+
# Return the dataframe adhering to DISDRODB L0 standards
|
|
117
|
+
return df
|
|
@@ -15,23 +15,21 @@
|
|
|
15
15
|
# You should have received a copy of the GNU General Public License
|
|
16
16
|
# along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
17
17
|
# -----------------------------------------------------------------------------.
|
|
18
|
-
|
|
18
|
+
"""Reader for RIETHOLZBACH 2011 campaign."""
|
|
19
|
+
import pandas as pd
|
|
20
|
+
|
|
19
21
|
from disdrodb.l0.l0_reader import is_documented_by, reader_generic_docstring
|
|
22
|
+
from disdrodb.l0.l0a_processing import read_raw_text_file
|
|
20
23
|
|
|
21
24
|
|
|
22
25
|
@is_documented_by(reader_generic_docstring)
|
|
23
26
|
def reader(
|
|
24
|
-
|
|
25
|
-
|
|
26
|
-
station_name,
|
|
27
|
-
# Processing options
|
|
28
|
-
force=False,
|
|
29
|
-
verbose=False,
|
|
30
|
-
parallel=False,
|
|
31
|
-
debugging_mode=False,
|
|
27
|
+
filepath,
|
|
28
|
+
logger=None,
|
|
32
29
|
):
|
|
30
|
+
"""Reader."""
|
|
33
31
|
##------------------------------------------------------------------------.
|
|
34
|
-
####
|
|
32
|
+
#### Define column names
|
|
35
33
|
column_names = [
|
|
36
34
|
"time",
|
|
37
35
|
"id",
|
|
@@ -58,7 +56,7 @@ def reader(
|
|
|
58
56
|
]
|
|
59
57
|
|
|
60
58
|
##------------------------------------------------------------------------.
|
|
61
|
-
####
|
|
59
|
+
#### Define reader options
|
|
62
60
|
reader_kwargs = {}
|
|
63
61
|
# - Define delimiter
|
|
64
62
|
reader_kwargs["delimiter"] = ","
|
|
@@ -74,49 +72,34 @@ def reader(
|
|
|
74
72
|
# - Available: gzip, bz2, zip
|
|
75
73
|
reader_kwargs["compression"] = "infer"
|
|
76
74
|
# - Strings to recognize as NA/NaN and replace with standard NA flags
|
|
77
|
-
# - Already included:
|
|
78
|
-
#
|
|
79
|
-
#
|
|
75
|
+
# - Already included: '#N/A', '#N/A N/A', '#NA', '-1.#IND', '-1.#QNAN',
|
|
76
|
+
# '-NaN', '-nan', '1.#IND', '1.#QNAN', '<NA>', 'N/A',
|
|
77
|
+
# 'NA', 'NULL', 'NaN', 'n/a', 'nan', 'null'
|
|
80
78
|
reader_kwargs["na_values"] = ["na", "", "error", "NA"]
|
|
81
79
|
|
|
82
80
|
##------------------------------------------------------------------------.
|
|
83
|
-
####
|
|
84
|
-
|
|
85
|
-
|
|
86
|
-
import pandas as pd
|
|
87
|
-
|
|
88
|
-
# - Convert time column to datetime
|
|
89
|
-
df["time"] = pd.to_datetime(df["time"], format="%Y-%m-%d %H:%M:%S", errors="coerce")
|
|
90
|
-
|
|
91
|
-
# - Drop columns not agreeing with DISDRODB L0 standards
|
|
92
|
-
columns_to_drop = [
|
|
93
|
-
"id",
|
|
94
|
-
"datalogger_temperature",
|
|
95
|
-
"datalogger_voltage",
|
|
96
|
-
"datalogger_debug",
|
|
97
|
-
"datalogger_error",
|
|
98
|
-
]
|
|
99
|
-
df = df.drop(columns=columns_to_drop)
|
|
100
|
-
return df
|
|
101
|
-
|
|
102
|
-
##------------------------------------------------------------------------.
|
|
103
|
-
#### - Define glob pattern to search data files in <raw_dir>/data/<station_name>
|
|
104
|
-
glob_patterns = "*.dat*"
|
|
105
|
-
|
|
106
|
-
####----------------------------------------------------------------------.
|
|
107
|
-
#### - Create L0A products
|
|
108
|
-
run_l0a(
|
|
109
|
-
raw_dir=raw_dir,
|
|
110
|
-
processed_dir=processed_dir,
|
|
111
|
-
station_name=station_name,
|
|
112
|
-
# Custom arguments of the reader for L0A processing
|
|
113
|
-
glob_patterns=glob_patterns,
|
|
81
|
+
#### Read the data
|
|
82
|
+
df = read_raw_text_file(
|
|
83
|
+
filepath=filepath,
|
|
114
84
|
column_names=column_names,
|
|
115
85
|
reader_kwargs=reader_kwargs,
|
|
116
|
-
|
|
117
|
-
# Processing options
|
|
118
|
-
force=force,
|
|
119
|
-
verbose=verbose,
|
|
120
|
-
parallel=parallel,
|
|
121
|
-
debugging_mode=debugging_mode,
|
|
86
|
+
logger=logger,
|
|
122
87
|
)
|
|
88
|
+
|
|
89
|
+
##------------------------------------------------------------------------.
|
|
90
|
+
#### Adapt the dataframe to adhere to DISDRODB L0 standards
|
|
91
|
+
# Convert time column to datetime
|
|
92
|
+
df["time"] = pd.to_datetime(df["time"], format="%Y-%m-%d %H:%M:%S", errors="coerce")
|
|
93
|
+
|
|
94
|
+
# Drop columns not agreeing with DISDRODB L0 standards
|
|
95
|
+
columns_to_drop = [
|
|
96
|
+
"id",
|
|
97
|
+
"datalogger_temperature",
|
|
98
|
+
"datalogger_voltage",
|
|
99
|
+
"datalogger_debug",
|
|
100
|
+
"datalogger_error",
|
|
101
|
+
]
|
|
102
|
+
df = df.drop(columns=columns_to_drop)
|
|
103
|
+
|
|
104
|
+
# Return the dataframe adhering to DISDRODB L0 standards
|
|
105
|
+
return df
|
|
@@ -0,0 +1,117 @@
|
|
|
1
|
+
#!/usr/bin/env python3
|
|
2
|
+
# -----------------------------------------------------------------------------.
|
|
3
|
+
# Copyright (c) 2021-2023 DISDRODB developers
|
|
4
|
+
#
|
|
5
|
+
# This program is free software: you can redistribute it and/or modify
|
|
6
|
+
# it under the terms of the GNU General Public License as published by
|
|
7
|
+
# the Free Software Foundation, either version 3 of the License, or
|
|
8
|
+
# (at your option) any later version.
|
|
9
|
+
#
|
|
10
|
+
# This program is distributed in the hope that it will be useful,
|
|
11
|
+
# but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
12
|
+
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
13
|
+
# GNU General Public License for more details.
|
|
14
|
+
#
|
|
15
|
+
# You should have received a copy of the GNU General Public License
|
|
16
|
+
# along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
17
|
+
# -----------------------------------------------------------------------------.
|
|
18
|
+
"""Reader for SAMOYLOV 2017 campaign."""
|
|
19
|
+
import numpy as np
|
|
20
|
+
import pandas as pd
|
|
21
|
+
|
|
22
|
+
from disdrodb.l0.l0_reader import is_documented_by, reader_generic_docstring
|
|
23
|
+
from disdrodb.l0.l0a_processing import read_raw_text_file
|
|
24
|
+
|
|
25
|
+
|
|
26
|
+
@is_documented_by(reader_generic_docstring)
|
|
27
|
+
def reader(
|
|
28
|
+
filepath,
|
|
29
|
+
logger=None,
|
|
30
|
+
):
|
|
31
|
+
"""Reader."""
|
|
32
|
+
##------------------------------------------------------------------------.
|
|
33
|
+
#### Define column names
|
|
34
|
+
column_names = [
|
|
35
|
+
"time",
|
|
36
|
+
"rainfall_rate_32bit",
|
|
37
|
+
"rainfall_accumulated_32bit",
|
|
38
|
+
"weather_code_synop_4680",
|
|
39
|
+
"weather_code_synop_4677",
|
|
40
|
+
"reflectivity_32bit",
|
|
41
|
+
"mor_visibility",
|
|
42
|
+
"laser_amplitude",
|
|
43
|
+
"number_particles",
|
|
44
|
+
"sensor_temperature",
|
|
45
|
+
"sensor_heating_current",
|
|
46
|
+
"sensor_battery_voltage",
|
|
47
|
+
"sensor_status",
|
|
48
|
+
"rainfall_amount_absolute_32bit",
|
|
49
|
+
"error_code",
|
|
50
|
+
"raw_drop_concentration",
|
|
51
|
+
"raw_drop_average_velocity",
|
|
52
|
+
"raw_drop_number",
|
|
53
|
+
]
|
|
54
|
+
|
|
55
|
+
##------------------------------------------------------------------------.
|
|
56
|
+
#### Define reader options
|
|
57
|
+
reader_kwargs = {}
|
|
58
|
+
|
|
59
|
+
# - Define delimiter
|
|
60
|
+
reader_kwargs["delimiter"] = ";"
|
|
61
|
+
|
|
62
|
+
# Skip first row as columns names
|
|
63
|
+
reader_kwargs["header"] = None
|
|
64
|
+
|
|
65
|
+
# - Avoid first column to become df index !!!
|
|
66
|
+
reader_kwargs["index_col"] = False
|
|
67
|
+
|
|
68
|
+
# - Define behaviour when encountering bad lines
|
|
69
|
+
reader_kwargs["on_bad_lines"] = "skip"
|
|
70
|
+
|
|
71
|
+
# - Avoid UnicodeDecodeError: 'utf-8' codec can't decode byte 0xf0
|
|
72
|
+
reader_kwargs["encoding"] = "latin-1"
|
|
73
|
+
|
|
74
|
+
# - Define reader engine
|
|
75
|
+
# - C engine is faster
|
|
76
|
+
# - Python engine is more feature-complete
|
|
77
|
+
reader_kwargs["engine"] = "python"
|
|
78
|
+
|
|
79
|
+
# - Define on-the-fly decompression of on-disk data
|
|
80
|
+
# - Available: gzip, bz2, zip
|
|
81
|
+
reader_kwargs["compression"] = "infer"
|
|
82
|
+
|
|
83
|
+
# - Strings to recognize as NA/NaN and replace with standard NA flags
|
|
84
|
+
# - Already included: '#N/A', '#N/A N/A', '#NA', '-1.#IND', '-1.#QNAN',
|
|
85
|
+
# '-NaN', '-nan', '1.#IND', '1.#QNAN', '<NA>', 'N/A',
|
|
86
|
+
# 'NA', 'NULL', 'NaN', 'n/a', 'nan', 'null'
|
|
87
|
+
reader_kwargs["na_values"] = ["na", "", "error"]
|
|
88
|
+
|
|
89
|
+
# Different encoding for this campaign
|
|
90
|
+
reader_kwargs["encoding"] = "latin-1"
|
|
91
|
+
|
|
92
|
+
##------------------------------------------------------------------------.
|
|
93
|
+
#### Read the data
|
|
94
|
+
df = read_raw_text_file(
|
|
95
|
+
filepath=filepath,
|
|
96
|
+
column_names=column_names,
|
|
97
|
+
reader_kwargs=reader_kwargs,
|
|
98
|
+
logger=logger,
|
|
99
|
+
)
|
|
100
|
+
|
|
101
|
+
##------------------------------------------------------------------------.
|
|
102
|
+
#### Adapt the dataframe to adhere to DISDRODB L0 standards
|
|
103
|
+
|
|
104
|
+
# Special parsing if 'Error in data reading' in rainfall_rate_32bit column
|
|
105
|
+
if np.any(df["rainfall_rate_32bit"].str.startswith("Error in data reading!", na=False)):
|
|
106
|
+
df["rainfall_rate_32bit"] = df["rainfall_rate_32bit"].str.replace("Error in data reading!", "")
|
|
107
|
+
df["raw_drop_number"] = df["raw_drop_average_velocity"]
|
|
108
|
+
df["raw_drop_average_velocity"] = df["raw_drop_concentration"]
|
|
109
|
+
df["raw_drop_concentration"] = df["error_code"]
|
|
110
|
+
df["error_code"] = df["rainfall_amount_absolute_32bit"].str[7:]
|
|
111
|
+
df["rainfall_amount_absolute_32bit"] = df["rainfall_amount_absolute_32bit"].str[:7]
|
|
112
|
+
|
|
113
|
+
# Convert time column to datetime
|
|
114
|
+
df["time"] = pd.to_datetime(df["time"], format="%d/%m/%Y %H:%M:%S", errors="coerce")
|
|
115
|
+
|
|
116
|
+
# Return the dataframe adhering to DISDRODB L0 standards
|
|
117
|
+
return df
|
|
@@ -0,0 +1,137 @@
|
|
|
1
|
+
#!/usr/bin/env python3
|
|
2
|
+
# -----------------------------------------------------------------------------.
|
|
3
|
+
# Copyright (c) 2021-2023 DISDRODB developers
|
|
4
|
+
#
|
|
5
|
+
# This program is free software: you can redistribute it and/or modify
|
|
6
|
+
# it under the terms of the GNU General Public License as published by
|
|
7
|
+
# the Free Software Foundation, either version 3 of the License, or
|
|
8
|
+
# (at your option) any later version.
|
|
9
|
+
#
|
|
10
|
+
# This program is distributed in the hope that it will be useful,
|
|
11
|
+
# but WITHOUT ANY WARRANTY; without even the implied warranty of
|
|
12
|
+
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
|
|
13
|
+
# GNU General Public License for more details.
|
|
14
|
+
#
|
|
15
|
+
# You should have received a copy of the GNU General Public License
|
|
16
|
+
# along with this program. If not, see <http://www.gnu.org/licenses/>.
|
|
17
|
+
# -----------------------------------------------------------------------------.
|
|
18
|
+
"""Reader for SAMOYLOV 2019 campaign."""
|
|
19
|
+
import pandas as pd
|
|
20
|
+
|
|
21
|
+
from disdrodb.l0.l0_reader import is_documented_by, reader_generic_docstring
|
|
22
|
+
from disdrodb.l0.l0a_processing import read_raw_text_file
|
|
23
|
+
|
|
24
|
+
|
|
25
|
+
@is_documented_by(reader_generic_docstring)
|
|
26
|
+
def reader(
|
|
27
|
+
filepath,
|
|
28
|
+
logger=None,
|
|
29
|
+
):
|
|
30
|
+
"""Reader."""
|
|
31
|
+
##------------------------------------------------------------------------.
|
|
32
|
+
#### Define column names
|
|
33
|
+
column_names = [
|
|
34
|
+
"id",
|
|
35
|
+
"latitude",
|
|
36
|
+
"longitude",
|
|
37
|
+
"time",
|
|
38
|
+
"datalogger_temperature",
|
|
39
|
+
"TO_BE_SPLITTED", # datalogger_voltage and rainfall_rate_32bit
|
|
40
|
+
"rainfall_accumulated_32bit",
|
|
41
|
+
"weather_code_synop_4680",
|
|
42
|
+
"weather_code_synop_4677",
|
|
43
|
+
"reflectivity_32bit",
|
|
44
|
+
"mor_visibility",
|
|
45
|
+
"laser_amplitude",
|
|
46
|
+
"number_particles",
|
|
47
|
+
"sensor_temperature",
|
|
48
|
+
"sensor_heating_current",
|
|
49
|
+
"sensor_battery_voltage",
|
|
50
|
+
"sensor_status",
|
|
51
|
+
"rainfall_amount_absolute_32bit",
|
|
52
|
+
"error_code",
|
|
53
|
+
"raw_drop_concentration",
|
|
54
|
+
"raw_drop_average_velocity",
|
|
55
|
+
"raw_drop_number",
|
|
56
|
+
"datalogger_error",
|
|
57
|
+
]
|
|
58
|
+
|
|
59
|
+
##------------------------------------------------------------------------.
|
|
60
|
+
#### Define reader options
|
|
61
|
+
reader_kwargs = {}
|
|
62
|
+
|
|
63
|
+
# - Define delimiter
|
|
64
|
+
reader_kwargs["delimiter"] = ";"
|
|
65
|
+
|
|
66
|
+
# Skip first row as columns names
|
|
67
|
+
reader_kwargs["header"] = None
|
|
68
|
+
|
|
69
|
+
# - Avoid first column to become df index !!!
|
|
70
|
+
reader_kwargs["index_col"] = False
|
|
71
|
+
|
|
72
|
+
# - Define behaviour when encountering bad lines
|
|
73
|
+
reader_kwargs["on_bad_lines"] = "skip"
|
|
74
|
+
|
|
75
|
+
# - Define encoding
|
|
76
|
+
# --> Avoid UnicodeDecodeError: 'utf-8' codec can't decode byte 0xf0
|
|
77
|
+
reader_kwargs["encoding"] = "latin-1"
|
|
78
|
+
|
|
79
|
+
# - Define reader engine
|
|
80
|
+
# - C engine is faster
|
|
81
|
+
# - Python engine is more feature-complete
|
|
82
|
+
reader_kwargs["engine"] = "python"
|
|
83
|
+
|
|
84
|
+
# - Define on-the-fly decompression of on-disk data
|
|
85
|
+
# - Available: gzip, bz2, zip
|
|
86
|
+
reader_kwargs["compression"] = "infer"
|
|
87
|
+
|
|
88
|
+
# - Strings to recognize as NA/NaN and replace with standard NA flags
|
|
89
|
+
# - Already included: '#N/A', '#N/A N/A', '#NA', '-1.#IND', '-1.#QNAN',
|
|
90
|
+
# '-NaN', '-nan', '1.#IND', '1.#QNAN', '<NA>', 'N/A',
|
|
91
|
+
# 'NA', 'NULL', 'NaN', 'n/a', 'nan', 'null'
|
|
92
|
+
reader_kwargs["na_values"] = [
|
|
93
|
+
"na",
|
|
94
|
+
"",
|
|
95
|
+
"error",
|
|
96
|
+
"NA",
|
|
97
|
+
]
|
|
98
|
+
|
|
99
|
+
##------------------------------------------------------------------------.
|
|
100
|
+
#### Read the data
|
|
101
|
+
df = read_raw_text_file(
|
|
102
|
+
filepath=filepath,
|
|
103
|
+
column_names=column_names,
|
|
104
|
+
reader_kwargs=reader_kwargs,
|
|
105
|
+
logger=logger,
|
|
106
|
+
)
|
|
107
|
+
|
|
108
|
+
##------------------------------------------------------------------------.
|
|
109
|
+
#### Adapt the dataframe to adhere to DISDRODB L0 standards
|
|
110
|
+
# Convert time column to datetime
|
|
111
|
+
df["time"] = pd.to_datetime(df["time"], format="%d-%m-%Y %H:%M:%S", errors="coerce")
|
|
112
|
+
|
|
113
|
+
# Drop rows with None value
|
|
114
|
+
df = df[~df["TO_BE_SPLITTED"].isna()]
|
|
115
|
+
|
|
116
|
+
# Drop rows when 'Error in data reading' in TO_BE_SPLITTED column
|
|
117
|
+
bad_indices = df[df["TO_BE_SPLITTED"].str.contains("Error in data reading!")].index
|
|
118
|
+
df = df.drop(bad_indices)
|
|
119
|
+
|
|
120
|
+
# Split TO_BE_SPLITTED columns
|
|
121
|
+
df_splitted = df["TO_BE_SPLITTED"].str.split(",", expand=True, n=1)
|
|
122
|
+
df_splitted.columns = ["datalogger_voltage", "rainfall_rate_32bit"]
|
|
123
|
+
df["rainfall_rate_32bit"] = df_splitted["rainfall_rate_32bit"]
|
|
124
|
+
|
|
125
|
+
# Drop columns not agreeing with DISDRODB L0 standards
|
|
126
|
+
columns_to_drop = [
|
|
127
|
+
"id",
|
|
128
|
+
"TO_BE_SPLITTED",
|
|
129
|
+
"datalogger_temperature",
|
|
130
|
+
"datalogger_error",
|
|
131
|
+
"latitude",
|
|
132
|
+
"longitude",
|
|
133
|
+
]
|
|
134
|
+
df = df.drop(columns=columns_to_drop)
|
|
135
|
+
|
|
136
|
+
# Return the dataframe adhering to DISDRODB L0 standards
|
|
137
|
+
return df
|