graphica-plot 2.0.0__py3-none-any.whl

graphica/gui/workers.py

@@ -0,0 +1,280 @@
+"""データファイルの読み込み(TaskRunner で別スレッドで動かす関数)と、文字コード・区切り文字の判定。"""
+import csv
+import dataclasses
+import re
+import warnings
+
+import pandas as pd
+
+# 順に試す文字コード。utf-8-sig は BOM の有無どちらの UTF-8 も読める。latin-1 は必ず読める最後の手段。
+# UTF-16 は入れない(cp932 や latin-1 が例外を出さずに化けた結果を返す)。BOM で別に見分ける。
+CSV_ENCODING_FALLBACKS = ['utf-8-sig', 'cp932', 'latin-1']
+
+# 組み込みで読める拡張子。開くダイアログ、ドラッグ&ドロップ、フォルダ一括取り込み、プレビューがすべてここを見る。
+# 区切り文字付きテキストは CSV と同じく文字コードと区切りを判定する。.xlsx は openpyxl、.xls は xlrd。
+DELIMITED_TEXT_EXTENSIONS = ('.csv', '.txt')
+EXCEL_EXTENSIONS = ('.xlsx', '.xls')
+BUILTIN_DATA_FILE_EXTENSIONS = DELIMITED_TEXT_EXTENSIONS + EXCEL_EXTENSIONS
+
+
+def is_delimited_text_file(file_path):
+    return bool(file_path) and file_path.lower().endswith(DELIMITED_TEXT_EXTENSIONS)
+
+
+def is_excel_file(file_path):
+    return bool(file_path) and file_path.lower().endswith(EXCEL_EXTENSIONS)
+
+
+def excel_engine_for(file_path):
+    """openpyxl は .xls を読めないので、拡張子で xlrd と使い分ける。"""
+    return 'xlrd' if file_path.lower().endswith('.xls') else 'openpyxl'
+
+
+def pandas_separator(delimiter):
+    """(sep, engine)。空白は空白の連続とみなす(桁揃えの空白で空の列が大量にできないように)。正規表現は python エンジン。"""
+    if delimiter == ' ':
+        return r'\s+', 'python'
+    if len(delimiter) == 1:
+        return delimiter, 'c'
+    return delimiter, 'python'
+
+
+def _detect_bom_encoding(file_path):
+    """BOM から文字コードを判定する。無ければ None。"""
+    with open(file_path, 'rb') as f:
+        head = f.read(4)
+    if head.startswith(b'\xff\xfe') or head.startswith(b'\xfe\xff'):
+        return 'utf-16'
+    if head.startswith(b'\xef\xbb\xbf'):
+        return 'utf-8-sig'
+    return None
+
+
+def _csv_encoding_candidates(file_path):
+    """BOM の判定を先頭に、CSV_ENCODING_FALLBACKS を並べた候補(読み込みと判定で順序を共有する)。"""
+    bom_encoding = _detect_bom_encoding(file_path)
+    candidates = [bom_encoding] if bom_encoding else []
+    candidates += [e for e in CSV_ENCODING_FALLBACKS if e not in candidates]
+    return candidates
+
+
+def detect_csv_encoding(file_path):
+    """テキストとしてデコードできる最初の文字コード。DataFrame にはしないので軽い。"""
+    last_error = None
+    for encoding in _csv_encoding_candidates(file_path):
+        try:
+            with open(file_path, encoding=encoding) as f:
+                f.read()
+            return encoding
+        except UnicodeDecodeError as e:
+            last_error = e
+            continue
+    raise ValueError(
+        f"CSVファイルの文字コードを判定できませんでした "
+        f"(試行: {', '.join(CSV_ENCODING_FALLBACKS)})。詳細: {last_error}"
+    )
+
+
+def _sniff_delimiter_from_text(sample_text):
+    """カンマ・タブ・セミコロン・空白のどれかを推測する。判定できなければカンマ。"""
+    if not sample_text.strip():
+        return ','
+    try:
+        dialect = csv.Sniffer().sniff(sample_text, delimiters=',\t; ')
+        return dialect.delimiter
+    except csv.Error:
+        return ','
+
+
+def detect_csv_delimiter(file_path, encoding, sample_lines=50):
+    lines = []
+    try:
+        with open(file_path, encoding=encoding, errors='replace') as f:
+            for _ in range(sample_lines):
+                line = f.readline()
+                if not line:
+                    break
+                lines.append(line)
+    except OSError:
+        return ','
+    return _sniff_delimiter_from_text(''.join(lines))
+
+
+def detect_clipboard_delimiter(text, sample_lines=50):
+    """クリップボードの区切り文字(Excel からのコピーはタブ、CSV の文字列ならカンマなど)。"""
+    sample = ''.join(text.splitlines(keepends=True)[:sample_lines])
+    return _sniff_delimiter_from_text(sample)
+
+
+# 測定装置の書き出し(JASCO の TXT など)は、数値の表の前後に測定条件の行が付く。そういうファイルでだけ数値の表を探す
+NUMERIC_TABLE_DELIMITERS = ('\t', ',', ';', ' ')
+MIN_NUMERIC_TABLE_ROWS = 3
+
+
+@dataclasses.dataclass
+class NumericTable:
+    first_line: int  # 0 始まり、表の最初の数値の行
+    stop_line: int  # 表の最後の行の次
+    delimiter: str
+    header: list | None  # 表の直前の行が同じ列数の文字なら列名
+
+
+def _split_numeric_row(line, delimiter):
+    """2 列以上がすべて数値として読めれば、その列の文字列。読めなければ None。"""
+    text = line.strip()
+    fields = text.split() if delimiter == ' ' else [f.strip() for f in text.split(delimiter)]
+    if len(fields) < 2:
+        return None
+    try:
+        for field in fields:
+            float(field)
+    except ValueError:
+        return None
+    return fields
+
+
+def find_numeric_table(lines):
+    """列数のそろった数値の行が最も長く続くところ。MIN_NUMERIC_TABLE_ROWS 行に満たなければ None。"""
+    best = None
+    for delimiter in NUMERIC_TABLE_DELIMITERS:
+        start, columns = None, 0
+        for index, line in enumerate([*lines, '']):
+            fields = _split_numeric_row(line, delimiter)
+            if fields is not None and start is not None and len(fields) == columns:
+                continue
+            if start is not None and (best is None or index - start > best.stop_line - best.first_line):
+                best = NumericTable(start, index, delimiter, None)
+            start, columns = (index, len(fields)) if fields is not None else (None, 0)
+    if best is None or best.stop_line - best.first_line < MIN_NUMERIC_TABLE_ROWS:
+        return None
+    if best.first_line > 0:
+        above = lines[best.first_line - 1].strip()
+        names = above.split() if best.delimiter == ' ' else [f.strip() for f in above.split(best.delimiter)]
+        column_count = len(_split_numeric_row(lines[best.first_line], best.delimiter))
+        if above and len(names) == column_count and _split_numeric_row(above, best.delimiter) is None:
+            best.header = names
+    return best
+
+
+def read_numeric_table(file_path, encoding):
+    """前後の説明の行を除いた数値の表を DataFrame にする。(DataFrame, NumericTable)、見つからなければ None。"""
+    import io
+
+    with open(file_path, encoding=encoding) as f:
+        lines = f.read().splitlines()
+    table = find_numeric_table(lines)
+    if table is None:
+        return None
+    sep, engine = pandas_separator(table.delimiter)
+    body = '\n'.join(line.strip() for line in lines[table.first_line:table.stop_line])
+    df = pd.read_csv(io.StringIO(body), sep=sep, engine=engine, header=None)
+    names = table.header or [f"列{i + 1}" for i in range(df.shape[1])]
+    # 同じ列名は read_csv と同じく .1, .2 を付けて分ける
+    seen = {}
+    unique = []
+    for name in names:
+        count = seen.get(name, 0)
+        unique.append(name if count == 0 else f"{name}.{count}")
+        seen[name] = count + 1
+    df.columns = unique
+    return df, table
+
+
+def has_numeric_column(df):
+    return any(pd.api.types.is_numeric_dtype(df[column]) for column in df.columns)
+
+
+# 年(4 桁の数字)を含む値だけを日付とみなす。"12:30" のような時刻だけの文字は、読んだ日の日付が補われてしまう
+_YEAR_PATTERN = re.compile(r"\d{4}")
+
+
+def parse_date_columns(df):
+    """文字の列のうち、空でない値がすべて日付として読める列を日付型にする(Excel の日付の列と同じ扱い)。"""
+    for column in df.columns:
+        series = df[column]
+        if series.dtype != object:
+            continue
+        text = series.dropna().astype(str).str.strip()
+        text = text[text != ""]
+        if text.empty or not text.map(lambda value: _YEAR_PATTERN.search(value) is not None).all():
+            continue
+        with warnings.catch_warnings():
+            # 書式を推測できない警告と、タイムゾーンが混ざった列の将来の仕様変更の警告(その列は下で日付にしない)
+            warnings.simplefilter("ignore", UserWarning)
+            warnings.simplefilter("ignore", FutureWarning)
+            try:
+                # 書式は値ごとに読む(先頭の値の書式に合わせると、時刻のある値と無い値が混ざった列を読めない)
+                parsed = pd.to_datetime(text, errors="coerce", format="mixed")
+            except (ValueError, TypeError):
+                continue
+        # タイムゾーンが混ざると日付型にならない(object のまま)
+        if parsed.isna().any() or not pd.api.types.is_datetime64_any_dtype(parsed):
+            continue
+        converted = pd.Series(pd.NaT, index=series.index, dtype=parsed.dtype)
+        converted.loc[parsed.index] = parsed
+        df[column] = converted
+    return df
+
+
+def read_data_file(file_path):
+    """データファイルを DataFrame にする。プラグインがその拡張子を登録していればそちらを使う。"""
+    ext = file_path.lower().split('.')[-1]
+
+    from graphica.core.plugin_api import get_plugin_api
+    from graphica.core.plugin_types import PluginExecutionError
+    api = get_plugin_api()
+    importer = api.get_importer_for_extension(ext) if api is not None else None
+    if importer is not None:
+        try:
+            result = importer.loader(file_path)
+        except Exception as e:
+            raise PluginExecutionError(importer.name, f"「{file_path}」の読み込みに失敗しました: {e}") from e
+        if not isinstance(result, pd.DataFrame):
+            raise PluginExecutionError(
+                importer.name,
+                "現在サポートされているのは単一のDataFrameを返すインポーターのみです"
+                "(複数シートを返す形式は未対応です)。"
+            )
+        return result
+
+    if is_delimited_text_file(file_path):
+        # 区切り文字もここで推測する(カンマで読むとタブ区切りが1列に潰れ、「2列以上必要」で弾かれる)
+        last_error = None
+        for encoding in _csv_encoding_candidates(file_path):
+            try:
+                sep, engine = pandas_separator(detect_csv_delimiter(file_path, encoding))
+                df = pd.read_csv(file_path, header=0, encoding=encoding, sep=sep, engine=engine)
+            except UnicodeDecodeError as e:
+                last_error = e
+                continue
+            except pd.errors.ParserError as e:
+                try:
+                    found = read_numeric_table(file_path, encoding)
+                except UnicodeDecodeError:
+                    found = None
+                if found is not None:
+                    return parse_date_columns(found[0])
+                last_error = e
+                continue
+            # 普通に読めて数値の列があるファイルは今までどおり。どの列も数値にならないときだけ、前後の説明の行を疑う
+            if not has_numeric_column(df):
+                found = read_numeric_table(file_path, encoding)
+                if found is not None:
+                    return parse_date_columns(found[0])
+            return parse_date_columns(df)
+        raise ValueError(
+            f"テキストファイルの文字コードを判定できませんでした "
+            f"(試行: {', '.join(CSV_ENCODING_FALLBACKS)})。詳細: {last_error}"
+        )
+    elif is_excel_file(file_path):
+        return pd.read_excel(file_path, engine=excel_engine_for(file_path))
+    else:
+        raise ValueError(f"未対応のファイル形式です: {ext}")
+
+
+def load_data_file_task(file_path, report_progress=None, is_cancelled=None):
+    """TaskRunner 用。読み込みは中断できないので、report_progress と is_cancelled は受け取るだけ。"""
+    df = read_data_file(file_path)
+    if len(df.columns) < 2:
+        raise ValueError("データには少なくとも2列必要です。")
+    return df

graphica/models/project.py

@@ -0,0 +1,197 @@
+import hashlib
+import json
+import os
+import logging
+
+from PySide6.QtCore import QObject, Signal
+
+from graphica.core.dataset import Dataset
+from graphica.core.json_utils import GraphicaJSONEncoder
+
+logger = logging.getLogger(__name__)
+
+# 保存は .gra。.graphica(以前の拡張子)も中身は同じ JSON なので読み書きできる
+PROJECT_FILE_EXTENSION = '.gra'
+PROJECT_FILE_EXTENSIONS = ('.gra', '.graphica')
+LEGACY_PICKLE_EXTENSION = '.pkl'
+
+
+# format_version が無いファイルは 0。構造を変えたら上げて、前の版からの変換を _MIGRATIONS に足す
+CURRENT_FORMAT_VERSION = 1
+
+
+def _migrate_v0_to_v1(data):
+    """0 から 1 は format_version を足しただけで、構造は同じ。"""
+    return data
+
+
+# from_version -> (from_version + 1) の形にする関数
+_MIGRATIONS = {
+    0: _migrate_v0_to_v1,
+}
+
+
+def _migrate_project_data(data):
+    """CURRENT_FORMAT_VERSION まで順に変換する。新しい版のファイルはエラーにする(黙って半分だけ読まない)。"""
+    version = data.get('format_version', 0)
+    if version > CURRENT_FORMAT_VERSION:
+        raise ValueError(
+            f"このプロジェクトファイルはバージョン{version}で保存されていますが、"
+            f"このアプリケーションが対応しているのはバージョン{CURRENT_FORMAT_VERSION}までです。"
+            "アプリケーションを最新版に更新してください。"
+        )
+    while version < CURRENT_FORMAT_VERSION:
+        migrate = _MIGRATIONS.get(version)
+        if migrate is None:
+            raise ValueError(f"バージョン{version}からの移行手順が見つかりません。")
+        data = migrate(data)
+        version += 1
+    return data
+
+
+class ProjectModel(QObject):
+    """1つの文書(データセット、軸ごとの設定、レイアウト)。
+
+    既存の変更箇所は PlotterApp._update_plot() を直接呼ぶ。changed は、そのメソッドを知らない新しい経路が
+    再描画を頼むためのもので、今は誰も発行しない。
+    """
+
+    changed = Signal()
+
+    def __init__(self):
+        super().__init__()
+        # 最後に保存・読み込みしたパス(オートセーブでも変わる。上書き保存先は PlotterApp._current_project_path)
+        self.current_filepath = ""
+
+        self.datasets = []
+        # {'name', 'children': [...]} の入れ子。子はフォルダか {'dataset': Dataset}。name='' のルートは表示しない
+        self.dataset_group_tree = {'name': '', 'children': []}
+        self.all_plot_settings = []     # 軸ごとの設定(core/axis_settings.py)
+        self.active_axis_index = 0
+
+        self.layout_rows = 1
+        self.layout_cols = 1
+        # 'grid' か 'free'。'free' では all_plot_settings の数がサブプロットの数で、各設定の
+        # 'free_rect' に (left, bottom, width, height) を 0〜1 で持つ
+        self.layout_mode = 'grid'
+
+        # (a)(b)(c)… は描くときに並びから決めるので、文字は保存しない
+        self.panel_labels_enabled = False
+
+        # 同じ行(sharex)・同じ列(sharey)のサブプロットの軸を束ね、内側の目盛りラベルを隠す。グリッドのときだけ
+        self.share_x_axis = False
+        self.share_y_axis = False
+
+    def notify_changed(self):
+        self.changed.emit()
+
+    def save_project(self, filepath):
+        """JSON で保存する(開いてもコードが実行されない)。拡張子は .gra か .graphica。"""
+        ext = os.path.splitext(filepath)[1].lower()
+        if ext not in PROJECT_FILE_EXTENSIONS:
+            raise ValueError(f"サポートされていない拡張子です: {ext}")
+        self._save_project_json(filepath)
+        self.current_filepath = filepath
+
+    def load_project(self, filepath):
+        if not os.path.exists(filepath):
+            raise FileNotFoundError(f"ファイルが見つかりません: {filepath}")
+
+        ext = os.path.splitext(filepath)[1].lower()
+        if ext == LEGACY_PICKLE_EXTENSION:
+            raise ValueError("旧形式(.pkl)のプロジェクトには対応していません。")
+        if ext not in PROJECT_FILE_EXTENSIONS:
+            raise ValueError(f"サポートされていない拡張子です: {ext}")
+        self._load_project_json(filepath)
+        self.current_filepath = filepath
+
+    @staticmethod
+    def _tree_to_json(node):
+        """Dataset の参照を dataset_id の文字列にする(JSON にするため)。"""
+        if 'dataset' in node:
+            return {'dataset_id': node['dataset'].dataset_id}
+        return {
+            'name': node.get('name', ''),
+            'children': [ProjectModel._tree_to_json(child) for child in node.get('children', [])],
+        }
+
+    @staticmethod
+    def _tree_from_json(node, dataset_map):
+        """_tree_to_json() の逆。見つからない ID は警告して除く(壊れたファイルでも読めるように)。"""
+        if 'dataset_id' in node:
+            ds = dataset_map.get(node['dataset_id'])
+            if ds is None:
+                logger.warning(
+                    "dataset_group_tree内に存在しないdataset_idがあるため、"
+                    "このリーフをスキップします: %s", node['dataset_id']
+                )
+                return None
+            return {'dataset': ds}
+
+        children = []
+        for child in node.get('children', []):
+            converted = ProjectModel._tree_from_json(child, dataset_map)
+            if converted is not None:
+                children.append(converted)
+        return {'name': node.get('name', ''), 'children': children}
+
+    def content_fingerprint(self):
+        """保存すると書き出される内容のハッシュ(未保存の変更の判定)。
+
+        Undo の clean 状態では判定しない(軸の設定・列の計算・フォルダ操作は Undo の対象外)。
+        呼ぶ前に UI にしか無い状態を反映しておくこと(PlotterApp._sync_project_from_ui)。
+        """
+        payload = self._json_payload()
+        # 編集対象のサブプロットを切り替えただけで保存の確認が出ないように
+        payload.pop('active_axis_index', None)
+        text = json.dumps(payload, cls=GraphicaJSONEncoder, sort_keys=True, ensure_ascii=False)
+        return hashlib.sha256(text.encode('utf-8')).hexdigest()
+
+    def _json_payload(self):
+        """.graphica に書き出す辞書(保存と content_fingerprint で共有)。"""
+        return {
+            'format_version': CURRENT_FORMAT_VERSION,
+            'datasets': [ds.to_dict() for ds in self.datasets],
+            'dataset_group_tree': self._tree_to_json(self.dataset_group_tree),
+            'all_plot_settings': self.all_plot_settings,
+            'active_axis_index': self.active_axis_index,
+            'layout_rows': self.layout_rows,
+            'layout_cols': self.layout_cols,
+            'layout_mode': self.layout_mode,
+            'panel_labels_enabled': self.panel_labels_enabled,
+            'share_x_axis': self.share_x_axis,
+            'share_y_axis': self.share_y_axis,
+        }
+
+    def _save_project_json(self, filepath):
+        data = self._json_payload()
+        # 名前に日本語が多いので \uXXXX にせず読める形で残す
+        with open(filepath, 'w', encoding='utf-8') as f:
+            json.dump(data, f, cls=GraphicaJSONEncoder, indent=2, ensure_ascii=False)
+
+    def _load_project_json(self, filepath):
+        with open(filepath, 'r', encoding='utf-8') as f:
+            data = json.load(f)
+
+        data = _migrate_project_data(data)
+
+        self.datasets = [Dataset.from_dict(d) for d in data.get('datasets', [])]
+        dataset_map = {ds.dataset_id: ds for ds in self.datasets}
+
+        tree_data = data.get('dataset_group_tree')
+        if tree_data:
+            self.dataset_group_tree = self._tree_from_json(tree_data, dataset_map)
+        else:
+            # キーが無ければ、.pkl と同じく全部をルートに置く
+            self.dataset_group_tree = {
+                'name': '', 'children': [{'dataset': ds} for ds in self.datasets]
+            }
+
+        self.all_plot_settings = data.get('all_plot_settings', [])
+        self.active_axis_index = data.get('active_axis_index', 0)
+        self.layout_rows = data.get('layout_rows', 1)
+        self.layout_cols = data.get('layout_cols', 1)
+        self.layout_mode = data.get('layout_mode', 'grid')
+        self.panel_labels_enabled = data.get('panel_labels_enabled', False)
+        self.share_x_axis = data.get('share_x_axis', False)
+        self.share_y_axis = data.get('share_y_axis', False)

graphica/plugin/__init__.py

@@ -0,0 +1,20 @@
+"""
+プラグイン向けの公開の入口(プラグイン API 2.0)。
+
+プラグインが本体から import してよいのは、ここと graphica.plugin.testing だけ。
+本体の内部(core / gui / models)は予告なく動くが、ここの名前はプラグイン API の版で守る。
+"""
+from graphica.core.dataset import Dataset
+from graphica.core.plugin_api import GraphicaPluginAPI
+from graphica.core.plugin_context import PluginContext
+from graphica.core.plugin_manifest import PLUGIN_API_VERSION
+from graphica.core.plugin_types import AnalysisResult, PluginExecutionError
+
+__all__ = [
+    "PLUGIN_API_VERSION",
+    "AnalysisResult",
+    "Dataset",
+    "GraphicaPluginAPI",
+    "PluginContext",
+    "PluginExecutionError",
+]

graphica/plugin/testing.py

@@ -0,0 +1,57 @@
+"""
+プラグインのテスト用の公開の入口。Graphica を起動せずに試せる。
+
+    from graphica.plugin.testing import FakeGraphicaPluginAPI, FakePluginContext, load_plugin_like_graphica
+"""
+import os
+import shutil
+import tempfile
+
+from graphica.core.plugin_api import GraphicaPluginAPI, PluginManager
+from graphica.core.plugin_install import PluginInstallError, install_plugin_zip
+from graphica.core.plugin_testing import FakeGraphicaPluginAPI, FakePluginContext
+from typing import Any
+
+__all__ = [
+    "FakeGraphicaPluginAPI",
+    "FakePluginContext",
+    "PluginInstallError",
+    "install_zip_like_graphica",
+    "load_plugin_like_graphica",
+]
+
+
+def load_plugin_like_graphica(plugin_folder: str, work_dir: str | None = None) -> tuple[GraphicaPluginAPI, dict[str, Any]]:
+    """
+    本体と同じ経路(manifest の検証 → import → register(api))でプラグインを読み込む。
+    相対 import の誤りや api_version の不一致など、偽物の API では見つからない問題を確かめる。
+
+    Args:
+        plugin_folder (str): plugin.json と __init__.py のあるフォルダ。
+        work_dir (str | None): 写しを置く場所(pytest の tmp_path など)。省略時は一時フォルダ。
+    Returns:
+        tuple[GraphicaPluginAPI, dict]: 登録を受けた API と、読み込み結果
+            ({"name", "info", "error", "disabled"}。成功なら error は None)。
+            フック単位の失敗は api.registration_errors にある。
+    """
+    plugin_folder = os.path.normpath(plugin_folder)
+    name = os.path.basename(plugin_folder)
+    # 本体は探索フォルダの中のサブフォルダをすべてプラグインとみなすので、対象だけを写した場所で読む。
+    plugins_dir = tempfile.mkdtemp(prefix="graphica_plugins_", dir=work_dir)
+    shutil.copytree(plugin_folder, os.path.join(plugins_dir, name),
+                    ignore=shutil.ignore_patterns("__pycache__"))
+    api = GraphicaPluginAPI()
+    records = PluginManager(plugins_dir).load_all(api)
+    return api, records[0]
+
+
+def install_zip_like_graphica(zip_path: str, target_dir: str) -> str:
+    """
+    本体の「プラグインをインストール」と同じ処理で zip を target_dir に展開する。
+
+    Returns:
+        str: インストールされたプラグインのフォルダ名。
+    Raises:
+        PluginInstallError: zip の中身が本体の期待する形でない。
+    """
+    return install_plugin_zip(zip_path, target_dir=target_dir)

graphica/sample_data/cooling_curve_sample.csv

@@ -0,0 +1,42 @@
+時間_s,温度_degC
+0.0,90.18
+0.5,83.78
+1.0,79.74
+1.5,75.18
+2.0,69.18
+2.5,65.66
+3.0,62.96
+3.5,59.43
+4.0,56.63
+4.5,53.4
+5.0,51.95
+5.5,49.62
+6.0,47.11
+6.5,45.85
+7.0,43.72
+7.5,41.34
+8.0,40.62
+8.5,38.5
+9.0,38.39
+9.5,36.73
+10.0,35.63
+10.5,34.41
+11.0,34.71
+11.5,33.11
+12.0,32.24
+12.5,31.64
+13.0,31.58
+13.5,30.94
+14.0,30.48
+14.5,30.04
+15.0,30.65
+15.5,28.75
+16.0,28.34
+16.5,27.85
+17.0,28.42
+17.5,28.46
+18.0,27.48
+18.5,26.82
+19.0,26.63
+19.5,27.33
+20.0,27.22