WO2014176340A1 - Improved methods of cancer detection - Google Patents
Improved methods of cancer detection Download PDFInfo
- Publication number
- WO2014176340A1 WO2014176340A1 PCT/US2014/035153 US2014035153W WO2014176340A1 WO 2014176340 A1 WO2014176340 A1 WO 2014176340A1 US 2014035153 W US2014035153 W US 2014035153W WO 2014176340 A1 WO2014176340 A1 WO 2014176340A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- tissue
- region
- fractal dimension
- view
- fractal
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 165
- 206010028980 Neoplasm Diseases 0.000 title description 81
- 201000011510 cancer Diseases 0.000 title description 48
- 238000001514 detection method Methods 0.000 title description 11
- 210000001519 tissue Anatomy 0.000 claims description 312
- 210000000481 breast Anatomy 0.000 claims description 72
- 230000015654 memory Effects 0.000 claims description 37
- 238000004458 analytical method Methods 0.000 claims description 35
- 210000005013 brain tissue Anatomy 0.000 claims description 6
- 230000002500 effect on skin Effects 0.000 claims description 6
- 210000005228 liver tissue Anatomy 0.000 claims description 6
- 210000004072 lung Anatomy 0.000 claims description 6
- 210000004923 pancreatic tissue Anatomy 0.000 claims description 6
- 210000005084 renal tissue Anatomy 0.000 claims description 6
- 208000004434 Calcinosis Diseases 0.000 description 57
- 230000003211 malignant effect Effects 0.000 description 55
- 230000003902 lesion Effects 0.000 description 52
- 230000006870 function Effects 0.000 description 47
- 238000009826 distribution Methods 0.000 description 31
- 208000026310 Breast neoplasm Diseases 0.000 description 17
- 238000004891 communication Methods 0.000 description 17
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 17
- 238000001228 spectrum Methods 0.000 description 17
- 238000004195 computer-aided diagnosis Methods 0.000 description 16
- 206010006187 Breast cancer Diseases 0.000 description 15
- 238000011161 development Methods 0.000 description 14
- 230000018109 developmental process Effects 0.000 description 14
- 239000000523 sample Substances 0.000 description 14
- 238000012360 testing method Methods 0.000 description 13
- 230000008569 process Effects 0.000 description 12
- 230000000875 corresponding effect Effects 0.000 description 10
- 208000035475 disorder Diseases 0.000 description 10
- 238000005192 partition Methods 0.000 description 10
- 241001465754 Metazoa Species 0.000 description 9
- 238000012216 screening Methods 0.000 description 9
- 210000000577 adipose tissue Anatomy 0.000 description 8
- 230000006399 behavior Effects 0.000 description 8
- 208000003174 Brain Neoplasms Diseases 0.000 description 7
- 238000012512 characterization method Methods 0.000 description 7
- 201000010099 disease Diseases 0.000 description 7
- 230000008520 organization Effects 0.000 description 7
- 238000011282 treatment Methods 0.000 description 7
- 206010018338 Glioma Diseases 0.000 description 6
- 210000004027 cell Anatomy 0.000 description 6
- 208000032839 leukemia Diseases 0.000 description 6
- 239000000203 mixture Substances 0.000 description 6
- 206010039491 Sarcoma Diseases 0.000 description 5
- 238000004590 computer program Methods 0.000 description 5
- 238000003708 edge detection Methods 0.000 description 5
- 238000009499 grossing Methods 0.000 description 5
- 230000036210 malignancy Effects 0.000 description 5
- 238000009607 mammography Methods 0.000 description 5
- 210000000056 organ Anatomy 0.000 description 5
- 206010003571 Astrocytoma Diseases 0.000 description 4
- 208000032612 Glial tumor Diseases 0.000 description 4
- 206010025323 Lymphomas Diseases 0.000 description 4
- 238000000692 Student's t-test Methods 0.000 description 4
- 230000001154 acute effect Effects 0.000 description 4
- 230000002308 calcification Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 4
- 230000001413 cellular effect Effects 0.000 description 4
- 230000001684 chronic effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000002267 hypothalamic effect Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000002560 therapeutic procedure Methods 0.000 description 4
- 238000010207 Bayesian analysis Methods 0.000 description 3
- 208000008839 Kidney Neoplasms Diseases 0.000 description 3
- 208000031422 Lymphocytic Chronic B-Cell Leukemia Diseases 0.000 description 3
- 241000699666 Mus <mouse, genus> Species 0.000 description 3
- 206010035226 Plasma cell myeloma Diseases 0.000 description 3
- 206010038389 Renal cancer Diseases 0.000 description 3
- 238000011869 Shapiro-Wilk test Methods 0.000 description 3
- 238000011319 anticancer therapy Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 208000002458 carcinoid tumor Diseases 0.000 description 3
- 230000000295 complement effect Effects 0.000 description 3
- 238000003745 diagnosis Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000003814 drug Substances 0.000 description 3
- 230000001747 exhibiting effect Effects 0.000 description 3
- 239000000284 extract Substances 0.000 description 3
- 230000012010 growth Effects 0.000 description 3
- 230000003993 interaction Effects 0.000 description 3
- 238000011835 investigation Methods 0.000 description 3
- 230000001788 irregular Effects 0.000 description 3
- 201000010982 kidney cancer Diseases 0.000 description 3
- 230000007246 mechanism Effects 0.000 description 3
- 201000005962 mycosis fungoides Diseases 0.000 description 3
- 208000025113 myeloid leukemia Diseases 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 238000001959 radiotherapy Methods 0.000 description 3
- 239000013074 reference sample Substances 0.000 description 3
- 210000000130 stem cell Anatomy 0.000 description 3
- 201000008205 supratentorial primitive neuroectodermal tumor Diseases 0.000 description 3
- 238000012353 t test Methods 0.000 description 3
- 208000008732 thymoma Diseases 0.000 description 3
- 210000000239 visual pathway Anatomy 0.000 description 3
- 230000004400 visual pathway Effects 0.000 description 3
- 208000024893 Acute lymphoblastic leukemia Diseases 0.000 description 2
- 208000014697 Acute lymphocytic leukaemia Diseases 0.000 description 2
- 208000031261 Acute myeloid leukaemia Diseases 0.000 description 2
- 206010060971 Astrocytoma malignant Diseases 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 2
- 208000010839 B-cell chronic lymphocytic leukemia Diseases 0.000 description 2
- 208000032791 BCR-ABL1 positive chronic myelogenous leukemia Diseases 0.000 description 2
- 206010004593 Bile duct cancer Diseases 0.000 description 2
- 206010005949 Bone cancer Diseases 0.000 description 2
- 208000018084 Bone neoplasm Diseases 0.000 description 2
- 241000283690 Bos taurus Species 0.000 description 2
- 206010006143 Brain stem glioma Diseases 0.000 description 2
- 206010007275 Carcinoid tumour Diseases 0.000 description 2
- 206010007953 Central nervous system lymphoma Diseases 0.000 description 2
- 208000010833 Chronic myeloid leukaemia Diseases 0.000 description 2
- AOJJSUZBOXZQNB-TZSSRYMLSA-N Doxorubicin Chemical compound O([C@H]1C[C@@](O)(CC=2C(O)=C3C(=O)C=4C=CC=C(C=4C(=O)C3=C(O)C=21)OC)C(=O)CO)[C@H]1C[C@H](N)[C@H](O)[C@H](C)O1 AOJJSUZBOXZQNB-TZSSRYMLSA-N 0.000 description 2
- 206010014967 Ependymoma Diseases 0.000 description 2
- 208000012468 Ewing sarcoma/peripheral primitive neuroectodermal tumor Diseases 0.000 description 2
- 241000282326 Felis catus Species 0.000 description 2
- 208000021309 Germ cell tumor Diseases 0.000 description 2
- 206010061252 Intraocular melanoma Diseases 0.000 description 2
- 206010058467 Lung neoplasm malignant Diseases 0.000 description 2
- 206010025557 Malignant fibrous histiocytoma of bone Diseases 0.000 description 2
- 241000124008 Mammalia Species 0.000 description 2
- 238000000585 Mann–Whitney U test Methods 0.000 description 2
- 208000000172 Medulloblastoma Diseases 0.000 description 2
- 206010027406 Mesothelioma Diseases 0.000 description 2
- 208000003445 Mouth Neoplasms Diseases 0.000 description 2
- 208000034578 Multiple myelomas Diseases 0.000 description 2
- 208000033761 Myelogenous Chronic BCR-ABL Positive Leukemia Diseases 0.000 description 2
- 208000033776 Myeloid Acute Leukemia Diseases 0.000 description 2
- 208000034176 Neoplasms, Germ Cell and Embryonal Diseases 0.000 description 2
- 241000283973 Oryctolagus cuniculus Species 0.000 description 2
- 206010061902 Pancreatic neoplasm Diseases 0.000 description 2
- 241001494479 Pecora Species 0.000 description 2
- 241000009328 Perro Species 0.000 description 2
- 208000006664 Precursor Cell Lymphoblastic Leukemia-Lymphoma Diseases 0.000 description 2
- 241000288906 Primates Species 0.000 description 2
- 241000700159 Rattus Species 0.000 description 2
- 208000006265 Renal cell carcinoma Diseases 0.000 description 2
- 201000000582 Retinoblastoma Diseases 0.000 description 2
- 241000282898 Sus scrofa Species 0.000 description 2
- NKANXQFJJICGDU-QPLCGJKRSA-N Tamoxifen Chemical compound C=1C=CC=CC=1C(/CC)=C(C=1C=CC(OCCN(C)C)=CC=1)/C1=CC=CC=C1 NKANXQFJJICGDU-QPLCGJKRSA-N 0.000 description 2
- 208000024770 Thyroid neoplasm Diseases 0.000 description 2
- 108060008682 Tumor Necrosis Factor Proteins 0.000 description 2
- 201000005969 Uveal melanoma Diseases 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 239000002246 antineoplastic agent Substances 0.000 description 2
- 201000007335 cerebellar astrocytoma Diseases 0.000 description 2
- 208000030239 cerebral astrocytoma Diseases 0.000 description 2
- 230000002490 cerebral effect Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 208000032852 chronic lymphocytic leukemia Diseases 0.000 description 2
- 229940079593 drug Drugs 0.000 description 2
- 208000024519 eye neoplasm Diseases 0.000 description 2
- 230000002496 gastric effect Effects 0.000 description 2
- 208000029824 high grade glioma Diseases 0.000 description 2
- 210000004153 islets of langerhan Anatomy 0.000 description 2
- 210000000244 kidney pelvis Anatomy 0.000 description 2
- 230000000670 limiting effect Effects 0.000 description 2
- 208000012987 lip and oral cavity carcinoma Diseases 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 201000005202 lung cancer Diseases 0.000 description 2
- 208000020816 lung neoplasm Diseases 0.000 description 2
- 208000030883 malignant astrocytoma Diseases 0.000 description 2
- 201000011614 malignant glioma Diseases 0.000 description 2
- 208000015486 malignant pancreatic neoplasm Diseases 0.000 description 2
- 230000005741 malignant process Effects 0.000 description 2
- 201000001441 melanoma Diseases 0.000 description 2
- GLVAUDGFNGKCSF-UHFFFAOYSA-N mercaptopurine Chemical compound S=C1NC=NC2=C1NC=N2 GLVAUDGFNGKCSF-UHFFFAOYSA-N 0.000 description 2
- 208000018795 nasal cavity and paranasal sinus carcinoma Diseases 0.000 description 2
- 201000008106 ocular cancer Diseases 0.000 description 2
- 201000002575 ocular melanoma Diseases 0.000 description 2
- 201000008968 osteosarcoma Diseases 0.000 description 2
- 201000002528 pancreatic cancer Diseases 0.000 description 2
- 208000008443 pancreatic carcinoma Diseases 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 208000016800 primary central nervous system lymphoma Diseases 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000011218 segmentation Effects 0.000 description 2
- 201000000849 skin cancer Diseases 0.000 description 2
- 230000006886 spatial memory Effects 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000004083 survival effect Effects 0.000 description 2
- 201000002510 thyroid cancer Diseases 0.000 description 2
- 206010044412 transitional cell carcinoma Diseases 0.000 description 2
- 102000003390 tumor necrosis factor Human genes 0.000 description 2
- 208000018417 undifferentiated high grade pleomorphic sarcoma of bone Diseases 0.000 description 2
- 210000000626 ureter Anatomy 0.000 description 2
- PXFBZOLANLWPMH-UHFFFAOYSA-N 16-Epiaffinine Natural products C1C(C2=CC=CC=C2N2)=C2C(=O)CC2C(=CC)CN(C)C1C2CO PXFBZOLANLWPMH-UHFFFAOYSA-N 0.000 description 1
- 208000030507 AIDS Diseases 0.000 description 1
- 208000002008 AIDS-Related Lymphoma Diseases 0.000 description 1
- 241000251468 Actinopterygii Species 0.000 description 1
- 206010061424 Anal cancer Diseases 0.000 description 1
- 208000007860 Anus Neoplasms Diseases 0.000 description 1
- 206010073360 Appendix cancer Diseases 0.000 description 1
- 102000015790 Asparaginase Human genes 0.000 description 1
- 108010024976 Asparaginase Proteins 0.000 description 1
- 241000271566 Aves Species 0.000 description 1
- 206010004146 Basal cell carcinoma Diseases 0.000 description 1
- 206010005003 Bladder cancer Diseases 0.000 description 1
- 108010006654 Bleomycin Proteins 0.000 description 1
- 206010006298 Breast pain Diseases 0.000 description 1
- 230000005653 Brownian motion process Effects 0.000 description 1
- 208000011691 Burkitt lymphomas Diseases 0.000 description 1
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 190000008236 Carboplatin Chemical compound 0.000 description 1
- 206010007279 Carcinoid tumour of the gastrointestinal tract Diseases 0.000 description 1
- 206010007281 Carcinoid tumour of the stomach Diseases 0.000 description 1
- 201000009030 Carcinoma Diseases 0.000 description 1
- DLGOEMSEDOSKAD-UHFFFAOYSA-N Carmustine Chemical compound ClCCNC(=O)N(N=O)CCCl DLGOEMSEDOSKAD-UHFFFAOYSA-N 0.000 description 1
- 241000282693 Cercopithecidae Species 0.000 description 1
- 206010008342 Cervix carcinoma Diseases 0.000 description 1
- 208000032544 Cicatrix Diseases 0.000 description 1
- 206010009944 Colon cancer Diseases 0.000 description 1
- 241000938605 Crocodylia Species 0.000 description 1
- CMSMOCZEIVJLDB-UHFFFAOYSA-N Cyclophosphamide Chemical compound ClCCN(CCCl)P1(=O)NCCCO1 CMSMOCZEIVJLDB-UHFFFAOYSA-N 0.000 description 1
- -1 Cytarabile Chemical compound 0.000 description 1
- 208000008743 Desmoplastic Small Round Cell Tumor Diseases 0.000 description 1
- 206010064581 Desmoplastic small round cell tumour Diseases 0.000 description 1
- 206010014733 Endometrial cancer Diseases 0.000 description 1
- 206010014759 Endometrial neoplasm Diseases 0.000 description 1
- 102000004190 Enzymes Human genes 0.000 description 1
- 108090000790 Enzymes Proteins 0.000 description 1
- 241000283073 Equus caballus Species 0.000 description 1
- 208000000461 Esophageal Neoplasms Diseases 0.000 description 1
- 208000006168 Ewing Sarcoma Diseases 0.000 description 1
- 208000017259 Extragonadal germ cell tumor Diseases 0.000 description 1
- 208000000571 Fibrocystic breast disease Diseases 0.000 description 1
- GHASVSINZRGABV-UHFFFAOYSA-N Fluorouracil Chemical compound FC1=CNC(=O)NC1=O GHASVSINZRGABV-UHFFFAOYSA-N 0.000 description 1
- 208000022072 Gallbladder Neoplasms Diseases 0.000 description 1
- 241000238631 Hexapoda Species 0.000 description 1
- 208000017604 Hodgkin disease Diseases 0.000 description 1
- 208000021519 Hodgkin lymphoma Diseases 0.000 description 1
- 208000010747 Hodgkins lymphoma Diseases 0.000 description 1
- 241000282412 Homo Species 0.000 description 1
- 206010020843 Hyperthermia Diseases 0.000 description 1
- 206010021042 Hypopharyngeal cancer Diseases 0.000 description 1
- 206010056305 Hypopharyngeal neoplasm Diseases 0.000 description 1
- 102000014150 Interferons Human genes 0.000 description 1
- 108010050904 Interferons Proteins 0.000 description 1
- 102000015696 Interleukins Human genes 0.000 description 1
- 108010063738 Interleukins Proteins 0.000 description 1
- 208000007766 Kaposi sarcoma Diseases 0.000 description 1
- FBOZXECLQNJBKD-ZDUSSCGKSA-N L-methotrexate Chemical compound C=1N=C2N=C(N)N=C(N)C2=NC=1CN(C)C1=CC=C(C(=O)N[C@@H](CCC(O)=O)C(O)=O)C=C1 FBOZXECLQNJBKD-ZDUSSCGKSA-N 0.000 description 1
- 206010023825 Laryngeal cancer Diseases 0.000 description 1
- 108010000817 Leuprolide Proteins 0.000 description 1
- 206010061523 Lip and/or oral cavity cancer Diseases 0.000 description 1
- GQYIWUVLTXOXAJ-UHFFFAOYSA-N Lomustine Chemical compound ClCCN(N=O)C(=O)NC1CCCCC1 GQYIWUVLTXOXAJ-UHFFFAOYSA-N 0.000 description 1
- 206010025312 Lymphoma AIDS related Diseases 0.000 description 1
- 208000006644 Malignant Fibrous Histiocytoma Diseases 0.000 description 1
- 208000030070 Malignant epithelial tumor of ovary Diseases 0.000 description 1
- 206010073059 Malignant neoplasm of unknown primary site Diseases 0.000 description 1
- 208000032271 Malignant tumor of penis Diseases 0.000 description 1
- 208000006662 Mastodynia Diseases 0.000 description 1
- 208000002030 Merkel cell carcinoma Diseases 0.000 description 1
- 229930192392 Mitomycin Natural products 0.000 description 1
- 201000003793 Myelodysplastic syndrome Diseases 0.000 description 1
- 208000014767 Myeloproliferative disease Diseases 0.000 description 1
- 201000007224 Myeloproliferative neoplasm Diseases 0.000 description 1
- NWIBSHFKIJFRCO-WUDYKRTCSA-N Mytomycin Chemical compound C1N2C(C(C(C)=C(N)C3=O)=O)=C3[C@@H](COC(N)=O)[C@@]2(OC)[C@@H]2[C@H]1N2 NWIBSHFKIJFRCO-WUDYKRTCSA-N 0.000 description 1
- 208000002454 Nasopharyngeal Carcinoma Diseases 0.000 description 1
- 206010061306 Nasopharyngeal cancer Diseases 0.000 description 1
- 206010029260 Neuroblastoma Diseases 0.000 description 1
- 206010029266 Neuroendocrine carcinoma of the skin Diseases 0.000 description 1
- 208000015914 Non-Hodgkin lymphomas Diseases 0.000 description 1
- 208000009893 Nonpenetrating Wounds Diseases 0.000 description 1
- 206010030155 Oesophageal carcinoma Diseases 0.000 description 1
- 206010031096 Oropharyngeal cancer Diseases 0.000 description 1
- 206010057444 Oropharyngeal neoplasm Diseases 0.000 description 1
- 208000007571 Ovarian Epithelial Carcinoma Diseases 0.000 description 1
- 206010033128 Ovarian cancer Diseases 0.000 description 1
- 206010061328 Ovarian epithelial cancer Diseases 0.000 description 1
- 206010033268 Ovarian low malignant potential tumour Diseases 0.000 description 1
- 206010061535 Ovarian neoplasm Diseases 0.000 description 1
- 229930012538 Paclitaxel Natural products 0.000 description 1
- 208000000821 Parathyroid Neoplasms Diseases 0.000 description 1
- 208000002471 Penile Neoplasms Diseases 0.000 description 1
- 206010034299 Penile cancer Diseases 0.000 description 1
- 208000009565 Pharyngeal Neoplasms Diseases 0.000 description 1
- 206010034811 Pharyngeal cancer Diseases 0.000 description 1
- 206010035052 Pineal germinoma Diseases 0.000 description 1
- 206010050487 Pinealoblastoma Diseases 0.000 description 1
- 208000007641 Pinealoma Diseases 0.000 description 1
- 208000007913 Pituitary Neoplasms Diseases 0.000 description 1
- 201000005746 Pituitary adenoma Diseases 0.000 description 1
- 206010061538 Pituitary tumour benign Diseases 0.000 description 1
- 201000008199 Pleuropulmonary blastoma Diseases 0.000 description 1
- 208000006994 Precancerous Conditions Diseases 0.000 description 1
- 206010060862 Prostate cancer Diseases 0.000 description 1
- 208000000236 Prostatic Neoplasms Diseases 0.000 description 1
- 208000015634 Rectal Neoplasms Diseases 0.000 description 1
- 241000283984 Rodentia Species 0.000 description 1
- 208000004337 Salivary Gland Neoplasms Diseases 0.000 description 1
- 206010061934 Salivary gland cancer Diseases 0.000 description 1
- 208000009359 Sezary Syndrome Diseases 0.000 description 1
- 208000021388 Sezary disease Diseases 0.000 description 1
- 208000000453 Skin Neoplasms Diseases 0.000 description 1
- 208000021712 Soft tissue sarcoma Diseases 0.000 description 1
- 231100000632 Spindle poison Toxicity 0.000 description 1
- 208000005718 Stomach Neoplasms Diseases 0.000 description 1
- 208000031673 T-Cell Cutaneous Lymphoma Diseases 0.000 description 1
- 206010042971 T-cell lymphoma Diseases 0.000 description 1
- 208000027585 T-cell non-Hodgkin lymphoma Diseases 0.000 description 1
- 210000001744 T-lymphocyte Anatomy 0.000 description 1
- 208000024313 Testicular Neoplasms Diseases 0.000 description 1
- 206010057644 Testis cancer Diseases 0.000 description 1
- 206010043515 Throat cancer Diseases 0.000 description 1
- 201000009365 Thymic carcinoma Diseases 0.000 description 1
- 206010044407 Transitional cell cancer of the renal pelvis and ureter Diseases 0.000 description 1
- 208000015778 Undifferentiated pleomorphic sarcoma Diseases 0.000 description 1
- 206010046431 Urethral cancer Diseases 0.000 description 1
- 206010046458 Urethral neoplasms Diseases 0.000 description 1
- 208000007097 Urinary Bladder Neoplasms Diseases 0.000 description 1
- 208000006105 Uterine Cervical Neoplasms Diseases 0.000 description 1
- 208000002495 Uterine Neoplasms Diseases 0.000 description 1
- JXLYSJRDGCGARV-WWYNWVTFSA-N Vinblastine Natural products O=C(O[C@H]1[C@](O)(C(=O)OC)[C@@H]2N(C)c3c(cc(c(OC)c3)[C@]3(C(=O)OC)c4[nH]c5c(c4CCN4C[C@](O)(CC)C[C@H](C3)C4)cccc5)[C@@]32[C@H]2[C@@]1(CC)C=CCN2CC3)C JXLYSJRDGCGARV-WWYNWVTFSA-N 0.000 description 1
- 206010047741 Vulval cancer Diseases 0.000 description 1
- 208000004354 Vulvar Neoplasms Diseases 0.000 description 1
- 208000033559 Waldenström macroglobulinemia Diseases 0.000 description 1
- 208000008383 Wilms tumor Diseases 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 208000020990 adrenal cortex carcinoma Diseases 0.000 description 1
- 208000007128 adrenocortical carcinoma Diseases 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000002152 alkylating effect Effects 0.000 description 1
- 239000003242 anti bacterial agent Substances 0.000 description 1
- 230000003474 anti-emetic effect Effects 0.000 description 1
- 230000000340 anti-metabolite Effects 0.000 description 1
- 229940088710 antibiotic agent Drugs 0.000 description 1
- 229940125683 antiemetic agent Drugs 0.000 description 1
- 239000002111 antiemetic agent Substances 0.000 description 1
- 229940100197 antimetabolite Drugs 0.000 description 1
- 239000002256 antimetabolite Substances 0.000 description 1
- 229940045719 antineoplastic alkylating agent nitrosoureas Drugs 0.000 description 1
- 201000011165 anus cancer Diseases 0.000 description 1
- 208000021780 appendiceal neoplasm Diseases 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 229960003272 asparaginase Drugs 0.000 description 1
- DCXYFEDJOCDNAF-UHFFFAOYSA-M asparaginate Chemical compound [O-]C(=O)C(N)CC(N)=O DCXYFEDJOCDNAF-UHFFFAOYSA-M 0.000 description 1
- 238000011888 autopsy Methods 0.000 description 1
- 238000013477 bayesian statistics method Methods 0.000 description 1
- 208000026900 bile duct neoplasm Diseases 0.000 description 1
- 230000002902 bimodal effect Effects 0.000 description 1
- 230000008512 biological response Effects 0.000 description 1
- 229960001561 bleomycin Drugs 0.000 description 1
- OYVAGSVQBOHSSS-UAPAGMARSA-O bleomycin A2 Chemical compound N([C@H](C(=O)N[C@H](C)[C@@H](O)[C@H](C)C(=O)N[C@@H]([C@H](O)C)C(=O)NCCC=1SC=C(N=1)C=1SC=C(N=1)C(=O)NCCC[S+](C)C)[C@@H](O[C@H]1[C@H]([C@@H](O)[C@H](O)[C@H](CO)O1)O[C@@H]1[C@H]([C@@H](OC(N)=O)[C@H](O)[C@@H](CO)O1)O)C=1N=CNC=1)C(=O)C1=NC([C@H](CC(N)=O)NC[C@H](N)C(N)=O)=NC(N)=C1C OYVAGSVQBOHSSS-UAPAGMARSA-O 0.000 description 1
- 210000001185 bone marrow Anatomy 0.000 description 1
- 201000008873 bone osteosarcoma Diseases 0.000 description 1
- 208000012172 borderline epithelial tumor of ovary Diseases 0.000 description 1
- 238000002725 brachytherapy Methods 0.000 description 1
- 208000011803 breast fibrocystic disease Diseases 0.000 description 1
- 201000002143 bronchus adenoma Diseases 0.000 description 1
- 238000005537 brownian motion Methods 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000005773 cancer-related death Effects 0.000 description 1
- 229960004562 carboplatin Drugs 0.000 description 1
- 231100000504 carcinogenesis Toxicity 0.000 description 1
- 229960005243 carmustine Drugs 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000010261 cell growth Effects 0.000 description 1
- 230000007248 cellular mechanism Effects 0.000 description 1
- 210000003169 central nervous system Anatomy 0.000 description 1
- 201000010881 cervical cancer Diseases 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 229940044683 chemotherapy drug Drugs 0.000 description 1
- 201000008522 childhood cerebral astrocytoma Diseases 0.000 description 1
- 208000011654 childhood malignant neoplasm Diseases 0.000 description 1
- 229960004630 chlorambucil Drugs 0.000 description 1
- JCKYGMPEJWAADB-UHFFFAOYSA-N chlorambucil Chemical compound OC(=O)CCCC1=CC=C(N(CCCl)CCCl)C=C1 JCKYGMPEJWAADB-UHFFFAOYSA-N 0.000 description 1
- 208000006990 cholangiocarcinoma Diseases 0.000 description 1
- 229960004316 cisplatin Drugs 0.000 description 1
- DQLATGHUWYMOKM-UHFFFAOYSA-L cisplatin Chemical compound N[Pt](N)(Cl)Cl DQLATGHUWYMOKM-UHFFFAOYSA-L 0.000 description 1
- 238000007621 cluster analysis Methods 0.000 description 1
- 208000029742 colonic neoplasm Diseases 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000004883 computer application Methods 0.000 description 1
- 230000001276 controlling effect Effects 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000315 cryotherapy Methods 0.000 description 1
- 201000007241 cutaneous T cell lymphoma Diseases 0.000 description 1
- 208000017763 cutaneous neuroendocrine carcinoma Diseases 0.000 description 1
- 229960004397 cyclophosphamide Drugs 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 239000002781 deodorant agent Substances 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 229960004679 doxorubicin Drugs 0.000 description 1
- 238000010894 electron beam technology Methods 0.000 description 1
- 238000009261 endocrine therapy Methods 0.000 description 1
- 229940034984 endocrine therapy antineoplastic and immunomodulating agent Drugs 0.000 description 1
- 230000002357 endometrial effect Effects 0.000 description 1
- 229940088598 enzyme Drugs 0.000 description 1
- 201000004101 esophageal cancer Diseases 0.000 description 1
- VJJPUSNTGOMMGY-MRVIYFEKSA-N etoposide Chemical compound COC1=C(O)C(OC)=CC([C@@H]2C3=CC=4OCOC=4C=C3[C@@H](O[C@H]3[C@@H]([C@@H](O)[C@@H]4O[C@H](C)OC[C@H]4O3)O)[C@@H]3[C@@H]2C(OC3)=O)=C1 VJJPUSNTGOMMGY-MRVIYFEKSA-N 0.000 description 1
- 229960005420 etoposide Drugs 0.000 description 1
- 201000008819 extrahepatic bile duct carcinoma Diseases 0.000 description 1
- 238000000799 fluorescence microscopy Methods 0.000 description 1
- 229960002949 fluorouracil Drugs 0.000 description 1
- 229960002074 flutamide Drugs 0.000 description 1
- MKXKFYHWDHIYRV-UHFFFAOYSA-N flutamide Chemical compound CC(C)C(=O)NC1=CC=C([N+]([O-])=O)C(C(F)(F)F)=C1 MKXKFYHWDHIYRV-UHFFFAOYSA-N 0.000 description 1
- 201000010175 gallbladder cancer Diseases 0.000 description 1
- 206010017758 gastric cancer Diseases 0.000 description 1
- 201000011243 gastrointestinal stromal tumor Diseases 0.000 description 1
- 229960005277 gemcitabine Drugs 0.000 description 1
- SDUQYLNIPVEERB-QPPQHZFASA-N gemcitabine Chemical compound O=C1N=C(N)C=CN1[C@H]1C(F)(F)[C@H](O)[C@@H](CO)O1 SDUQYLNIPVEERB-QPPQHZFASA-N 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 201000007116 gestational trophoblastic neoplasm Diseases 0.000 description 1
- 208000035474 group of disease Diseases 0.000 description 1
- 201000009277 hairy cell leukemia Diseases 0.000 description 1
- 201000010536 head and neck cancer Diseases 0.000 description 1
- 208000014829 head and neck neoplasm Diseases 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 201000010235 heart cancer Diseases 0.000 description 1
- 208000024348 heart neoplasm Diseases 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 230000036031 hyperthermia Effects 0.000 description 1
- 238000009217 hyperthermia therapy Methods 0.000 description 1
- 201000006866 hypopharynx cancer Diseases 0.000 description 1
- 229960001101 ifosfamide Drugs 0.000 description 1
- HOMGKSMUEGBAAB-UHFFFAOYSA-N ifosfamide Chemical compound ClCCNP1(=O)OCCCN1CCCl HOMGKSMUEGBAAB-UHFFFAOYSA-N 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000007943 implant Substances 0.000 description 1
- 229910001410 inorganic ion Inorganic materials 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 229940047124 interferons Drugs 0.000 description 1
- 229940047122 interleukins Drugs 0.000 description 1
- 229960004768 irinotecan Drugs 0.000 description 1
- UWKQSNNFCGGAFS-XIFFEERXSA-N irinotecan Chemical compound C1=C2C(CC)=C3CN(C(C4=C([C@@](C(=O)OC4)(O)CC)C=4)=O)C=4C3=NC2=CC=C1OC(=O)N(CC1)CCC1N1CCCCC1 UWKQSNNFCGGAFS-XIFFEERXSA-N 0.000 description 1
- 206010023841 laryngeal neoplasm Diseases 0.000 description 1
- GFIJNRVAKGFPGQ-LIJARHBVSA-N leuprolide Chemical compound CCNC(=O)[C@@H]1CCCN1C(=O)[C@H](CCCNC(N)=N)NC(=O)[C@H](CC(C)C)NC(=O)[C@@H](CC(C)C)NC(=O)[C@@H](NC(=O)[C@H](CO)NC(=O)[C@H](CC=1C2=CC=CC=C2NC=1)NC(=O)[C@H](CC=1N=CNC=1)NC(=O)[C@H]1NC(=O)CC1)CC1=CC=C(O)C=C1 GFIJNRVAKGFPGQ-LIJARHBVSA-N 0.000 description 1
- 229960004338 leuprorelin Drugs 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 206010024627 liposarcoma Diseases 0.000 description 1
- 210000004185 liver Anatomy 0.000 description 1
- 201000007270 liver cancer Diseases 0.000 description 1
- 208000014018 liver neoplasm Diseases 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000011551 log transformation method Methods 0.000 description 1
- 229960002247 lomustine Drugs 0.000 description 1
- 230000000527 lymphocytic effect Effects 0.000 description 1
- 201000000564 macroglobulinemia Diseases 0.000 description 1
- 238000002595 magnetic resonance imaging Methods 0.000 description 1
- 208000026045 malignant tumor of parathyroid gland Diseases 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 229960004961 mechlorethamine Drugs 0.000 description 1
- HAWPXGHAZFHHAD-UHFFFAOYSA-N mechlorethamine Chemical compound ClCCN(C)CCCl HAWPXGHAZFHHAD-UHFFFAOYSA-N 0.000 description 1
- 229960001786 megestrol Drugs 0.000 description 1
- RQZAXGRLVPAYTJ-GQFGMJRRSA-N megestrol acetate Chemical compound C1=C(C)C2=CC(=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@@](C(C)=O)(OC(=O)C)[C@@]1(C)CC2 RQZAXGRLVPAYTJ-GQFGMJRRSA-N 0.000 description 1
- 229960001924 melphalan Drugs 0.000 description 1
- SGDBTWWWUNNDEQ-LBPRGKRZSA-N melphalan Chemical compound OC(=O)[C@@H](N)CC1=CC=C(N(CCCl)CCCl)C=C1 SGDBTWWWUNNDEQ-LBPRGKRZSA-N 0.000 description 1
- 229960001428 mercaptopurine Drugs 0.000 description 1
- 210000000716 merkel cell Anatomy 0.000 description 1
- 230000001394 metastastic effect Effects 0.000 description 1
- 206010061289 metastatic neoplasm Diseases 0.000 description 1
- 208000037970 metastatic squamous neck cancer Diseases 0.000 description 1
- 229960000485 methotrexate Drugs 0.000 description 1
- 229960004857 mitomycin Drugs 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 206010051747 multiple endocrine neoplasia Diseases 0.000 description 1
- 201000006938 muscular dystrophy Diseases 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 208000017869 myelodysplastic/myeloproliferative disease Diseases 0.000 description 1
- 201000000050 myeloid neoplasm Diseases 0.000 description 1
- 201000011216 nasopharynx carcinoma Diseases 0.000 description 1
- 230000009826 neoplastic cell growth Effects 0.000 description 1
- 201000008026 nephroblastoma Diseases 0.000 description 1
- 230000006855 networking Effects 0.000 description 1
- 210000002445 nipple Anatomy 0.000 description 1
- 208000002154 non-small cell lung carcinoma Diseases 0.000 description 1
- 238000011275 oncology therapy Methods 0.000 description 1
- 208000022982 optic pathway glioma Diseases 0.000 description 1
- 201000006958 oropharynx cancer Diseases 0.000 description 1
- 230000000399 orthopedic effect Effects 0.000 description 1
- 230000002611 ovarian Effects 0.000 description 1
- 208000021284 ovarian germ cell tumor Diseases 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 229960001592 paclitaxel Drugs 0.000 description 1
- 201000002530 pancreatic endocrine carcinoma Diseases 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 238000003909 pattern recognition Methods 0.000 description 1
- 230000002085 persistent effect Effects 0.000 description 1
- 208000028591 pheochromocytoma Diseases 0.000 description 1
- 201000007315 pineal gland astrocytoma Diseases 0.000 description 1
- 201000004838 pineal region germinoma Diseases 0.000 description 1
- 201000003113 pineoblastoma Diseases 0.000 description 1
- 208000021310 pituitary gland adenoma Diseases 0.000 description 1
- 208000010626 plasma cell neoplasm Diseases 0.000 description 1
- 210000004180 plasmocyte Anatomy 0.000 description 1
- 208000025638 primary cutaneous T-cell non-Hodgkin lymphoma Diseases 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 238000002661 proton therapy Methods 0.000 description 1
- 239000000649 purine antagonist Substances 0.000 description 1
- 239000003790 pyrimidine antagonist Substances 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000002285 radioactive effect Effects 0.000 description 1
- 206010038038 rectal cancer Diseases 0.000 description 1
- 201000001275 rectum cancer Diseases 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 208000015347 renal cell adenocarcinoma Diseases 0.000 description 1
- 208000030859 renal pelvis/ureter urothelial carcinoma Diseases 0.000 description 1
- 201000009410 rhabdomyosarcoma Diseases 0.000 description 1
- 231100000241 scar Toxicity 0.000 description 1
- 230000037387 scars Effects 0.000 description 1
- 238000005204 segregation Methods 0.000 description 1
- 230000001953 sensory effect Effects 0.000 description 1
- 210000003491 skin Anatomy 0.000 description 1
- 201000008261 skin carcinoma Diseases 0.000 description 1
- 201000002314 small intestine cancer Diseases 0.000 description 1
- 210000004872 soft tissue Anatomy 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 206010041823 squamous cell carcinoma Diseases 0.000 description 1
- 208000037969 squamous neck cancer Diseases 0.000 description 1
- 238000000551 statistical hypothesis test Methods 0.000 description 1
- 238000000528 statistical test Methods 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 201000011549 stomach cancer Diseases 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 208000011580 syndromic disease Diseases 0.000 description 1
- 230000009885 systemic effect Effects 0.000 description 1
- 229960001603 tamoxifen Drugs 0.000 description 1
- RCINICONZNJXQF-MZXODVADSA-N taxol Chemical compound O([C@@H]1[C@@]2(C[C@@H](C(C)=C(C2(C)C)[C@H](C([C@]2(C)[C@@H](O)C[C@H]3OC[C@]3([C@H]21)OC(C)=O)=O)OC(=O)C)OC(=O)[C@H](O)[C@@H](NC(=O)C=1C=CC=CC=1)C=1C=CC=CC=1)O)C(=O)C1=CC=CC=C1 RCINICONZNJXQF-MZXODVADSA-N 0.000 description 1
- 201000003120 testicular cancer Diseases 0.000 description 1
- 229960000303 topotecan Drugs 0.000 description 1
- UCFGDBYHRUNTLO-QHCPKHFHSA-N topotecan Chemical compound C1=C(O)C(CN(C)C)=C2C=C(CN3C4=CC5=C(C3=O)COC(=O)[C@]5(O)CC)C4=NC2=C1 UCFGDBYHRUNTLO-QHCPKHFHSA-N 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000009261 transgenic effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
- 208000029387 trophoblastic neoplasm Diseases 0.000 description 1
- 239000000717 tumor promoter Substances 0.000 description 1
- 208000029729 tumor suppressor gene on chromosome 11 Diseases 0.000 description 1
- 230000002100 tumorsuppressive effect Effects 0.000 description 1
- 201000005112 urinary bladder cancer Diseases 0.000 description 1
- 206010046766 uterine cancer Diseases 0.000 description 1
- 208000037965 uterine sarcoma Diseases 0.000 description 1
- 206010046885 vaginal cancer Diseases 0.000 description 1
- 208000013139 vaginal neoplasm Diseases 0.000 description 1
- 229960003048 vinblastine Drugs 0.000 description 1
- JXLYSJRDGCGARV-XQKSVPLYSA-N vincaleukoblastine Chemical compound C([C@@H](C[C@]1(C(=O)OC)C=2C(=CC3=C([C@]45[C@H]([C@@]([C@H](OC(C)=O)[C@]6(CC)C=CCN([C@H]56)CC4)(O)C(=O)OC)N3C)C=2)OC)C[C@@](C2)(O)CC)N2CCC2=C1NC1=CC=CC=C21 JXLYSJRDGCGARV-XQKSVPLYSA-N 0.000 description 1
- 229960004528 vincristine Drugs 0.000 description 1
- OGWKCGZFUXNPDA-XQKSVPLYSA-N vincristine Chemical compound C([N@]1C[C@@H](C[C@]2(C(=O)OC)C=3C(=CC4=C([C@]56[C@H]([C@@]([C@H](OC(C)=O)[C@]7(CC)C=CCN([C@H]67)CC5)(O)C(=O)OC)N4C=O)C=3)OC)C[C@@](C1)(O)CC)CC1=C2NC2=CC=CC=C12 OGWKCGZFUXNPDA-XQKSVPLYSA-N 0.000 description 1
- OGWKCGZFUXNPDA-UHFFFAOYSA-N vincristine Natural products C1C(CC)(O)CC(CC2(C(=O)OC)C=3C(=CC4=C(C56C(C(C(OC(C)=O)C7(CC)C=CCN(C67)CC5)(O)C(=O)OC)N4C=O)C=3)OC)CN1CCC1=C2NC2=CC=CC=C12 OGWKCGZFUXNPDA-UHFFFAOYSA-N 0.000 description 1
- GBABOYUKABKIAF-GHYRFKGUSA-N vinorelbine Chemical compound C1N(CC=2C3=CC=CC=C3NC=22)CC(CC)=C[C@H]1C[C@]2(C(=O)OC)C1=CC([C@]23[C@H]([C@]([C@H](OC(C)=O)[C@]4(CC)C=CCN([C@H]34)CC2)(O)C(=O)OC)N2C)=C2C=C1OC GBABOYUKABKIAF-GHYRFKGUSA-N 0.000 description 1
- 229960002066 vinorelbine Drugs 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
- 201000005102 vulva cancer Diseases 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4312—Breast evaluation or disorder diagnosis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/7257—Details of waveform analysis characterised by using transforms using Fourier transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7235—Details of waveform analysis
- A61B5/7253—Details of waveform analysis characterised by using transforms
- A61B5/726—Details of waveform analysis characterised by using transforms using Wavelet transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/502—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of breast, i.e. mammography
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/40—Analysis of texture
- G06T7/41—Analysis of texture based on statistical description of texture
- G06T7/48—Analysis of texture based on statistical description of texture using fractals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2576/00—Medical imaging apparatus involving image processing or analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/08—Detecting, measuring or recording devices for evaluating the respiratory organs
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/20—Measuring for diagnostic purposes; Identification of persons for measuring urological functions restricted to the evaluation of the urinary system
- A61B5/201—Assessing renal or kidney functions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/40—Detecting, measuring or recording for evaluating the nervous system
- A61B5/4058—Detecting, measuring or recording for evaluating the nervous system for evaluating the central nervous system
- A61B5/4064—Evaluating the brain
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
- A61B5/4222—Evaluating particular parts, e.g. particular organs
- A61B5/4244—Evaluating particular parts, e.g. particular organs liver
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/42—Detecting, measuring or recording for evaluating the gastrointestinal, the endocrine or the exocrine systems
- A61B5/4222—Evaluating particular parts, e.g. particular organs
- A61B5/425—Evaluating particular parts, e.g. particular organs pancreas
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/43—Detecting, measuring or recording for evaluating the reproductive systems
- A61B5/4306—Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
- A61B5/4318—Evaluation of the lower reproductive system
- A61B5/4325—Evaluation of the lower reproductive system of the uterine cavities, e.g. uterus, fallopian tubes, ovaries
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/44—Detecting, measuring or recording for evaluating the integumentary system, e.g. skin, hair or nails
- A61B5/441—Skin evaluation, e.g. for skin disorder diagnosis
- A61B5/444—Evaluating skin marks, e.g. mole, nevi, tumour, scar
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20064—Wavelet transform [DWT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30024—Cell structures in vitro; Tissue sections in vitro
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H30/00—ICT specially adapted for the handling or processing of medical images
- G16H30/40—ICT specially adapted for the handling or processing of medical images for processing medical images, e.g. editing
Definitions
- the human body is composed of about ten trillion cells, and throughout a lifetime cells go through many damaging processes that concern their genetic programming. However, not all of these mutations result in cancer. As it turns out, many individuals actually have occult tumors throughout their bodies that are only discovered through microscopic investigation during autopsy. Therefore, there must be some mechanism that prohibits the growth and development of tumors.
- Recent research has pointed to the microenvironment of potential tumors to help suppress malignant phenotype and instruct otherwise malignant cells to participate in normal development.
- the microenvironment may actually provide tumor suppressive signals as long as the tissue architecture is controlled, but as the structure of healthy tissue is lost, the tissue has the potential to become a tumor promoter.
- the tissue surrounding the tumor, along with the lesion may be important in identifying potential cancerous tissue.
- the present invention is based, in part, on the realization that disorder of the microenvironment surrounding a portion of tissue in a subject may be indicative of the presence of, or pre-disposition toward developing, cancer in that tissue.
- the present invention provides, for the first time, methods for analyzing the amount and distribution of disorder across tissue in order to predict the presence of, or propensity of tissue to develop, cancer.
- provided methods include analysis of images of target tissue, as opposed to analysis of the tissue itself.
- provided methods include analysis of tissue or images of tissue adjacent to a target tissue/tissue region of interest, as opposed to the target tissue itself.
- the present invention provides, among other things, methods of identifying cancerous or pre-cancerous tissue including providing a first region of tissue from a subject, calculating a roughness exponent (e.g., a Hurst exponent) for the first region of tissue, and comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 between the roughness exponent of the first region of tissue and 0.5 indicates that the tissue is cancerous or pre-cancerous.
- the providing, calculating, and comparing steps are each performed a plurality of times. In some embodiments, the plurality of times is at least 10, at least 100, or at least 1,000. In some embodiments, performing the providing, calculating and comparing steps a plurality of times allows for a higher resolution characterization of the morphology of the tissue surrounding a region of tissue suspected to be cancerous or pre-cancerous.
- 0.5 may be less than 0.2. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.15. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.1. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.05.
- a roughness exponent is calculated from one or more images of a target tissue. In some embodiments, a roughness exponent is calculated from two or more images of a first region of tissue. In some embodiments, each of the two or more images are taken from different angles and/or points of view. In some embodiments, a roughness exponent is calculated for two or more regions of tissue. In some embodiments, the two or more regions of tissue are from a single organ. In some embodiments, the two or more regions of tissue are from different organs.
- a roughness exponent may be generated via any appropriate multi-scale analytical method.
- the roughness exponent is calculated using one or more multi-scale methods selected from a wavelet-transform modulus maxima, a wavelet leader, detrended fluctuation, and Fourier analysis.
- a roughness exponent may also be calculated using any mathematically related and/or similar multi-scale density fluctuation assessment method yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent.
- a roughness exponent is calculated from two- dimensional data.
- a roughness exponent is calculated from three dimensional data (e.g., 3D data cubes).
- the present invention additionally provides methods including providing a first view of a region of tissue; providing a second view of the region of tissue; calculating a first fractal dimension for the first view of the region of tissue; and calculating a second fractal dimension for the second view of the region of tissue; wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
- provided methods further comprise treating the region of tissue if it is cancerous.
- the fractal zone is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
- the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8); first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8); second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2); third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
- the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9); first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9); second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1); third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
- the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7); first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7); second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3); third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
- provided methods further include steps of providing a third view of the region of tissue; and calculating a third fractal dimension for the third view of the region of tissue; wherein if the fractal dimension of at least one of the first fractal dimension, the second fractal dimension, and the third fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
- a fractal dimension is calculated from two-dimensional data.
- a fractal dimension is calculated from three dimensional data (e.g., 3D data cubes).
- the methods provided by the present invention may be performed on any tissue.
- the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
- animal refers to any member of the animal kingdom. In some embodiments, “animal” refers to humans, at any stage of development. In some embodiments, “animal” refers to non-human animals, at any stage of development. In certain embodiments, the non-human animal is a mammal ⁇ e.g., a rodent, a mouse, a rat, a rabbit, a monkey, a dog, a cat, a sheep, cattle, a primate, and/or a pig). In some embodiments, animals include, but are not limited to, mammals, birds, reptiles, amphibians, fish, insects, and/or worms. In some embodiments, an animal may be a transgenic animal, genetically-engineered animal, and/or a clone.
- cancer refers to a group of diseases, all involving unregulated cell growth.
- exemplary cancers include, without limitation: Acute lymphoblastic leukemia, Acute myeloid leukemia, Adrenocortical carcinoma; AIDS-related cancers; AIDS-related lymphoma; Anal cancer; Appendix cancer; Astrocytoma,childhood cerebellar or cerebral; Basal cell carcinoma; Bile duct cancer, extrahepatic; Bladder cancer; Bone cancer, Osteosarcoma/Malignant fibrous histiocytoma; Brainstem glioma; Brain tumor; Brain tumor, cerebellar astrocytoma; Brain tumor, cerebral astrocytoma/malignant glioma; Brain tumor, ependymoma; Brain tumor, medulloblastoma; Brain tumor, supratentorial primitive neuroectodermal tumors; Brain tumor, visual pathway and hypothalamic glioma;
- Extragonadal Germ cell tumor Extrahepatic bile duct cancer
- Eye Cancer Intraocular melanoma
- Eye Cancer Retinoblastoma
- Gallbladder cancer Gastric (Stomach) cancer
- Gastrointestinal Carcinoid Tumor Gastrointestinal stromal tumor (GIST); Germ cell tumor: extracranial, extragonadal, or ovarian; Gestational trophoblastic tumor; Glioma of the brain stem; Glioma, Childhood Cerebral Astrocytoma; Glioma, Childhood Visual Pathway and
- hypothalamic Gastric carcinoid; Hairy cell leukemia; Head and neck cancer; Heart cancer; Hepatocellular (liver) cancer; Hodgkin lymphoma; Hypopharyngeal cancer; Hypothalamic and visual pathway glioma, childhood; Intraocular Melanoma; Islet Cell Carcinoma (Endocrine Pancreas); Kaposi sarcoma; Kidney cancer (renal cell cancer); Laryngeal Cancer; Leukemias; Leukemia, acute lymphoblastic (also called acute lymphocytic leukemia); Leukemia, acute myeloid (also called acute myelogenous leukemia); Leukemia, chronic lymphocytic (also called chronic lymphocytic leukemia); Leukemia, chronic myelogenous (also called chronic myeloid leukemia); Leukemia, hairy cell; Lip and Oral Cavity Cancer; Liposarcoma; Liver Cancer (Primary); Lung Cancer, Non-Small Cell; Lung Cancer Small
- Myelodysplasia Syndromes Myelodysplastic/Myeloproliferative Diseases; Myelogenous Leukemia, Chronic; Myeloid Leukemia, Adult Acute; Myeloid Leukemia, Childhood Acute; Myeloma, Multiple (Cancer of the Bone-Marrow); Myeloproliferative Disorders, Chronic; Nasal cavity and paranasal sinus cancer; Nasopharyngeal carcinoma; Neuroblastoma; Non-Hodgkin lymphoma; Non-small cell lung cancer; Oral Cancer; Oropharyngeal cancer;
- Sarcoma Kaposi; Sarcoma, soft tissue; Sarcoma, uterine; Sezary syndrome; Skin cancer (nonmelanoma); Skin carcinoma, Merkel cell; Small intestine cancer; Soft tissue sarcoma;
- Pre-Cancer As used herein the term “pre-cancer” or “pre-cancerous” is used to refer to tissue in or from a subject that is not yet cancerous, but has a higher chance of becoming cancerous that normal tissue.
- a "risk" of a disease, disorder, and/or condition comprises a likelihood that a particular individual will develop a disease, disorder, and/or condition ⁇ e.g., a cancer).
- risk is expressed as a percentage.
- risk is from 0,1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90 up to 100%.
- risk is expressed as a risk relative to a risk associated with a reference sample or group of reference samples.
- a reference sample or group of reference samples have a known risk of a disease, disorder, condition and/or event (e.g., a muscular dystrophy).
- a reference sample or group of reference samples are from individuals comparable to a particular individual.
- relative risk is 0,1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more.
- Subject refers to a human or any non-human animal (e.g., mouse, rat, rabbit, dog, cat, cattle, swine, sheep, horse or primate).
- a human includes pre and post natal forms.
- a subject is a human being.
- a subject can be a patient, which refers to a human presenting to a medical provider for diagnosis or treatment of a disease.
- the term "subject” is used herein interchangeably with “individual” or "patient.”
- a subject can be afflicted with or is susceptible to a disease or disorder but may or may not display symptoms of the disease or disorder.
- the term “substantially” refers to the qualitative condition of exhibiting total or near-total extent or degree of a characteristic or property of interest.
- One of ordinary skill in the biological arts will understand that biological and chemical phenomena rarely, if ever, go to completion and/or proceed to completeness or achieve or avoid an absolute result.
- the term “substantially” is therefore used herein to capture the potential lack of completeness inherent in many biological and chemical phenomena.
- FIG. 1 shows a representative schematic of imaging breast tissue from
- MLO mediolateral oblique
- CC cranio-caudal
- FIG. 2 depicts a graph of the minimum fractal dimension by maximum fractal dimension, with tissue classified according to known methods as benign in blue, and tissue classified according to known methods as malignant tissue, in red.
- FIG. 3 shows exemplary results of a cluster analysis of scatter plots obtained from graphing minimum fractal dimension by maximum fractal dimension of (a) benign tissue, and (b) malignant tissue.
- FIG. 4 shows images of breast tissue classified as (a) benign and (b) malignant as classified by provided methods. Note the presence of a majority of disrupted tissue in the image of malignant tissue (represented by yellow circles) and the lack of disrupted tissue in the benign tissue (denoted by a lack of yellow circles).
- FIG. 5 shows a boxplot of the distribution of ⁇ for both benign and malignant
- FIG. 6 shows exemplary photographs of breast tissue classified according to provided methods. Specifically, the left panel shows control breast tissue made up of healthy fatty tissue (H ⁇ 0.45); the center panel shows a breast with a benign lesion lying at the interface between fatty and dense tissue, showing a few surrounding sub-regions of disrupted tissue (yellow, 0.45 ⁇ H ⁇ 0.55); and the right panel shows a breast with a malignant lesion with several disrupted tissue neighboring regions.
- H ⁇ 0.45 healthy fatty tissue
- the center panel shows a breast with a benign lesion lying at the interface between fatty and dense tissue, showing a few surrounding sub-regions of disrupted tissue (yellow, 0.45 ⁇ H ⁇ 0.55)
- the right panel shows a breast with a malignant lesion with several disrupted tissue neighboring regions.
- S (X)
- 7B shows the gradient of the image in 7A is obtained as the modulus of the wavelet transformed using equation 2.10 below;
- 7C shows maxima chains in blue defined as positions where the WT modulus is locally maximal (i.e. the WTMM) and along these chains in 7C further local maxima are found in red (i.e. the WTMMM); this is then repeated as several different scales, three of which are shown in 7D, 7E, and 7F; the WTMMM were then connected vertically through scales to define the WT skeleton shown in 7G, with the gray-scale coding going from black (minimum) to white (maximum).
- FIG. 8A shows an original image obtained from the DDSM database
- 8B shows a zoomed in image with a suspicious region encircled by a radiologist
- FIG. 9A-G shows the same analysis as in FIG8A-G, only on a different case.
- FIG. 10 shows the frequency distributions of fractal dimensions D calculated for benign cranio-caudal (CC) and mediolateral oblique (MLO) views (top panels) as compared to the fractal dimensions D calculated for cancer CC and MLO views (bottom panels); it is of note that the distributions are drastically different between benign and cancer samples.
- CC cranio-caudal
- MLO mediolateral oblique
- FIG. 11 shows a fractal dimension plot of several samples; specifically, each case is represented by a single dot and is plotted with the fractal dimension obtained from the mediolateral oblique (MLO) view as a function of the fractal dimension obtained from the cranio-caudal (CC) view.
- MLO mediolateral oblique
- CC cranio-caudal
- a polygonal region is outlined, the inside of which is defined as the "fractal zone” while the outside is defined as the "Euclidean zone”.
- the dots represent malignant (dark red) and benign (light green) cases.
- FIG. 12 shows exemplary random point distribution models to describe how different kinds of objects can have only a limited number of possible fractal dimension as a function of the projection angle; the left panel (12A) shows the possible fractal dimensions of a line; the center panel (12B) shows the possible fractal dimensions of a surface, and the right panel (12C) shows the possible fractal dimensions of a fractal cluster (here a simulated diffusion- limited aggregate).
- FIG. 13 shows a block diagram of an exemplary cloud computing environment suitable for use with provided methods.
- FIG. 14 is a block diagram of an exemplary computing device and a mobile computing device suitable for use with provided methods.
- the present invention provides, among other things, methods and systems for identifying cancerous or pre-cancerous tissue in a subject.
- the present invention is based, in part, on the surprising discovery that analysis of tissue surrounding a region of interest rather than tissue from the region of interest itself, may be used to determine if a particular region of tissue is cancerous or pre-cancerous.
- the present invention provides, among other things, methods of identifying cancerous or pre-cancerous tissue including providing a first region of tissue from a subject, calculating a roughness exponent for the first region of tissue, and comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 (e.g., less than 0.1, less than 0.05) between the roughness exponent of the first region of tissue and 0.5 indicates that the tissue is cancerous or precancerous.
- provided methods further comprise treating the tissue adjacent to the first region of tissue if it is cancerous or pre-cancerous.
- the present invention also provides methods including providing a first view of a region of tissue, providing a second view of a region of tissue, calculating a first fractal dimension for the first view of the region of tissue, and calculating a second fractal dimension for the second view of the region of tissue, wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
- the providing, calculating, and comparing steps may be advantageous to perform the providing, calculating, and/or comparing steps of provided methods multiple times.
- performing the providing, calculating and comparing steps multiple times will allow for superior characterization of the tissue or the tissue surrounding a tissue region of interest as compared to performing each step a single time according to various embodiments.
- the providing, calculating, and comparing steps are each performed a plurality of times.
- the plurality of times is at least 10, at least 20, at least 50, at least 100, at least 300, at least 500, or at least 1,000 times.
- performing the providing, calculating, and comparing steps multiple times occurs for different regions of the same organ or tissue.
- performing the providing, calculating, and comparing steps multiple times occurs for substantially the same region(s) of the same organ or tissue.
- performing the providing, calculating, and comparing steps multiple times occurs for different organs or tissues.
- MC microcalcifications
- mammograms are currently the most effective way of detecting breast cancer, it remains difficult to interpret mammograms due to high tissue variability and 3D to 2D projection effects.
- Current practice is to have two expert radiologists read the mammograms to reduce interpretation errors.
- CAD Computer Aided Diagnosis
- FDA Food and Drug Administration
- many of the current CAD methods rely on tissue homogeneity characterized by the Hurst exponent, meaning the tissue is assumed to be monofractal and uncorrected. Healthy fatty and dense tissue, shown by dark and light intensities respectively on mammograms, demonstrate how tissue may not be statistically homogeneous and there are many fluctuations of intensity throughout the image. By neglecting to accurately characterize the environment of possible breast lesions, important information is lost.
- the two scores provided by the 2D Wavelet-Transform Modulus Maxima (WTMM) method i.e. the fractal dimension of the MC and the roughness of the tissue surrounding the breast lesion given by the Hurst exponent, have known physical properties and can provide an insight to the invasiveness and severity of the possible lesion.
- WTMM Wavelet-Transform Modulus Maxima
- tissue composition of the breast By taking into account the tissue composition of the breast, radiologists can more accurately identify suspicious regions. Since tissue composition should be symmetric between both breasts, differing densities maybe a potential sight for a tumor, causing
- the databank contains over 2,500 studies made up of normal, benign and malignant mammograms all categorized by an expert radiologist. Each study has two images of each breast, consisting of a mediolateral oblique (MLO) view and cranio- caudal (CC) view with any suspicious region circled by a radiologist.
- MLO mediolateral oblique
- CC cranio- caudal
- the suspicious region could contain a mass and/or micro calcifications (MC), but only the cases that were classified as normal, i.e. no suspicious area was identified, were looked at in the Kestener et al. study.
- Various embodiments of the present invention apply wavelet-transform modulus maxima (WTMM) to one or more images of a tissue in order to characterize one or more attributes of the tissue and/or tissue adjacent to the tissue (i.e., the tissue's microenvironment).
- WTMM wavelet-transform modulus maxima
- the present invention evidences the potential of the two dimensional WTMM method to become a powerful tool in interpreting mammograms.
- the WTMM method has proven to be successful in several fields of applied science, including geology, astrophysics, cellular biology and orthopedic medicine.
- the method was originally developed as a multifractal formalism to analyze highly complex ID signals, 2D images, and 3D images.
- the wavelet transform acts as a mathematical microscope to characterize spatial image information over a continuous range of size scales. It is the gradient vector of the image smoothed by dilated versions of a Gaussian filter.
- WTMM wavelet transform modulus maxima
- WTMMM WTMM maxima
- provided methods make use of an adaptation of a two dimensional WTMM method, specifically, the continuous wavelet transform method, as a mathematical microscope used to characterize the fractal geometry of clusters of
- microcalcifications in human tissue (e.g., breast tissue) and/or to determine the roughness of the background tissue (i.e. tissue adjacent to the tissue of interest) seen in mammograms.
- the WTMM method yields the so-called singularity spectrum, D(h), i.e. the fractal dimension D, of points having a Holder exponent of h.
- the WTMM method is used to perform a segmentation of tissues, in some embodiments, the breast tumor tissue, based on the strength of the singularities composing the mammogram images, and to simultaneously quantify their fractal dimension.
- the fractal dimensions of benign and malignant breast lesions are significantly different, with benign having an integer dimension corresponding to a non-invasive Euclidean object and cancer having a non- integer dimension, representing an invasive structure.
- the microenvironments characterized by the roughness of the tissue in which the lesions are embedded are different for benign and malignant tumors, and provides an insight into the onset and development of breast cancer.
- One score provided by the two dimensional WTMM method is the roughness of the tissue surrounding the breast lesion, given by the Hurst exponent, which has known physical properties and can provide an insight to the invasiveness and severity of the possible lesion.
- WTMM multi-scale methods
- wavelets leaders detrended fluctuation analysis, Fourier transform methodologies, and/or any mathematically related and/or similar multi-scale density fluctuation assessment methods yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent.
- the Hurst exponent represents, in part, the physical status of the microenvironment of the tissue (e.g., tissue density as determined from roughness fluctuations) and the level of correlation in that physical status across adjacent tissue.
- H ⁇ 0.5 the roughness fluctuation are considered anti-correlated
- H > 0.5 means that the roughness fluctuations are positively correlated.
- the physical system of the analyzed tissue is considered to have some form of spatial memory.
- the present invention recognizes that the results obtained by using the WTMM method on screening mammograms could not only be used as a possible computer-aided detection (CAD) method, but also as a method to further study the biophysics of tumor onset and progression. As described in the Examples below, after all images were analyzed using provided methods, a statistical analysis was performed on all data, which provided us with information on the critical differences between the organization, behavior, and biological and physical processes of both benign and malignant tumors.
- CAD computer-aided detection
- This methodology has the ability to be applied to a variety of rough surfaces, such as mammogram landscapes.
- rough surface for an irregular surface. This means the surface can be accurately described by a single -valued, self-affme function satisfying the following:
- the Hurst exponent is used to characterize the global roughness of the function under investigation. If H ⁇ 1 , then the function, f, is nowhere
- edges of the different structures that appear in an image are often the most important features for pattern recognition. Therefore, many computer tools for edge detection look for points where the gradient of the image intensity has a modulus locally maximum in its direction. With an appropriate choice of analyzing wavelets, one can redefine the Canny's multi- scale edge detection in terms of a 2D wavelet transform.
- a multi-scale edge detection methodology is desirable to assist in determining the borders of cancerous or pre-cancerous tissue.
- exemplary multi-edge detection techniques that may be used in accordance with the present invention.
- equation 2.7 amounts to define the 2D wavelet transform as the gradient vector of f(x) smoothed by dilated versions of ⁇ (a _1 x) of this filter.
- the wavelet transform modulus maxima are defined as the points b where the wavelet transform modulus, M [f](b, a), is locally
- WTMM WTMM maxima
- the 2D WTMM method relies upon the space-scale partitioning given by the WT skeleton. Let us define L(a) as the set of maxima lines that exist at the scale a and which contain maxima at any scale a ' ⁇ a.
- the WTMM method consists in defining the following partition functions directly from the WTMMM that belong to the wavelet transform skeleton:
- q ER Compared to the classical box-counting techniques, the analyzing wavelet, ⁇ , plays the role of a generalized "oscillating box", the scale, a, defines its size, while the WT skeleton indicates how to position our boxes to obtain a partition at the considered scale. From the analogy that links the multifractal formalism to thermodynamics, one can define the exponent q) from the power law behavior of the partition function: [0071]
- q and r( ) are the inverse temperature and the free energy respectively.
- WTMMM sup,-., .,, -r ⁇ ,' ⁇ -, - ⁇ ⁇ f ifx. a' ⁇ ⁇ q
- the roughness analysis of tissue is not performed with two dimensional images.
- the WTMM method as well as all other similar space-scale techniques that allow the calculation of the Hurst exponent (or any analogous exponent), such as, but not limited to wavelet leaders, Fourier techniques, etc are also directly applicable to 3D data cubes.
- the assessment of cancerous or pre-cancerous tissue regions through a roughness analysis is applicable to 3D data.
- 3D data the distance between the estimated Hurst exponent (or other analogous exponent) and the critical value 0.5 is the same as for the 2D case.
- the strategy of analyzing square sub-regions neighboring a tumor or suspicious tissue area such as presented in FIG. 6 is directly generalizable to the analysis of cubic sub-regions in 3D.
- the present invention provides methods and systems for detecting microcalcifications (MC) in subject tissue and characterizing their fractal dimension.
- provided methods use two dimensional (or, in some embodiments, three dimensional) WTMM to provide this information, though any mathematically related and/or similar multi-scale density fluctuation assessment method yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent and/or fractal dimension may be used.
- MC MC
- MCs typically appear as bright spots on mammograms and the number, size, type (Table 2) and distribution (Table 3) provide the radiologist with information about the potential malignancy of the lesion and allow an assessment score (Table 4) to be assigned to each case.
- radiologists can predict the potential malignancy of the identified tumor on the mammogram. For example, tumors that contain pleomorphic or fine linear branching calcification types have a higher risk of being malignant than tumors that are classified as punctate. However, it is not uncommon for tumors to be made up of components that make it difficult for radiologists to the identify the tumor as cancerous or a benign disease such as fibrocystic breast condition or a blunt injury.
- One benefit of a CAD is the objectivity of the software of having a range of indices associated with benign and/or malignant tumors. Identifying Microcalcifications (MC)
- the present invention provides computer aided diagnostic tools, using the WT methodology to detect MC by inspecting the WT maxima chains.
- WT the WT methodology
- 7 pixels, MC which can be considered at strong singularities are contour-shaped by some maxima chains. Since the average size of MC is about 200 ⁇ , or 5 pixels, these singularities are identified by the mathematical microscope as Dirac singularities.
- the MC wavelet-transform skeleton can be used to compute the partition functions, and thus fully characterize the fractal geometry of the MC cluster.
- the two dimensional WTMM method has the ability to be a powerful CAD.
- the information obtained from the wavelet-transform analysis can help quantify data that would otherwise be considered subjective.
- provided methods use this information to aid in the development of an objective mammographic classification, which can be used to judge the potential malignancy of a radiologist-identified suspicious region.
- characterization of the fractal geometry of MC clusters from a plurality of separate two-dimensional views of the same region of tissue allows for an assessment of the overall fractal dimension of the tissue region of interest, such as through application of the WT methods described herein.
- a plurality of separate two-dimensional views of the same region of tissue one can assess the fractal dimension of the three dimensional tissue region of interest.
- non-cancerous or non-pre-cancerous tissue will have a
- Euclidean fractal dimension i.e. an integer fractal dimension
- a cancerous or pre-cancerous tissue will have a non-integer fractal dimension indicating an invasive morphology.
- the present invention also provides methods including providing a first view of a region of tissue, providing a second view of a region of tissue, calculating a first fractal dimension for the first view of the region of tissue, and calculating a second fractal dimension for the second view of the region of tissue, wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
- provided methods further comprise treating the region of tissue if it is cancerous or pre-cancerous.
- the "fractal zone” is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
- a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8); first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8); second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2); third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
- a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9); first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9); second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1); third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
- a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7); first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7); second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3); third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
- the precise boundaries of the fractal zone may vary according to tissue type and/or tumor type (or subtype). Additionally, without wishing to be held to a particular theory it is possible that the precise borders of the fractal zone may change over time as tissue changes from normal tissue to pre-cancerous or cancerous tissue.
- tissue type and/or tumor type or subtype
- the precise borders of the fractal zone may change over time as tissue changes from normal tissue to pre-cancerous or cancerous tissue.
- One of skill in the art will be able to ascertain from the present disclosure how to determine the appropriate boundaries of the fractal zone in a particular clinical scenario using no more than routine experimentation.
- provided methods may be applied to more than two views of a tissue region of interest. In some embodiments, provided methods may be applied to three, four, five, six, seven, eight, nine, ten, fifteen, twenty or more views of a particular tissue region of interest.
- One of skill will understand how the addition of additional two-dimension images to the analysis of provided methods will alter the definition of the fractal zone. For example, if three views of a tissue region of interest are used, the fractal zone may then be defined as the union of seven juxtaposed regions: one central cube and six extending pyramidal or conical regions.
- the borders of the fractal zone may comprise other forms and definitions.
- the borders of the fractal zone may be smooth curves, jagged lines, or even fractals.
- any fractal technique or techniques yielding a fractal dimension such as, but not limited to: box-counting techniques, the perimeter-area relationship, packing dimension, Hausdorff dimension, capacity dimension, correlation dimension, the generalized dimensions (multifractal), and/or any space-scale technique that yields a value quantifying the structural complexity of a tumor or object of interest as a function of scale may be used.
- tissue that may contain or comprise cancer.
- Non-limiting exemplary such tissues include breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue. Treatment of Cancerous or Pre-Cancerous Tissue
- provided methods may be used to detect and/or characterize the presence of cancerous or pre-cancerous tissue in a subject. It is contemplated that if cancerous or pre-cancerous tissue is discovered, that treatment of such tissue would commence in accordance with sound medical judgment. It is further contemplated that treatment may occur with any anti-cancer therapy, whether currently known or discovered in the future.
- anti-cancer therapy may be appropriate for use in some embodiments, exemplary types of anti-cancer therapy are described below in order to better illustrate some of the principles of the present invention.
- brachytherapy and systemic radioactive isotopes, to name a few
- endocrine therapy biologic response modifiers (interferons, interleukins, and tumor necrosis factor (TNF) to name a few)
- TNF tumor necrosis factor
- hyperthermia and cryotherapy agents to attenuate any adverse effects (e.g., antiemetics), and other approved chemotherapeutic drugs, including, but not limited to, alkylating drugs
- the present invention also provides systems and devices, such as computing devices, implementing provided methods.
- a computing device maybe a cloud computing device, a mobile computing device, or any other application-appropriate computing device.
- the cloud computing environment 1300 may include one or more resource providers 1302a, 1302b, 1302c (collectively, 1302).
- Each resource provider 1302 may include computing resources.
- computing resources may include any hardware and/or software used to process data.
- computing resources may include hardware and/or software capable of executing algorithms, computer programs, and/or computer applications.
- exemplary computing resources may include application servers and/or databases with storage and retrieval capabilities.
- Each resource provider 1302 may be connected to any other resource provider 1302 in the cloud computing environment 1300.
- the resource providers 1302 may be connected over a computer network 1308.
- Each resource provider 1302 may be connected to one or more computing device 1304a, 1304b, 1304c (collectively, 1304), over the computer network 1308.
- the cloud computing environment 1300 may include a resource manager 1306.
- the resource manager 1306 may be connected to the resource providers 1302 and the computing devices 1304 over the computer network 1308. In some implementations, the resource manager
- 1306 may facilitate the provision of computing resources by one or more resource providers
- the resource manager 1306 may receive a request for a computing resource from a particular computing device 1304.
- the resource manager 1306 may identify one or more resource providers 1302 capable of providing the computing resource requested by the computing device 1304.
- the resource manager 1306 may select a resource provider 1302 to provide the computing resource.
- the resource manager 1306 may facilitate a connection between the resource provider 1302 and a particular computing device 1304.
- the resource manager 1306 may establish a connection between a particular resource provider 1302 and a particular computing device 1304.
- the resource manager 1306 may redirect a particular computing device 1304 to a particular resource provider 1302 with the requested computing resource.
- FIG. 14 shows an example of a computing device 1400 and a mobile computing device 1450 that can be used to implement the techniques described in this disclosure.
- the computing device 1400 is intended to represent various forms of digital computers, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers.
- the mobile computing device 1450 is intended to represent various forms of mobile devices, such as personal digital assistants, cellular telephones, smart- phones, and other similar computing devices.
- the components shown here, their connections and relationships, and their functions, are meant to be examples only, and are not meant to be limiting.
- the computing device 1400 includes a processor 1402, a memory 1404, a storage device 1406, a high-speed interface 1408 connecting to the memory 1404 and multiple highspeed expansion ports 1410, and a low-speed interface 1412 connecting to a low-speed expansion port 1414 and the storage device 1406.
- Each of the processor 1402, the memory 1404, the storage device 1406, the high-speed interface 1408, the high-speed expansion ports 1410, and the low-speed interface 1412 are interconnected using various busses, and may be mounted on a common motherboard or in other manners as appropriate.
- the processor 1402 can process instructions for execution within the computing device 1400, including instructions stored in the memory 1404 or on the storage device 1406 to display graphical information for a GUI on an external input/output device, such as a display 1416 coupled to the high-speed interface 1408.
- an external input/output device such as a display 1416 coupled to the high-speed interface 1408.
- multiple processors and/or multiple buses may be used, as appropriate, along with multiple memories and types of memory.
- multiple computing devices may be connected, with each device providing portions of the necessary operations (e.g., as a server bank, a group of blade servers, or a multi-processor system).
- the memory 1404 stores information within the computing device 1400.
- the memory 1404 is a volatile memory unit or units.
- the memory 1404 is a non-volatile memory unit or units.
- the memory 1404 may also be another form of computer-readable medium, such as a magnetic or optical disk.
- the storage device 1406 is capable of providing mass storage for the computing device 1400.
- the storage device 1406 may be or contain a computer- readable medium, such as a floppy disk device, a hard disk device, an optical disk device, or a tape device, a flash memory or other similar solid state memory device, or an array of devices, including devices in a storage area network or other configurations.
- Instructions can be stored in an information carrier.
- the instructions when executed by one or more processing devices (for example, processor 1402), perform one or more methods, such as those described above.
- the instructions can also be stored by one or more storage devices such as computer- or machine- readable mediums (for example, the memory 1404, the storage device 1406, or memory on the processor 1402).
- the high-speed interface 1408 manages bandwidth-intensive operations for the computing device 1400, while the low-speed interface 1412 manages lower bandwidth-intensive operations. Such allocation of functions is an example only.
- the highspeed interface 1408 is coupled to the memory 1404, the display 1416 (e.g., through a graphics processor or accelerator), and to the high-speed expansion ports 1410, which may accept various expansion cards (not shown).
- the low-speed interface 1412 is coupled to the storage device 1406 and the low-speed expansion port 1414.
- the low-speed expansion port 1414 which may include various communication ports (e.g., USB, Bluetooth®, Ethernet, wireless Ethernet) may be coupled to one or more input/output devices, such as a keyboard, a pointing device, a scanner, or a networking device such as a switch or router, e.g., through a network adapter.
- input/output devices such as a keyboard, a pointing device, a scanner, or a networking device such as a switch or router, e.g., through a network adapter.
- the computing device 1400 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a standard server 1420, or multiple times in a group of such servers. In addition, it may be implemented in a personal computer such as a laptop computer 1422. It may also be implemented as part of a rack server system 1424. Alternatively, components from the computing device 1400 may be combined with other components in a mobile device (not shown), such as a mobile computing device 1450. Each of such devices may contain one or more of the computing device 1400 and the mobile computing device 1450, and an entire system may be made up of multiple computing devices communicating with each other.
- the mobile computing device 1450 includes a processor 1452, a memory 1464, an input/output device such as a display 1454, a communication interface 1466, and a transceiver 1468, among other components.
- the mobile computing device 1450 may also be provided with a storage device, such as a micro-drive or other device, to provide additional storage.
- a storage device such as a micro-drive or other device, to provide additional storage.
- Each of the processor 1452, the memory 1464, the display 1454, the communication interface 1466, and the transceiver 1468, are interconnected using various buses, and several of the components may be mounted on a common motherboard or in other manners as appropriate.
- the processor 1452 can execute instructions within the mobile computing device
- the processor 1452 may be any circuitry 1450, including instructions stored in the memory 1464.
- the processor 1452 may be any circuitry 1452.
- the processor 1452 may provide, for example, for coordination of the other components of the mobile computing device 1450, such as control of user interfaces,
- applications run by the mobile computing device 1450, and wireless communication by the mobile computing device 1450.
- the processor 1452 may communicate with a user through a control interface
- the display 1454 may be, for example, a TFT (Thin-Film-Transistor Liquid Crystal Display) display or an OLED (Organic Light Emitting Diode) display, or other appropriate display technology.
- the display interface 1456 may comprise appropriate circuitry for driving the display 1454 to present graphical and other information to a user.
- the control interface 1458 may receive commands from a user and convert them for submission to the processor 1452.
- an external interface 1462 may provide communication with the processor 1452, so as to enable near area communication of the mobile computing device 1450 with other devices.
- the external interface 1462 may provide, for example, for wired communication in some implementations, or for wireless communication in other implementations, and multiple interfaces may also be used.
- the memory 1464 stores information within the mobile computing device 1450.
- the memory 1464 can be implemented as one or more of a computer-readable medium or media, a volatile memory unit or units, or a non-volatile memory unit or units.
- expansion memory 1474 may also be provided and connected to the mobile computing device 1450 through an expansion interface 1472, which may include, for example, a SIMM (Single In Line Memory Module) card interface.
- the expansion memory 1474 may provide extra storage space for the mobile computing device 1450, or may also store applications or other information for the mobile computing device 1450.
- the expansion memory 1474 may include instructions to carry out or supplement the processes described above, and may include secure information also.
- the expansion memory 1474 may be provided as a security module for the mobile computing device 1450, and may be programmed with instructions that permit secure use of the mobile computing device 1450.
- secure applications may be provided via the SIMM cards, along with additional information, such as placing identifying information on the SIMM card in a non-hackable manner.
- the memory may include, for example, flash memory and/or NVRAM memory
- instructions are stored in an information carrier and, when executed by one or more processing devices (for example, processor 1452), perform one or more methods, such as those described above.
- the instructions can also be stored by one or more storage devices, such as one or more computer- or machine -readable mediums (for example, the memory 1464, the expansion memory 1474, or memory on the processor 1452).
- the instructions can be received in a propagated signal, for example, over the transceiver 1468 or the external interface 1462.
- the mobile computing device 1450 may communicate wirelessly through the communication interface 1466, which may include digital signal processing circuitry where necessary.
- the communication interface 1466 may provide for communications under various modes or protocols, such as GSM voice calls (Global System for Mobile communications), SMS (Short Message Service), EMS (Enhanced Messaging Service), or MMS messaging (Multimedia Messaging Service), CDMA (code division multiple access), TDMA (time division multiple access), PDC (Personal Digital Cellular), WCDMA (Wideband Code Division Multiple Access), CDMA2000, or GPRS (General Packet Radio Service), among others.
- GSM voice calls Global System for Mobile communications
- SMS Short Message Service
- EMS Enhanced Messaging Service
- MMS messaging Multimedia Messaging Service
- CDMA code division multiple access
- TDMA time division multiple access
- PDC Personal Digital Cellular
- WCDMA Wideband Code Division Multiple Access
- CDMA2000 Code Division Multiple Access
- GPRS General Packet Radio Service
- GPS Global Positioning System
- short-range communication may occur, such as using a Bluetooth®, Wi-FiTM, or other such transceiver (not shown).
- a GPS (Global Positioning System) receiver module 1470 may provide additional navigation- and location-related wireless data to the mobile computing device 1450, which may be used as appropriate by applications running on the mobile computing device 1450.
- the mobile computing device 1450 may also communicate audibly using an audio codec 1460, which may receive spoken information from a user and convert it to usable digital information.
- the audio codec 1460 may likewise generate audible sound for a user, such as through a speaker, e.g., in a handset of the mobile computing device 1450.
- Such sound may include sound from voice telephone calls, may include recorded sound (e.g., voice messages, music files, etc.) and may also include sound generated by applications operating on the mobile computing device 1450.
- the mobile computing device 1450 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a cellular telephone 1480. It may also be implemented as part of a smart-phone 1482, personal digital assistant, or other similar mobile device.
- Various implementations of the systems and techniques described here can be realized in digital electronic circuitry, integrated circuitry, specially designed ASICs (application specific integrated circuits), computer hardware, firmware, software, and/or combinations thereof.
- ASICs application specific integrated circuits
- These various implementations can include implementation in one or more computer programs that are executable and/or interpretable on a programmable system including at least one programmable processor, which may be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.
- machine-readable medium and computer-readable medium refer to any computer program product, apparatus and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a machine- readable medium that receives machine instructions as a machine -readable signal.
- machine-readable signal refers to any signal used to provide machine instructions and/or data to a programmable processor.
- the systems and techniques described here can be implemented on a computer having a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the computer.
- a display device e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor
- a keyboard and a pointing device e.g., a mouse or a trackball
- Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form, including acoustic, speech, or tactile input.
- the systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components.
- the components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), and the Internet.
- LAN local area network
- WAN wide area network
- the Internet the global information network
- the computing system can include clients and servers.
- a client and server are generally remote from each other and typically interact through a communication network.
- the relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
- This Example shows that the 2D WTMM method has the ability to characterize breast tumors and their microenvironment.
- Data was obtained from the Digital Database for Screening Mammography (DDSM), described above. Cases having a suspicious region, i.e. containing a benign or malignant tumor, were considered. From these cases, only cases having one set of microcalcifications were kept for further investigation and masses were disregarded for this analysis. Using the methods described above, we were able to fully characterize the lesion by the fractal dimension of the breast tumor and the roughness of the microenvironment given by the Hurst exponent.
- DDSM Digital Database for Screening Mammography
- Benign breast lesions' three components are centered at (.84,1.04), (.85,1.78) and (1.78,1.88), while malignant breast lesions' two components are centered at (1.07,1.62) and (1.49,1.52).
- the center of the benign categories demonstrate how these tumors are noninvasive since the populations fractal dimensions lay in Euclidean space and are represented by lines, sheets and spheres. However, both malignant categories are shown by invasive fractal dimensions and represented by fractal structures in breast tissue.
- This example illustrates how provided methods may be used to better understand the mechanisms that drive the differing organization of benign and malignant tumors, and specifically the microenvironment of a radiologist-identified suspicious region.
- each image of a breast lesion was segmented into 9 subimages, corresponding to the same size and shape of the radiologist encircled tumor. Once the images were properly segmented, only the central part of the neighboring images were analyzed to counter any edge effects that may disrupt the analysis. The rough surface images were analyzed via the above methodology.
- FIG. 4 the results allowed for the characterization of abnormal tissue corresponding to a benign tumor and malignant tumor.
- red (0.55 - 1), or dense tissue
- yellow (0.45; 0.55).
- H (0.45 - 0.55
- the boxplot in FIG. 5 shows the distribution of ⁇ for both benign and malignant microenvironments.
- the max and min values are shown on the plot as 0.33 and 0.07 for benign and 0.26 and 0.08 for cancer respectively.
- three outliers were identified for benign and one for cancer, meaning these values for ⁇ are 1.5 times more than the maximum ⁇ value.
- m (mi; im; .. . ; ⁇ ) ⁇ and mi is expected value of the order statistic and V is the covariance matrix.
- Example 2 As indicated above in Example 2, obtaining one average ⁇ value per subject by averaging the calculated
- Example 4 Characterization of Tissue Microenvironment
- data from analyzed breast tissue was characterized on a per neighboring region basis, as opposed to a per patient basis.
- a total of 1,131 images of tissue classified as malignant, 1,131 images of tissue classified as benign, and 64 images of tissue classified as normal were analyzed.
- the number of regions in analyzed microenvironments characterized as fatty, dense or disrupted were calculated from the samples examined in Example 1.
- Table 8 shows the results of this analysis:
- the percentage of disrupted regions was highest in malignant tissue, followed by benign tissue, with normal tissue having the lowest levels of disruption across the tissue. Exemplary photographs of the differences in observed disruption may be found in FIG. 6. It is important to note that the percentage of disrupted regions for normal breast tissue is not 0%>. Without wishing to be held to a particular theory, it is possible these disrupted regions may be representative of a transitional state of the tissue that may be prone to the eventual development of a tumor. Longitudinal studies will help to determine the nature of disrupted regions in breast tissue classified as normal by current diagnostic
- DDSM Database for Screening Mammography
- the databank contains over 2,600 studies made up of normal, benign, benign without call back and malignant mammograms all categorized by an expert radiologist. Each study has two images of each breast, consisting of a mediolateral-oblique (MLO) view and cranio-caudal (CC) view with any suspicious region circled by a radiologist.
- MLO mediolateral-oblique
- CC cranio-caudal
- the suspicious region could contain a mass and/or microcalcifications (MC), but only the cases that were classified as having exactly one tumor composed of only MC in the benign and malignant categories were looked at in this particular study.
- any mammographic images that contained artifacts inside the radiologist's encircled region were discarded due to the impact it has on the analysis.
- artifacts could include scratches, hair, deodorant, patient movement, scanner artifacts (rollers slipped), pacemaker, breast implants, skin markers (for scars, moles, and nipples, as well as marked lumps of breast pain), metallic foreign bodies, and fingerprints.
- Some (but not all) of these effects were recorded under notes in the DDSM website.
- a total of 59 cases were considered, corresponding to 118 images of size greater than 2562 pixels, 34 of which are benign (68 images) and 25 of which are malignant (50 images).
- the 2D WTMM method requires us to define a smoothing function, 0(x, y), in two dimensions that is a well-localized isotropic function around the origin.
- a smoothing function 0(x, y)
- the Gaussian function we used the Gaussian function and define the wavelets as:
- the wavelet transform modulus maxima are defined as the locations b where M ⁇ , ⁇
- the WTMM maxima, or WTMMM are defined as the points along the maxima chains where the ⁇ ⁇ [/] (b,a) is locally maximum.
- WTMMM are linked through scales to form the space-scale skeleton (Fig. 7G).
- Fig. 7G space-scale skeleton
- h and D(h) as mean quantities defined in a canonical ensemble, i.e. with respect to their Boltzmann weights computed from the WTMMM:
- H are characterized by a linear ⁇ (q) curve of slope H.
- D F (noted simply D throughout the text) is the fractal dimension of the support of singularities of f.
- the Wilcoxon rank-sum test is a non-parametric statistical hypothesis test that is used as an alternative to Student's t-test when the population cannot be assumed to be normally distributed. It was used here to calculate the / ⁇ -values comparing the CC and MLO fractal dimensions and benign and malignant cases images since the benign data followed a bimodal distribution (FIG. 10). The calculations were done using the Wilcox test in R.
- ⁇ ( ⁇ ⁇ 0) p(G ⁇ 9) * ⁇ ( ⁇ ) / p(G)
- ⁇ ( ⁇ ) represents the strength of our belief in malignant lesions p(M) or benign lesions p(B) out of those that have been diagnosed by a radiologist.
- the posterior, p(9IG) represents the strength of our belief, having accounted for the geometrical evidence, G, where G represents the position of the lesion in the fractal dimension plot, either fractal (F) or Euclidean (E).
- the quotient of the likelihood over the evidence, p(9IG)/p(G) represents the support the evidence, G, provides for ⁇ . Since the prior reflects uncertainty in the parameter value ⁇ , ⁇ ( ⁇ ) was based on a Beta distribution with a specified mean and standard deviation. The Beta distribution is as follows:
- ⁇ ,- ⁇ , ⁇ ) ⁇ ( ⁇ + ⁇ ) ⁇ ⁇ - ⁇ - ⁇ ⁇
- the Probability Model To estimate the mean of the Beta distribution for malignant cases, p(M), the prevalence of mammograms with a BI-RADS assessment score of 3, 4, and 5 were used as determined by the radiologists diagnostics multiplied by the historical probability of mammograms receiving those assessment scores result in malignant MC clusters respectively.
- the probability of malignant lesions given an assessment score of 3 is 2%
- an assessment score of 4 is 26.5% (taken as the midpoint of the reported interval of [23%>-30%>]]
- an assessment score of 5 is 95%. Therefore the Beta distribution for p(M) was chosen with a mean of
- Beta distribution for p(B) with a mean of 1-0.3885 0.6115.
- the likelihood p(FlM), where F represents breast lesions characterized as being in the fractal 465 zone (and E likewise represents those in the Euclidean zone), is based on the 23 of the 25 malignant cases that were in the fractal zone; this likelihood is the probability that the data could be generated with parameter values ⁇ .
- the likelihood p(ElB) is based on the 30 out of 34 benign cases that were in the Euclidean zone.
- the wavelet transform acts as a mathematical microscope to characterize spatial image information over a continuous range of size scales. It is the gradient vector of the image smoothed by dilated versions of a Gaussian filter.
- WTMM wavelet transform modulus maxima
- each space-scale line obtained from the WT skeleton is represented by plotting the evolution of the WT modulus, M (see the 2D WTMM section above), as a function of the scale parameter, a, in a log-log plot. This relationship between M and a is characterized by a power-law behavior via the equation
- the background breast tissue is dense, which makes the contrast between background and MC weak (i.e. causing a low value for the WT modulus of red curves at the smallest scale in FIG. 8D).
- the roughness fluctuations of dense breast tissue are characterized by a relatively high smoothness level, which translates to blue curves with a large slope (i.e., a high h value, -2/3) for scales 10 ⁇ a ⁇ 200 pixels as compared to the red curves with negative slopes for scales a > 10 pixels that correspond to WTMMM lines that point to MC at small scale (a ⁇ 0+) (FIG. 8D).
- a high h value a large slope
- the background breast tissue is fatty, which is characterized by a higher roughness level (i.e. a lower h value -1/3, although still positive), that reduces the discriminatory power of h.
- the contrast is high, which translates to a high value of log( T). Therefore, applying a threshold on both parameters, h and log(K), is key to segregating MC from their background regardless of the density (fatty or dense) of the composition of the breast tissue.
- Euclidean zones are avoided and the data are very clearly centered in the fractal range for both views, with the peak of the
- Example 5 a fractal dimension plot between the MLO and CC views of a breast is calculated. Unless otherwise specified, all samples, sample handling, and analysis was as described above in Example 5. As shown in Example 5, the significance of the difference between benign and malignant is quite interesting. However, it is still only based on statistics of populations. The histograms in FIG. 10 show that, when each view is taken independently (CC or MLO), it is still possible, though unlikely, for a malignant lesion to have a Euclidean dimension, and vice-versa, for benign lesions to have a (non-integer) fractal dimension.
- Bayesian models are designed to be appropriate to the data structure without having to make approximation assumptions typical in null hypothesis significance testing.
- the results reported above show that the vast majority of malignant breast lesions are fractal, and that the vast majority of benign breast lesions are Euclidean.
- the condition of interest is how breast lesions in the fractal zone can indicate malignancy, and how breast lesions in the Euclidean zone can indicate benignancy.
- Bayesian inference derives the posterior probability as a consequence of two antecedents, a prior probability and a likelihood function derived from a probability model for the data to be observed.
- the model is based on historical radiology assessment scores using the BI-RADS system.
- Bayesian inference then becomes a computation of the posterior probability according to Bayes' rule.
- the interpretable output of this Bayesian analysis is the so-called 95% highest density interval (HDI), which is analogous to the 95% confidence interval in frequentist statistics.
- HDI 95% highest density interval
- the 95% HDI from the resulting posterior distribution indicates that the percentage of breast lesions in the fractal zone that are malignant is between 74.2% and 97.5%.
- the percentage of breast lesions in the Euclidean zone that are benign is between 75.7% and 96.2%.
- LINES i.e. those that are approximately in the area, which are therefore seen as one-dimensional objects from both views (FIG. 12A, left panel)
- [1.2,1.8] x[1.2,1.8] square have a fractal signature that is seen from both views, whereas those that are in the triangular areas would represent fractal clusters that grow onto a 2D plane, i.e. seen as a fractal from one view, but seen either as a line (bottom or left triangular regions) or a plane (top or right triangular regions) from the other view (FIG. 12C, right panel).
- a diffusion-limited aggregate embedded in 3D space and for which 2 ⁇ D ⁇ 3 will have a 2D projection with D 2.
- Examples 5 and 6 illustrate embodiments of provided methodologies which offer a way to accurately classify benign and malignant tumors based on their invasiveness as determined by the geometrical structure.
- current mammographic practice may be improved by increasing accuracy, and potentially decreasing recall rates and costs.
- the inferred 3-dimensional geometry of the breast lesions based on the analysis of the mammographic images using the 2D WTMM methodology allows one to explore the invasiveness of the breast tumors and provide an interpretation of the severity of the lesion.
- a score similar to the Breast Imaging-Reporting and Data System (BI-RADS) assessment score may be assigned to each case.
- BIOS Breast Imaging-Reporting and Data System
- this tool may also provide insight into the underlying mechanisms that drive overall growth and structure at the time of the screening mammogram.
- benign lesions likely being smooth Euclidean objects and malignant lesions being branching objects (and possibly, for both cases, being restricted to growing along 2D manifolds within the 3D breast tissue environment)
- branching objects and possibly, for both cases, being restricted to growing along 2D manifolds within the 3D breast tissue environment
- Use of provided methods may be helpful in revealing these relationships and improving detection and treatment.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Veterinary Medicine (AREA)
- Animal Behavior & Ethology (AREA)
- Public Health (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Surgery (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Artificial Intelligence (AREA)
- Psychiatry (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Probability & Statistics with Applications (AREA)
- Mathematical Physics (AREA)
- Quality & Reliability (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Optics & Photonics (AREA)
- Gynecology & Obstetrics (AREA)
- Reproductive Health (AREA)
- Apparatus For Radiation Diagnosis (AREA)
- Investigating Or Analysing Biological Materials (AREA)
- Neurology (AREA)
- Pharmaceuticals Containing Other Organic And Inorganic Compounds (AREA)
- Medicines That Contain Protein Lipid Enzymes And Other Medicines (AREA)
- Psychology (AREA)
- Neurosurgery (AREA)
Abstract
The present invention provides, among other things, methods of identifying cancerous or pre-cancerous tissue including providing a first region of tissue from a subject, calculating a roughness exponent for the first region of tissue, and comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 between the roughness exponent of the first region of tissue and 0.5 indicates that the tissue is cancerous or pre-cancerous.
Description
IMPROVED METHODS OF CANCER DETECTION
BACKGROUND
[0001] The human body is composed of about ten trillion cells, and throughout a lifetime cells go through many damaging processes that concern their genetic programming. However, not all of these mutations result in cancer. As it turns out, many individuals actually have occult tumors throughout their bodies that are only discovered through microscopic investigation during autopsy. Therefore, there must be some mechanism that prohibits the growth and development of tumors.
[0002] Recent research has pointed to the microenvironment of potential tumors to help suppress malignant phenotype and instruct otherwise malignant cells to participate in normal development. The microenvironment may actually provide tumor suppressive signals as long as the tissue architecture is controlled, but as the structure of healthy tissue is lost, the tissue has the potential to become a tumor promoter. Thus, considering the tissue surrounding the tumor, along with the lesion, may be important in identifying potential cancerous tissue.
SUMMARY OF THE INVENTION
[0003] The present invention is based, in part, on the realization that disorder of the microenvironment surrounding a portion of tissue in a subject may be indicative of the presence of, or pre-disposition toward developing, cancer in that tissue. Among other things, the present invention provides, for the first time, methods for analyzing the amount and distribution of disorder across tissue in order to predict the presence of, or propensity of tissue to develop, cancer. In some embodiments, provided methods include analysis of images of target tissue, as opposed to analysis of the tissue itself. In some embodiments, provided methods include analysis of tissue or images of tissue adjacent to a target tissue/tissue region of interest, as opposed to the target tissue itself.
[0004] The present invention provides, among other things, methods of identifying cancerous or pre-cancerous tissue including providing a first region of tissue from a subject, calculating a roughness exponent (e.g., a Hurst exponent) for the first region of tissue, and
comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 between the roughness exponent of the first region of tissue and 0.5 indicates that the tissue is cancerous or pre-cancerous. In some embodiments, the providing, calculating, and comparing steps are each performed a plurality of times. In some embodiments, the plurality of times is at least 10, at least 100, or at least 1,000. In some embodiments, performing the providing, calculating and comparing steps a plurality of times allows for a higher resolution characterization of the morphology of the tissue surrounding a region of tissue suspected to be cancerous or pre-cancerous.
[0005] In some embodiments, the difference between the roughness exponent values and
0.5 may be less than 0.2. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.15. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.1. In some embodiments, the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.05.
[0006] In some embodiments, a roughness exponent is calculated from one or more images of a target tissue. In some embodiments, a roughness exponent is calculated from two or more images of a first region of tissue. In some embodiments, each of the two or more images are taken from different angles and/or points of view. In some embodiments, a roughness exponent is calculated for two or more regions of tissue. In some embodiments, the two or more regions of tissue are from a single organ. In some embodiments, the two or more regions of tissue are from different organs.
[0007] A roughness exponent may be generated via any appropriate multi-scale analytical method. In some embodiments, the roughness exponent is calculated using one or more multi-scale methods selected from a wavelet-transform modulus maxima, a wavelet leader, detrended fluctuation, and Fourier analysis. A roughness exponent may also be calculated using any mathematically related and/or similar multi-scale density fluctuation assessment method yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent. In some embodiments, a roughness exponent is calculated from two- dimensional data. In some embodiments, a roughness exponent is calculated from three dimensional data (e.g., 3D data cubes).
[0008] In some embodiments, the present invention additionally provides methods including providing a first view of a region of tissue; providing a second view of the region of tissue; calculating a first fractal dimension for the first view of the region of tissue; and calculating a second fractal dimension for the second view of the region of tissue; wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous. In some embodiments, provided methods further comprise treating the region of tissue if it is cancerous.
[0009] In some embodiments, the fractal zone is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
[0010] In some embodiments, the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8); first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8); second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2); third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
[0011] In some embodiments, the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9); first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9); second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1); third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
[0012] In some embodiments, the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7); first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7); second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3); third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
[0013] In some embodiments, provided methods further include steps of providing a third view of the region of tissue; and calculating a third fractal dimension for the third view of the region of tissue; wherein if the fractal dimension of at least one of the first fractal dimension, the second fractal dimension, and the third fractal dimension is in the fractal zone, the region of tissue is considered cancerous. In some embodiments, a fractal dimension is calculated from two-dimensional data. In some embodiments, a fractal dimension is calculated from three dimensional data (e.g., 3D data cubes).
[0014] The methods provided by the present invention may be performed on any tissue.
In some embodiments, the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
[0015] As used in this application, the terms "about" and "approximately" are used as equivalents. Any numerals used in this application with or without about/approximately are meant to cover any normal fluctuations appreciated by one of ordinary skill in the relevant art.
[0016] Other features, objects, and advantages of the present invention are apparent in the detailed description that follows. It should be understood, however, that the detailed description, while indicating embodiments of the present invention, is given by way of illustration only, not limitation. Various changes and modifications within the scope of the invention will become apparent to those skilled in the art from the detailed description.
DEFINITIONS
[0017] In order for the present invention to be more readily understood, certain terms are first defined below. Additional definitions for the following terms and other terms are set forth throughout the specification.
[0018] Animal: As used herein, the term "animal" refers to any member of the animal kingdom. In some embodiments, "animal" refers to humans, at any stage of development. In some embodiments, "animal" refers to non-human animals, at any stage of development. In certain embodiments, the non-human animal is a mammal {e.g., a rodent, a mouse, a rat, a rabbit, a monkey, a dog, a cat, a sheep, cattle, a primate, and/or a pig). In some embodiments, animals include, but are not limited to, mammals, birds, reptiles, amphibians, fish, insects, and/or worms. In some embodiments, an animal may be a transgenic animal, genetically-engineered animal, and/or a clone.
[0019] Approximately or about: As used herein, the term "approximately" or "about," as applied to one or more values of interest, refers to a value that is similar to a stated reference value. In certain embodiments, the term "approximately" or "about" refers to a range of values that fall within 25%, 20%, 19%, 18%, 17%, 16%, 15%, 14%, 13%, 12%, 11%, 10%, 9%, 8%, 7%), 6%), 5%, 4%, 3%, 2%, 1%, or less in either direction (greater than or less than) of the stated reference value unless otherwise stated or otherwise evident from the context (except where such number would exceed 100% of a possible value).
[0020] Cancer: As used herein, the term "cancer" refers to a group of diseases, all involving unregulated cell growth. Exemplary cancers include, without limitation: Acute lymphoblastic leukemia, Acute myeloid leukemia, Adrenocortical carcinoma; AIDS-related cancers; AIDS-related lymphoma; Anal cancer; Appendix cancer; Astrocytoma,childhood cerebellar or cerebral; Basal cell carcinoma; Bile duct cancer, extrahepatic; Bladder cancer; Bone cancer, Osteosarcoma/Malignant fibrous histiocytoma; Brainstem glioma; Brain tumor; Brain tumor, cerebellar astrocytoma; Brain tumor, cerebral astrocytoma/malignant glioma; Brain tumor, ependymoma; Brain tumor, medulloblastoma; Brain tumor, supratentorial primitive neuroectodermal tumors; Brain tumor, visual pathway and hypothalamic glioma; Breast cancer;
Bronchial adenomas/carcinoids; Burkitt lymphoma; Carcinoid tumor, childhood; Carcinoid tumor, gastrointestinal; Carcinoma of unknown primary; Central nervous system lymphoma, primary; Cerebellar astrocytoma, childhood; Cerebral astrocytoma/Malignant glioma, childhood; Cervical cancer; Childhood cancers; Chronic lymphocytic leukemia; Chronic myelogenous leukemia; Chronic myeloproliferative disorders; Colon Cancer; Cutaneous T-cell lymphoma; Desmoplastic small round cell tumor; Endometrial cancer; Ependymoma; Esophageal cancer; Ewing's sarcoma in the Ewing family of tumors; Extracranial germ cell tumor, Childhood;
Extragonadal Germ cell tumor; Extrahepatic bile duct cancer; Eye Cancer, Intraocular melanoma; Eye Cancer, Retinoblastoma; Gallbladder cancer; Gastric (Stomach) cancer;
Gastrointestinal Carcinoid Tumor; Gastrointestinal stromal tumor (GIST); Germ cell tumor: extracranial, extragonadal, or ovarian; Gestational trophoblastic tumor; Glioma of the brain stem; Glioma, Childhood Cerebral Astrocytoma; Glioma, Childhood Visual Pathway and
Hypothalamic; Gastric carcinoid; Hairy cell leukemia; Head and neck cancer; Heart cancer; Hepatocellular (liver) cancer; Hodgkin lymphoma; Hypopharyngeal cancer; Hypothalamic and visual pathway glioma, childhood; Intraocular Melanoma; Islet Cell Carcinoma (Endocrine Pancreas); Kaposi sarcoma; Kidney cancer (renal cell cancer); Laryngeal Cancer; Leukemias; Leukemia, acute lymphoblastic (also called acute lymphocytic leukemia); Leukemia, acute myeloid (also called acute myelogenous leukemia); Leukemia, chronic lymphocytic (also called chronic lymphocytic leukemia); Leukemia, chronic myelogenous (also called chronic myeloid leukemia); Leukemia, hairy cell; Lip and Oral Cavity Cancer; Liposarcoma; Liver Cancer (Primary); Lung Cancer, Non-Small Cell; Lung Cancer Small Cell Lymphomas; Lymphoma, Burkitt; Lymphoma, cutaneous T-Cell; Lymphoma, Hodgkin; Lymphoma, Primary Central Nervous System; Macroglobulinemia, Waldenstrom; Malignant Fibrous Histiocytoma of Bone/Osteosarcoma; Medulloblastoma, Childhood; Melanoma; Melanoma, Intraocular (Eye); Merkel Cell Carcinoma; Mesothelioma, Adult Malignant; Mesothelioma, Childhood; Metastatic Squamous Neck Cancer with Occult Primary; Mouth Cancer; Multiple Endocrine Neoplasia Syndrome, Childhood; Multiple Myeloma/Plasma Cell Neoplasm; Mycosis Fungoides;
Myelodysplasia Syndromes; Myelodysplastic/Myeloproliferative Diseases; Myelogenous Leukemia, Chronic; Myeloid Leukemia, Adult Acute; Myeloid Leukemia, Childhood Acute; Myeloma, Multiple (Cancer of the Bone-Marrow); Myeloproliferative Disorders, Chronic; Nasal
cavity and paranasal sinus cancer; Nasopharyngeal carcinoma; Neuroblastoma; Non-Hodgkin lymphoma; Non-small cell lung cancer; Oral Cancer; Oropharyngeal cancer;
Osteosarcoma/malignant fibrous histiocytoma of bone; Ovarian cancer; Ovarian epithelial cancer (Surface epithelial-stromal tumor); Ovarian germ cell tumor; Ovarian low malignant potential tumor; Pancreatic cancer; Pancreatic cancer, islet cell; Paranasal sinus and nasal cavity cancer; Parathyroid cancer; Penile cancer; Pharyngeal cancer; Pheochromocytoma; Pineal astrocytoma; Pineal germinoma; Pineoblastoma and supratentorial primitive neuroectodermal tumors, childhood; Pituitary adenoma; Plasma cell neoplasia/Multiple myeloma; Pleuropulmonary blastoma; Primary central nervous system lymphoma; Prostate cancer; Rectal cancer; Renal cell carcinoma (kidney cancer); Renal pelvis and ureter, transitional cell cancer; Retinoblastoma; Rhabdomyosarcoma, childhood; Salivary gland cancer; Sarcoma, Ewing family of tumors;
Sarcoma, Kaposi; Sarcoma, soft tissue; Sarcoma, uterine; Sezary syndrome; Skin cancer (nonmelanoma); Skin carcinoma, Merkel cell; Small intestine cancer; Soft tissue sarcoma;
Squamous cell carcinoma; Squamous neck cancer with occult primary, metastatic; Stomach cancer; Supratentorial primitive neuroectodermal tumor, childhood; T-Cell lymphoma, cutaneous; Testicular cancer; Throat cancer; Thymoma, childhood; Thymoma and Thymic carcinoma; Thyroid cancer; Thyroid cancer, childhood; Transitional cell cancer of the renal pelvis and ureter; Trophoblastic tumor, gestational; Unknown primary site, carcinoma of, adult; Unknown primary site, cancer of, childhood; Ureter and renal pelvis, transitional cell cancer; Urethral cancer; Uterine cancer, endometrial; Uterine sarcoma; Vaginal cancer; Visual pathway and hypothalamic glioma, childhood; Vulvar cancer; and Wilms tumor (kidney cancer), childhood.
[0021] Pre-Cancer: As used herein the term "pre-cancer" or "pre-cancerous" is used to refer to tissue in or from a subject that is not yet cancerous, but has a higher chance of becoming cancerous that normal tissue.
[0022] Risk: As will be understood from context, a "risk" of a disease, disorder, and/or condition comprises a likelihood that a particular individual will develop a disease, disorder, and/or condition {e.g., a cancer). In some embodiments, risk is expressed as a percentage. In some embodiments, risk is from 0,1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 20, 30, 40, 50, 60, 70, 80, 90 up to 100%. In some embodiments risk is expressed as a risk relative to a risk associated with a
reference sample or group of reference samples. In some embodiments, a reference sample or group of reference samples have a known risk of a disease, disorder, condition and/or event (e.g., a muscular dystrophy). In some embodiments a reference sample or group of reference samples are from individuals comparable to a particular individual. In some embodiments, relative risk is 0,1, 2, 3, 4, 5, 6, 7, 8, 9, 10, or more.
[0023] Subject: As used herein, the term "subject" refers to a human or any non-human animal (e.g., mouse, rat, rabbit, dog, cat, cattle, swine, sheep, horse or primate). A human includes pre and post natal forms. In many embodiments, a subject is a human being. A subject can be a patient, which refers to a human presenting to a medical provider for diagnosis or treatment of a disease. The term "subject" is used herein interchangeably with "individual" or "patient." A subject can be afflicted with or is susceptible to a disease or disorder but may or may not display symptoms of the disease or disorder.
[0024] Substantially: As used herein, the term "substantially" refers to the qualitative condition of exhibiting total or near-total extent or degree of a characteristic or property of interest. One of ordinary skill in the biological arts will understand that biological and chemical phenomena rarely, if ever, go to completion and/or proceed to completeness or achieve or avoid an absolute result. The term "substantially" is therefore used herein to capture the potential lack of completeness inherent in many biological and chemical phenomena.
BRIEF DESCRIPTION OF THE DRAWING
[0025] FIG. 1 shows a representative schematic of imaging breast tissue from
mediolateral oblique (MLO) and cranio-caudal (CC) angles.
[0026] FIG. 2 depicts a graph of the minimum fractal dimension by maximum fractal dimension, with tissue classified according to known methods as benign in blue, and tissue classified according to known methods as malignant tissue, in red.
[0027] FIG. 3 shows exemplary results of a cluster analysis of scatter plots obtained from graphing minimum fractal dimension by maximum fractal dimension of (a) benign tissue, and (b) malignant tissue.
[0028] FIG. 4 shows images of breast tissue classified as (a) benign and (b) malignant as classified by provided methods. Note the presence of a majority of disrupted tissue in the image of malignant tissue (represented by yellow circles) and the lack of disrupted tissue in the benign tissue (denoted by a lack of yellow circles).
[0029] FIG. 5 shows a boxplot of the distribution of ΔΗ for both benign and malignant
(cancer) tissues.
[0030] FIG. 6 shows exemplary photographs of breast tissue classified according to provided methods. Specifically, the left panel shows control breast tissue made up of healthy fatty tissue (H < 0.45); the center panel shows a breast with a benign lesion lying at the interface between fatty and dense tissue, showing a few surrounding sub-regions of disrupted tissue (yellow, 0.45 < H < 0.55); and the right panel shows a breast with a malignant lesion with several disrupted tissue neighboring regions.
[0031] FIG. 7 A shows a sample simulated fractional Brownian motion image BH=O.S(X);
7B shows the gradient of the image in 7A is obtained as the modulus of the wavelet transformed using equation 2.10 below; 7C shows maxima chains in blue defined as positions where the WT modulus is locally maximal (i.e. the WTMM) and along these chains in 7C further local maxima are found in red (i.e. the WTMMM); this is then repeated as several different scales, three of which are shown in 7D, 7E, and 7F; the WTMMM were then connected vertically through scales to define the WT skeleton shown in 7G, with the gray-scale coding going from black (minimum) to white (maximum).
[0032] FIG. 8A shows an original image obtained from the DDSM database, and 8B shows a zoomed in image with a suspicious region encircled by a radiologist; by selecting appropriate values for the slope of h of the WT modulus as a function of scale in a logarithmic representation, and the log of the pre-factor, log(k) in 8D, the WTMMM (blue) from the tissue background in 8C are distinguished from the WTMMM (red) that belong to the MC in 8E; from that point the WT skeleton can be calculated from the WTMMM that belong to the lesion from those that belong to the background tissue; the corresponding WTMM chains at the smallest scale are shown in 8F and 8G for the background and lesion, respectively.
[0033] FIG. 9A-G shows the same analysis as in FIG8A-G, only on a different case.
[0034] FIG. 10 shows the frequency distributions of fractal dimensions D calculated for benign cranio-caudal (CC) and mediolateral oblique (MLO) views (top panels) as compared to the fractal dimensions D calculated for cancer CC and MLO views (bottom panels); it is of note that the distributions are drastically different between benign and cancer samples.
[0035] FIG. 11 shows a fractal dimension plot of several samples; specifically, each case is represented by a single dot and is plotted with the fractal dimension obtained from the mediolateral oblique (MLO) view as a function of the fractal dimension obtained from the cranio-caudal (CC) view. A polygonal region is outlined, the inside of which is defined as the "fractal zone" while the outside is defined as the "Euclidean zone". The dots represent malignant (dark red) and benign (light green) cases.
[0036] FIG. 12 shows exemplary random point distribution models to describe how different kinds of objects can have only a limited number of possible fractal dimension as a function of the projection angle; the left panel (12A) shows the possible fractal dimensions of a line; the center panel (12B) shows the possible fractal dimensions of a surface, and the right panel (12C) shows the possible fractal dimensions of a fractal cluster (here a simulated diffusion- limited aggregate).
[0037] FIG. 13 shows a block diagram of an exemplary cloud computing environment suitable for use with provided methods.
[0038] FIG. 14 is a block diagram of an exemplary computing device and a mobile computing device suitable for use with provided methods.
DETAILED DESCRIPTION OF CERTAIN EMBODIMENTS
[0039] The present invention provides, among other things, methods and systems for identifying cancerous or pre-cancerous tissue in a subject. The present invention is based, in part, on the surprising discovery that analysis of tissue surrounding a region of interest rather than tissue from the region of interest itself, may be used to determine if a particular region of tissue is cancerous or pre-cancerous. In some embodiments, the present invention provides, among other things, methods of identifying cancerous or pre-cancerous tissue including providing a first region of tissue from a subject, calculating a roughness exponent for the first
region of tissue, and comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 (e.g., less than 0.1, less than 0.05) between the roughness exponent of the first region of tissue and 0.5 indicates that the tissue is cancerous or precancerous. In some embodiments, provided methods further comprise treating the tissue adjacent to the first region of tissue if it is cancerous or pre-cancerous.
[0040] In some embodiments, the present invention also provides methods including providing a first view of a region of tissue, providing a second view of a region of tissue, calculating a first fractal dimension for the first view of the region of tissue, and calculating a second fractal dimension for the second view of the region of tissue, wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
[0041] In some embodiments, it may be advantageous to perform the providing, calculating, and/or comparing steps of provided methods multiple times. In some embodiments, it is contemplated that performing the providing, calculating and comparing steps multiple times will allow for superior characterization of the tissue or the tissue surrounding a tissue region of interest as compared to performing each step a single time according to various embodiments. In some embodiments, the providing, calculating, and comparing steps are each performed a plurality of times. In some embodiments, the plurality of times is at least 10, at least 20, at least 50, at least 100, at least 300, at least 500, or at least 1,000 times. In some embodiments, performing the providing, calculating, and comparing steps multiple times occurs for different regions of the same organ or tissue. In some embodiments, performing the providing, calculating, and comparing steps multiple times occurs for substantially the same region(s) of the same organ or tissue. In some embodiments, performing the providing, calculating, and comparing steps multiple times occurs for different organs or tissues.
Prior Detection Methods
[0042] Breast cancer is the most common cancer worldwide according to the World
Health Organization (WHO) and the second leading cause of cancer related death among women in the United States. Despite the recent advances in the medical field, the breast cancer rate has
continued to increase over the last 30 years. Cancer is easiest to treat when it is found in the early stages of development, making it critical for women to have regular screenings as recommended by the American Cancer Society (ACS). Mammograms are currently one of the most accepted screening processes and have been proven to be successful in detecting microcalcifications (MC), which are small deposits of calcium in breast tissue (about 200 microns in size), and can be an indicator of the early development of breast cancer.
[0043] By identifying early breast cancer, the survival rate for the patient increases.
Thus, early detection is key for the patient. However, detecting tumors in the early stages of development may not always be easy due to complex tissue composition and the small size of the tumor. Therefore, it's important for radiologists to identify any asymmetric changes in tissue composition in the breast, as it could be a sign of a potential tumor environment. Doctors rely heavily on the use of mammography, especially for screening of older women, making it a very widely accepted form of breast cancer detection.
[0044] Although mammograms are currently the most effective way of detecting breast cancer, it remains difficult to interpret mammograms due to high tissue variability and 3D to 2D projection effects. Current practice is to have two expert radiologists read the mammograms to reduce interpretation errors. However, it is not always possible to have access to two
radiologist's interpretation due the size of the hospital or the cost. Computer Aided Diagnosis (CAD) systems were developed to assist doctors in analyzing medical images such as mammograms and since the approval of CAD by the Food and Drug Administration (FDA) in 1998 there has been much attention to developing such software. However, many of the current CAD methods rely on tissue homogeneity characterized by the Hurst exponent, meaning the tissue is assumed to be monofractal and uncorrected. Healthy fatty and dense tissue, shown by dark and light intensities respectively on mammograms, demonstrate how tissue may not be statistically homogeneous and there are many fluctuations of intensity throughout the image. By neglecting to accurately characterize the environment of possible breast lesions, important information is lost. Therefore, many CAD methods are not offering the expected performance, which leads to increase recall rates and can cause false-positives on up to 70 % of normal cases, resulting in an increase of unnecessary stress on women.
[0045] A key insight to developing a successful CAD is given by recent research which has shown the importance the micro environment can have on the growth and development of tumors. The environment of cancerous lesions has shown to be able to create a niche around the stem cells, favoring the survival of cancerous stem cells and protecting cells from any treatment or therapy. Thus, differences between normal and cancer stem cells and their interactions with the neighboring tissue may exist and be able to be detected by a CAD at a much larger scale. Therefore, it is critical for the neighboring tissue to be carefully examined and taken into consideration when using or developing a CAD to evaluate mammograms, something not done prior to the present invention.
[0046] One limitation of existing CAD is the inability of the software to rate the significance of the findings, similar to the Breast Imaging-Reporting and Data-System
(BI-RADS) assessment score given by the radiologist at the time of the mammogram
interpretation. As described herein, the two scores provided by the 2D Wavelet-Transform Modulus Maxima (WTMM) method, i.e. the fractal dimension of the MC and the roughness of the tissue surrounding the breast lesion given by the Hurst exponent, have known physical properties and can provide an insight to the invasiveness and severity of the possible lesion. Indeed, one of the several advantages provided by the present invention is the recognition that by examining not only the tumor, but also the microenvironment, important aspects of benign and malignant processes can be evaluated and possibly lead to detection of early cancer or even precancerous tissue.
Characterization of Dense and Fatty Tissues
[0047] An important part of interpreting mammograms relies on the ability of radiologists to identify the composition of tissue determined by dense and fatty components. Genetics play a role in determining breast density and radiologists can detect dense tissue as light intensities on an x-ray and fatty tissue as dark intensities. Density decreases with age as a normal process. To help provide radiologists with a uniform scoring system, the American College of Radiology developed an index which ranks density from 1 to 4, ranging from fatty to dense as shown in Table 1.
Table 1 - Breast Density Score
[0048] By taking into account the tissue composition of the breast, radiologists can more accurately identify suspicious regions. Since tissue composition should be symmetric between both breasts, differing densities maybe a potential sight for a tumor, causing
radiologists to pay close attention to these differing areas.
[0049] There have been many statistical studies devoted to mammography analysis by using fractal techniques. One such study was conducted by P. Kestener and colleagues using the 2D WTMM method to analyze normal mammary parenchyma. The goal of the Kestener et al. study was to accurately classify mammographic tissue as dense or fatty. There, the images were obtained from the Digital Database for Screening Mammography (DDSM) at the
University of South Florida. The databank contains over 2,500 studies made up of normal, benign and malignant mammograms all categorized by an expert radiologist. Each study has two images of each breast, consisting of a mediolateral oblique (MLO) view and cranio- caudal (CC) view with any suspicious region circled by a radiologist. The suspicious region could contain a mass and/or micro calcifications (MC), but only the cases that were classified as normal, i.e. no suspicious area was identified, were looked at in the Kestener et al. study.
[0050] After Kestener analyzed 10 images, 5 fatty breasts and 5 dense breasts, from the
DDSM website using the 2D WTMM method, their analysis characterized fatty tissue as H = 0.25 ±0.05 and dense tissue as H = 0.65 ±0.05.
[0051] Both fatty and dense tissue display monofractal scaling behavior, with fatty tissue, H = [0.20 - 0.35] beingthe signature of anti-persistent roughness fluctuations and dense tissue,H = [0.55 - 0.75], the signature of persistent long-range correlations. Note that no tissue classification index exists in H = (0.35 - 0.65) and healthy breast tissue is composed of only fatty or dense tissue.
Application of Wavelet-transform Modulus Maxima to Subject Tissue
[0052] Various embodiments of the present invention apply wavelet-transform modulus maxima (WTMM) to one or more images of a tissue in order to characterize one or more attributes of the tissue and/or tissue adjacent to the tissue (i.e., the tissue's microenvironment). The present invention evidences the potential of the two dimensional WTMM method to become a powerful tool in interpreting mammograms. The WTMM method has proven to be successful in several fields of applied science, including geology, astrophysics, cellular biology and orthopedic medicine. The method was originally developed as a multifractal formalism to analyze highly complex ID signals, 2D images, and 3D images. As described in detail herein, the wavelet transform (WT) acts as a mathematical microscope to characterize spatial image information over a continuous range of size scales. It is the gradient vector of the image smoothed by dilated versions of a Gaussian filter. At each size scale, the wavelet transform modulus maxima (WTMM) are defined by the positions where the modulus of the WT is locally maximal. These WTMM are automatically organized as maxima chains at the considered scale. Along each of these chains, further local maxima are found, the WTMM maxima (WTMMM). This process is repeated for all size scales and the WTMMM from each scale are then linked to form the WT skeleton.
[0053] The ability to consider vertical lines in the WT skeleton individually is significant, since it allows one to objectively discriminate between lines pointing to the tissue background from those pointing to the lesion by considering how the WT modulus varies as a function of the scale parameter along each line. One can then calculate the so-called singularity spectrum separately for each subset, which then allows consideration of the roughness exponent
H, characterizing the tissue background, and the fractal dimension D of the lesion, characterizing the architecture of the lesion, and their microenvironment.
[0054] In some embodiments, provided methods make use of an adaptation of a two dimensional WTMM method, specifically, the continuous wavelet transform method, as a mathematical microscope used to characterize the fractal geometry of clusters of
microcalcifications (MC) in human tissue (e.g., breast tissue) and/or to determine the roughness of the background tissue (i.e. tissue adjacent to the tissue of interest) seen in mammograms. The WTMM method yields the so-called singularity spectrum, D(h), i.e. the fractal dimension D, of points having a Holder exponent of h. The MC are seen as Dirac singularities by the WTMM method, therefore having Holder exponent value of h = -1. This allows the MC with h ~ -1 to be abstracted from the background tissue which has h ~ 0.30 for fatty breast tissue and h ~ 0.65 for dense tissue. Thus, the WTMM method is used to perform a segmentation of tissues, in some embodiments, the breast tumor tissue, based on the strength of the singularities composing the mammogram images, and to simultaneously quantify their fractal dimension.
[0055] Comparing the results of provided methods applied on several hundred images from a digital databank of mammograms with known radiologist diagnostics, the fractal dimensions of benign and malignant breast lesions are significantly different, with benign having an integer dimension corresponding to a non-invasive Euclidean object and cancer having a non- integer dimension, representing an invasive structure. In addition, the microenvironments characterized by the roughness of the tissue in which the lesions are embedded (i.e., adjacent to) are different for benign and malignant tumors, and provides an insight into the onset and development of breast cancer.
[0056] One score provided by the two dimensional WTMM method is the roughness of the tissue surrounding the breast lesion, given by the Hurst exponent, which has known physical properties and can provide an insight to the invasiveness and severity of the possible lesion. By examining not only the tumor but also the microenvironment, important aspects of benign and malignant processes can be evaluated and possibly lead to early detection of cancer or even precancerous tissue. In various embodiments, other multi-scale methods may also be used as an alternative to WTMM including, but not limited to, wavelets leaders, detrended fluctuation analysis, Fourier transform methodologies, and/or any mathematically related and/or similar
multi-scale density fluctuation assessment methods yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent.
[0057] The Hurst exponent (or equivalent roughness exponents obtained by similar multi-scale analytical methods) represents, in part, the physical status of the microenvironment of the tissue (e.g., tissue density as determined from roughness fluctuations) and the level of correlation in that physical status across adjacent tissue. In general, if H < 0.5, then the roughness fluctuation are considered anti-correlated, while H > 0.5 means that the roughness fluctuations are positively correlated. In either case, the physical system of the analyzed tissue is considered to have some form of spatial memory. An H = 0.5 means that the roughness fluctuations are uncorrected and are thus considered to have no spatial memory.
[0058] The present invention encompasses the recognition that an H = 0.5, or an H value within some range surrounding 0.5, in some embodiments between about 0.3-0.65, between about 0.35-0.6, between about 0.4-0.6, between about 0.4-0.55, or between about 0.45 - 0.55, indicates a region of tissue that is potentially cancerous or pre-cancerous. As described above, around H = 0.5 the system has a reduced or loss of correlation between roughness fluctuations, indicating a breakdown in order in the system. It is the recognition that this disorder likely indicates a cancerous or pre-cancerous state that is an important aspect of the methods and systems of the present invention.
[0059] The present invention recognizes that the results obtained by using the WTMM method on screening mammograms could not only be used as a possible computer-aided detection (CAD) method, but also as a method to further study the biophysics of tumor onset and progression. As described in the Examples below, after all images were analyzed using provided methods, a statistical analysis was performed on all data, which provided us with information on the critical differences between the organization, behavior, and biological and physical processes of both benign and malignant tumors.
Exemplary Mechanics of Characterizing the Local Regularity of Rough Surfaces with the Wavelet Transform Modulus Maxima Method
[0060] Most of the fractal methods used to analyze mammograms rely on the estimate of the fractal dimension which is related to the Hurst exponent H which statistically characterizes the global roughness of the mammogram landscape. The two dimensional WTMM method accounts for possible fluctuations of the local regularity of a rough surface as defined by the Holder exponent, h of the function / whose graph defines the rough surface under study. The 2D WTMM method provides a way to estimate the D(h) singularity spectrum which provides us with the Hausdorff dimension of the set of points x where the local roughness of the exponent h(x) is h.
[0061] This methodology has the ability to be applied to a variety of rough surfaces, such as mammogram landscapes. We will use the term "rough surface" for an irregular surface. This means the surface can be accurately described by a single -valued, self-affme function satisfying the following:
Vjx — (;CQ, Ϊ/Ο) E S2. Vx— { x, y ) € 2 ff e E
in the neighborhood of xo, such that, for any λ > 0: f {xQ ÷ ι/ο ÷ Xay) - / (;rn . ¾¾ ) ~ λΗ [f(x0 ÷ χ . 0 ÷ y)— /(χ 0, ?/0)| . ( 2.1)
[0062] In various embodiments, the Hurst exponent is used to characterize the global roughness of the function under investigation. If H < 1 , then the function, f, is nowhere
differentiable and the smaller the value for H, the more singular and irregular f. Several methods have been used to estimate the Hurst exponent of self-affme functions, but one has to be careful due to the fact that some functions, such as fractal functions, generally display multi- affine properties in the sense that their roughness fluctuates from point to point. To describe these multifractal functions, one needs to change the definition of the Hurst exponent off to become a local value, h(xo). Here, h is called the Holder exponent and it provides the strength of the singularity off at xo. A rigorous definition of the Holder exponent is given by the largest exponent h(x0) such that there exists a polynomial of degree n < h(x0) and a constant C > 0, such that for any point x in the neighborhood of xo one has:
(x) — Pn(x— X()) | < C|x /;.(x{) }
— Xr ( 2.2 '
[0063] If is n times continuously differentiable at the point xo, then one can use the order n Taylor series of/ at xo for the polynomial Pn(x - xO) and prove h(xo) > n. Thus h(x0) measures how irregular the function f is.
Multi-Scale Edge Detection
[0064] The edges of the different structures that appear in an image are often the most important features for pattern recognition. Therefore, many computer tools for edge detection look for points where the gradient of the image intensity has a modulus locally maximum in its direction. With an appropriate choice of analyzing wavelets, one can redefine the Canny's multi- scale edge detection in terms of a 2D wavelet transform.
[0065] In several embodiments, use of a multi-scale edge detection methodology is desirable to assist in determining the borders of cancerous or pre-cancerous tissue. Below is a brief discussion of exemplary multi-edge detection techniques that may be used in accordance with the present invention.
[0066] First, define two wavelets that are the partial derivatives with respect to x and y, respectively: dOix. tf ) ΰθίχ, ν )
V i xxu ) =——— , -VO i ^ v) =— x , " ' , ( 2.3 s where 9{x, y) is a 2D smoothing function that is well localized (around x = y = 0) and isotropic.
For any function/ (x, y) that exists in L (R) the wavelet transform has two components with respect to ψ\ and ψ2 and can be expressed in vectorial form as:
ίτ T^, [f] = or2 f ί χ¾(α_1( - b))/(x)
T,;, /j{b, a.)— (2,4;
After performing an integration by arts, one can obtain
2.7)
Wo; will take <vix) to be the Gaussian function:
[0067] If o(x) is simply a smoothing function like the Gaussian, then equation 2.7 amounts to define the 2D wavelet transform as the gradient vector of f(x) smoothed by dilated versions of ø (a_1x) of this filter.
[0068] We will be referring to the wavelet-transform in terms of the modulus and argument:
T^[/](b?«) - (M {f}(b,a),A {f]{b, )). (2.9^
Ablf}( t ) = i h[f](h, ) 4- iT,h[f](h, )) 12.11 )
[0069] As originally proposed by Mallat and collaborators, a very efficient way to
perform point- wise regularity analysis is to examine the wavelet transform modulus maxima. At any given scale a, the wavelet transform modulus maxima (WTMM) are defined as the points b where the wavelet transform modulus, M [f](b, a), is locally
maximum along the gradient direction given by the wavelet transform argument,
A [f](b, a). These WTMM, are inflection points of f *0a(x) and lie on connected chains called maxima chains. One only needs to record the position of the local maxima of M along the maxima chains and the values of M [fj and A [fj at the corresponding location. At each scale a, the wavelet analysis reduces to record the WTMM maxima (WTMMM) only. They indicate locally the direction where the signal has the sharpest variation and are disposed along connected curves across scales named maxima lines. The wavelet transform skeleton will be defined as the set of maxima lines that converge to the (x; y)-plane in the limit a 0+.
Discussion of Two Dimensional WTMM Method Mechanics
[0070] The 2D WTMM method relies upon the space-scale partitioning given by the WT skeleton. Let us define L(a) as the set of maxima lines that exist at the scale a and which contain maxima at any scale a '< a. The WTMM method consists in defining the following partition functions directly from the WTMMM that belong to the wavelet transform skeleton:
where q ER. Compared to the classical box-counting techniques, the analyzing wavelet, ψ, plays the role of a generalized "oscillating box", the scale, a, defines its size, while the WT skeleton indicates how to position our boxes to obtain a partition at the considered scale. From the analogy that links the multifractal formalism to thermodynamics, one can define the exponent q) from the power law behavior of the partition function:
[0071] Here, q and r( ) are the inverse temperature and the free energy respectively.
This formalism replaces these variables with the Holder exponent, A, and the singularity spectrum, D(h). This means that the D(h) singularity spectrum off can be determined from the Legendre transform of the partition function scaling exponent q):
D[ h) = mm{qh— r(q} ) (2.14)
q
[0072] From the properties of the Legendre transform, the homogeneous fractal functions that involve singularities of unique Holder exponent h =
are characterized by x{q)
spectrum which is a linear function of q. A nonlinear q) curve is the signature of
nonhomogeneous functions that display multifractal properties in the sense that the Holder exponent h(x) is a fluctuating quantity that depends upon the spatial position x.
[0073] The computation of the D(h) singularity spectrum, by the Legendre transform, requires first a smoothing of the r(q) curve. This smoothing process may loss any
interesting physics of phase transitions in the scaling properties of fractal functions. One can avoid directly performing the Legendre transform by considering the quantities h and D(h) as mean quantities defined with respect to their Boltzmann weights computed from the
WTMMM: sup,-., .,, -r■,'^ -, -Μ Ι f ifx. a' } \q
πν ί Π α, L. α ) = — - . < 2. l where Ziq. is the partition function. Th n one computes the expectation values:
h(q, a) = ln j sup .M.J/Kx. a ) \ W^ [f l {q, £. a) , (2. 10] tx .''}Ei
and
Diq, a) = V
C, a) liiiWU/lfV/- a)) ( 2.17) from which one extracts h{q) = iim h {q . a} / in a. (2.18) s-H>÷
and therefore the l i s ) singularity spectrum.
[0074] It is important to note that while several provided methods use a two dimensional analysis, it is contemplated as within the scope of the present invention that three dimensional analyses may also be accomplished. For example, provided methods may be used to analyze a three dimensional image generated via tomosynthesis, fluorescence microscopy, computed tomography, magnetic resonance imaging, or other digital methodology.
[0075] In some embodiments, the roughness analysis of tissue is not performed with two dimensional images. The WTMM method as well as all other similar space-scale techniques that allow the calculation of the Hurst exponent (or any analogous exponent), such as, but not limited to wavelet leaders, Fourier techniques, etc are also directly applicable to 3D data cubes.
Therefore, the assessment of cancerous or pre-cancerous tissue regions through a roughness analysis is applicable to 3D data. For 3D data the distance between the estimated Hurst exponent (or other analogous exponent) and the critical value 0.5 is the same as for the 2D case. The strategy of analyzing square sub-regions neighboring a tumor or suspicious tissue area such as presented in FIG. 6 is directly generalizable to the analysis of cubic sub-regions in 3D.
Detecting Microcalcifications Through Wavelet Transform Skeleton Segmentation
[0076] In addition to determination of roughness exponents, the present invention provides methods and systems for detecting microcalcifications (MC) in subject tissue and characterizing their fractal dimension. In some embodiments, provided methods use two
dimensional (or, in some embodiments, three dimensional) WTMM to provide this information, though any mathematically related and/or similar multi-scale density fluctuation assessment method yielding an exponent or spectrum of exponents that is similar and/or can complement the use of the Hurst exponent and/or fractal dimension may be used.
[0077] The presence of MC is one of the most important and sometimes the only sign of cancer, for example, breast cancer in a screening mammogram. MCs typically appear as bright spots on mammograms and the number, size, type (Table 2) and distribution (Table 3) provide the radiologist with information about the potential malignancy of the lesion and allow an assessment score (Table 4) to be assigned to each case.
Table 2 - Calcification Types
Table 3 - Calcification Distributions
tissue
Table 4 - BI-RADS Assessment Score
[0078] By using the information listed in the above tables, radiologists can predict the potential malignancy of the identified tumor on the mammogram. For example, tumors that contain pleomorphic or fine linear branching calcification types have a higher risk of being malignant than tumors that are classified as punctate. However, it is not uncommon for tumors to be made up of components that make it difficult for radiologists to the identify the tumor as cancerous or a benign disease such as fibrocystic breast condition or a blunt injury. One benefit of a CAD is the objectivity of the software of having a range of indices associated with benign and/or malignant tumors.
Identifying Microcalcifications (MC)
[0079] In some embodiments, the present invention provides computer aided diagnostic tools, using the WT methodology to detect MC by inspecting the WT maxima chains. At the smallest scale resolved by the WT microscope, 7 pixels, MC which can be considered at strong singularities are contour-shaped by some maxima chains. Since the average size of MC is about 200 μτη, or 5 pixels, these singularities are identified by the mathematical microscope as Dirac singularities. Thus, the corresponding maxima lines pointing to the MC are likely to display scaling properties with a Holder exponent h = -1 down to scales of the order of the MC size where one can observe a cross-over to the value h = 0 as the signature of the MC boundary. Therefore, one can perform a classification of these lines according to the behavior of the Μψ [f ] and then separate MC, h = -1, from dense tissue, h ~ 0.65, and fatty tissue, h ~ 0.30.
[0080] As pointed out above, the MC wavelet-transform skeleton can be used to compute the partition functions, and thus fully characterize the fractal geometry of the MC cluster.
[0081] Separating the MC from the subject tissue (e.g., breast tissue) in this way is key, since it allows one to consider the properties of the background tissue separate from the properties of the MC. Using provided methods, one is able to quantify the arrangement and organization of the MC and consider the distribution of the tumor, making the similar
information currently used in Table 3 less subjective. At the same time, one can classify the tissue and analyze any interesting physical or structural changes from previous mammograms or surrounding tissue. Therefore, the two dimensional WTMM method has the ability to be a powerful CAD. The information obtained from the wavelet-transform analysis can help quantify data that would otherwise be considered subjective. In some embodiments, provided methods use this information to aid in the development of an objective mammographic classification, which can be used to judge the potential malignancy of a radiologist-identified suspicious region.
Assessment of Three Dimensional Fractal Geometry from Multiple Two Dimensional Views
[0082] In some embodiments, characterization of the fractal geometry of MC clusters from a plurality of separate two-dimensional views of the same region of tissue (e.g. a breast),
allows for an assessment of the overall fractal dimension of the tissue region of interest, such as through application of the WT methods described herein. In other words, according to various embodiments, by using a plurality of separate two-dimensional views of the same region of tissue one can assess the fractal dimension of the three dimensional tissue region of interest. These concepts are described in more detail in the Examples below.
[0083] In some embodiments, non-cancerous or non-pre-cancerous tissue will have a
Euclidean fractal dimension (i.e. an integer fractal dimension). In some embodiments, a cancerous or pre-cancerous tissue will have a non-integer fractal dimension indicating an invasive morphology.
[0084] In some embodiments, the present invention also provides methods including providing a first view of a region of tissue, providing a second view of a region of tissue, calculating a first fractal dimension for the first view of the region of tissue, and calculating a second fractal dimension for the second view of the region of tissue, wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous. In some embodiments, provided methods further comprise treating the region of tissue if it is cancerous or pre-cancerous.
[0085] In some embodiments, the "fractal zone" is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
[0086] In order to provide greater clarity around the concept of fractal zones, several non-limiting exemplary embodiments are provided. In some embodiments, such as that shown in FIG. 11, a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8); first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8);
second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2); third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
[0087] In some embodiments, a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9); first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9); second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1); third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
[0088] In some embodiments, a fractal zone may be defined as consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue with the following dimensions: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7); first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7); second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3); third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
[0089] It is contemplated that the precise boundaries of the fractal zone may vary according to tissue type and/or tumor type (or subtype). Additionally, without wishing to be held to a particular theory it is possible that the precise borders of the fractal zone may change over time as tissue changes from normal tissue to pre-cancerous or cancerous tissue. One of skill in the art will be able to ascertain from the present disclosure how to determine the appropriate
boundaries of the fractal zone in a particular clinical scenario using no more than routine experimentation.
[0090] In some embodiments, provided methods may be applied to more than two views of a tissue region of interest. In some embodiments, provided methods may be applied to three, four, five, six, seven, eight, nine, ten, fifteen, twenty or more views of a particular tissue region of interest. One of skill will understand how the addition of additional two-dimension images to the analysis of provided methods will alter the definition of the fractal zone. For example, if three views of a tissue region of interest are used, the fractal zone may then be defined as the union of seven juxtaposed regions: one central cube and six extending pyramidal or conical regions. In some embodiments, wherein four or more views of a tissue region of interest are used, it may be more convenient to define the fractal zone purely mathematically, rather than in terms of geometric shapes. Additionally, the borders of the fractal zone, though described herein (including in the below Examples) as lines, may comprise other forms and definitions. For example, in some embodiments, the borders of the fractal zone may be smooth curves, jagged lines, or even fractals.
[0091] In addition to the use of the WTMM method to calculate the fractal dimension of clusters of microcalcifications , other methods may also be used according to various embodiments. In some embodiments, any fractal technique or techniques yielding a fractal dimension, such as, but not limited to: box-counting techniques, the perimeter-area relationship, packing dimension, Hausdorff dimension, capacity dimension, correlation dimension, the generalized dimensions (multifractal), and/or any space-scale technique that yields a value quantifying the structural complexity of a tumor or object of interest as a function of scale may be used.
[0092] As described elsewhere herein, provided methods may be applied to any tissue that may contain or comprise cancer. Non-limiting exemplary such tissues include breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
Treatment of Cancerous or Pre-Cancerous Tissue
[0093] According to various embodiments, provided methods may be used to detect and/or characterize the presence of cancerous or pre-cancerous tissue in a subject. It is contemplated that if cancerous or pre-cancerous tissue is discovered, that treatment of such tissue would commence in accordance with sound medical judgment. It is further contemplated that treatment may occur with any anti-cancer therapy, whether currently known or discovered in the future.
[0094] While any anti-cancer therapy may be appropriate for use in some embodiments, exemplary types of anti-cancer therapy are described below in order to better illustrate some of the principles of the present invention.
[0095] Traditional therapies or anticancer agents include surgery, radiotherapy (γ- radiation, neutron beam radiotherapy, electron beam radiotherapy, proton therapy,
brachytherapy, and systemic radioactive isotopes, to name a few), endocrine therapy, biologic response modifiers (interferons, interleukins, and tumor necrosis factor (TNF) to name a few), hyperthermia and cryotherapy, agents to attenuate any adverse effects (e.g., antiemetics), and other approved chemotherapeutic drugs, including, but not limited to, alkylating drugs
(mechlorethamine, chlorambucil, Cyclophosphamide, Melphalan, Ifosfamide), antimetabolites (Methotrexate), purine antagonists and pyrimidine antagonists (6-Mercaptopurine, 5- Fluorouracil, Cytarabile, Gemcitabine), spindle poisons (Vinblastine, Vincristine, Vinorelbine, Paclitaxel), podophyllotoxins (Etoposide, Irinotecan, Topotecan), antibiotics (Doxorubicin, Bleomycin, Mitomycin), nitrosoureas (Carmustine, Lomustine), inorganic ions (Cisplatin, Carboplatin), enzymes (Asparaginase), and hormones (Tamoxifen, Leuprolide, Flutamide, and Megestrol), to name a few. Any and all of these exemplary therapies, as well as other cancer therapies known by one of skill in the art, may be used in connection with the present invention.
[0096] It is contemplated that the specific treatment, including dose, dosing regimen, mode of administration, and timing of the onset and termination of therapy will be determined by a medical practitioner in accordance with sound medical judgment.
Systems and Devices for Implementing Provided Methods
[0097] The present invention also provides systems and devices, such as computing devices, implementing provided methods. In some embodiments, a computing device maybe a cloud computing device, a mobile computing device, or any other application-appropriate computing device.
[0098] As shown in FIG. 13, an exemplary implementation of a network environment
1300 for use with provided methods is shown and described. In brief overview, referring now to FIG. 13, a block diagram of an exemplary cloud computing environment 1300 is shown and described. The cloud computing environment 1300 may include one or more resource providers 1302a, 1302b, 1302c (collectively, 1302). Each resource provider 1302 may include computing resources. In some implementations, computing resources may include any hardware and/or software used to process data. For example, computing resources may include hardware and/or software capable of executing algorithms, computer programs, and/or computer applications. In some implementations, exemplary computing resources may include application servers and/or databases with storage and retrieval capabilities. Each resource provider 1302 may be connected to any other resource provider 1302 in the cloud computing environment 1300. In some implementations, the resource providers 1302 may be connected over a computer network 1308. Each resource provider 1302 may be connected to one or more computing device 1304a, 1304b, 1304c (collectively, 1304), over the computer network 1308.
[0099] The cloud computing environment 1300 may include a resource manager 1306.
The resource manager 1306 may be connected to the resource providers 1302 and the computing devices 1304 over the computer network 1308. In some implementations, the resource manager
1306 may facilitate the provision of computing resources by one or more resource providers
1302 to one or more computing devices 1304. The resource manager 1306 may receive a request for a computing resource from a particular computing device 1304. The resource manager 1306 may identify one or more resource providers 1302 capable of providing the computing resource requested by the computing device 1304. The resource manager 1306 may select a resource provider 1302 to provide the computing resource. The resource manager 1306 may facilitate a connection between the resource provider 1302 and a particular computing device 1304. In some implementations, the resource manager 1306 may establish a connection between a
particular resource provider 1302 and a particular computing device 1304. In some implementations, the resource manager 1306 may redirect a particular computing device 1304 to a particular resource provider 1302 with the requested computing resource.
[0100] FIG. 14 shows an example of a computing device 1400 and a mobile computing device 1450 that can be used to implement the techniques described in this disclosure. The computing device 1400 is intended to represent various forms of digital computers, such as laptops, desktops, workstations, personal digital assistants, servers, blade servers, mainframes, and other appropriate computers. The mobile computing device 1450 is intended to represent various forms of mobile devices, such as personal digital assistants, cellular telephones, smart- phones, and other similar computing devices. The components shown here, their connections and relationships, and their functions, are meant to be examples only, and are not meant to be limiting.
[0101] The computing device 1400 includes a processor 1402, a memory 1404, a storage device 1406, a high-speed interface 1408 connecting to the memory 1404 and multiple highspeed expansion ports 1410, and a low-speed interface 1412 connecting to a low-speed expansion port 1414 and the storage device 1406. Each of the processor 1402, the memory 1404, the storage device 1406, the high-speed interface 1408, the high-speed expansion ports 1410, and the low-speed interface 1412, are interconnected using various busses, and may be mounted on a common motherboard or in other manners as appropriate. The processor 1402 can process instructions for execution within the computing device 1400, including instructions stored in the memory 1404 or on the storage device 1406 to display graphical information for a GUI on an external input/output device, such as a display 1416 coupled to the high-speed interface 1408. In other implementations, multiple processors and/or multiple buses may be used, as appropriate, along with multiple memories and types of memory. Also, multiple computing devices may be connected, with each device providing portions of the necessary operations (e.g., as a server bank, a group of blade servers, or a multi-processor system).
[0102] The memory 1404 stores information within the computing device 1400. In some implementations, the memory 1404 is a volatile memory unit or units. In some implementations, the memory 1404 is a non-volatile memory unit or units. The memory 1404 may also be another form of computer-readable medium, such as a magnetic or optical disk.
[0103] The storage device 1406 is capable of providing mass storage for the computing device 1400. In some implementations, the storage device 1406 may be or contain a computer- readable medium, such as a floppy disk device, a hard disk device, an optical disk device, or a tape device, a flash memory or other similar solid state memory device, or an array of devices, including devices in a storage area network or other configurations. Instructions can be stored in an information carrier. The instructions, when executed by one or more processing devices (for example, processor 1402), perform one or more methods, such as those described above. The instructions can also be stored by one or more storage devices such as computer- or machine- readable mediums (for example, the memory 1404, the storage device 1406, or memory on the processor 1402).
[0104] The high-speed interface 1408 manages bandwidth-intensive operations for the computing device 1400, while the low-speed interface 1412 manages lower bandwidth-intensive operations. Such allocation of functions is an example only. In some implementations, the highspeed interface 1408 is coupled to the memory 1404, the display 1416 (e.g., through a graphics processor or accelerator), and to the high-speed expansion ports 1410, which may accept various expansion cards (not shown). In the implementation, the low-speed interface 1412 is coupled to the storage device 1406 and the low-speed expansion port 1414. The low-speed expansion port 1414, which may include various communication ports (e.g., USB, Bluetooth®, Ethernet, wireless Ethernet) may be coupled to one or more input/output devices, such as a keyboard, a pointing device, a scanner, or a networking device such as a switch or router, e.g., through a network adapter.
[0105] The computing device 1400 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a standard server 1420, or multiple times in a group of such servers. In addition, it may be implemented in a personal computer such as a laptop computer 1422. It may also be implemented as part of a rack server system 1424. Alternatively, components from the computing device 1400 may be combined with other components in a mobile device (not shown), such as a mobile computing device 1450. Each of such devices may contain one or more of the computing device 1400 and the mobile computing device 1450, and an entire system may be made up of multiple computing devices communicating with each other.
[0106] The mobile computing device 1450 includes a processor 1452, a memory 1464, an input/output device such as a display 1454, a communication interface 1466, and a transceiver 1468, among other components. The mobile computing device 1450 may also be provided with a storage device, such as a micro-drive or other device, to provide additional storage. Each of the processor 1452, the memory 1464, the display 1454, the communication interface 1466, and the transceiver 1468, are interconnected using various buses, and several of the components may be mounted on a common motherboard or in other manners as appropriate.
[0107] The processor 1452 can execute instructions within the mobile computing device
1450, including instructions stored in the memory 1464. The processor 1452 may be
implemented as a chipset of chips that include separate and multiple analog and digital processors. The processor 1452 may provide, for example, for coordination of the other components of the mobile computing device 1450, such as control of user interfaces,
applications run by the mobile computing device 1450, and wireless communication by the mobile computing device 1450.
[0108] The processor 1452 may communicate with a user through a control interface
1458 and a display interface 1456 coupled to the display 1454. The display 1454 may be, for example, a TFT (Thin-Film-Transistor Liquid Crystal Display) display or an OLED (Organic Light Emitting Diode) display, or other appropriate display technology. The display interface 1456 may comprise appropriate circuitry for driving the display 1454 to present graphical and other information to a user. The control interface 1458 may receive commands from a user and convert them for submission to the processor 1452. In addition, an external interface 1462 may provide communication with the processor 1452, so as to enable near area communication of the mobile computing device 1450 with other devices. The external interface 1462 may provide, for example, for wired communication in some implementations, or for wireless communication in other implementations, and multiple interfaces may also be used.
[0109] The memory 1464 stores information within the mobile computing device 1450.
The memory 1464 can be implemented as one or more of a computer-readable medium or media, a volatile memory unit or units, or a non-volatile memory unit or units. An expansion memory
1474 may also be provided and connected to the mobile computing device 1450 through an expansion interface 1472, which may include, for example, a SIMM (Single In Line Memory
Module) card interface. The expansion memory 1474 may provide extra storage space for the mobile computing device 1450, or may also store applications or other information for the mobile computing device 1450. Specifically, the expansion memory 1474 may include instructions to carry out or supplement the processes described above, and may include secure information also. Thus, for example, the expansion memory 1474 may be provided as a security module for the mobile computing device 1450, and may be programmed with instructions that permit secure use of the mobile computing device 1450. In addition, secure applications may be provided via the SIMM cards, along with additional information, such as placing identifying information on the SIMM card in a non-hackable manner.
[0110] The memory may include, for example, flash memory and/or NVRAM memory
(non-volatile random access memory), as discussed below. In some implementations, instructions are stored in an information carrier and, when executed by one or more processing devices (for example, processor 1452), perform one or more methods, such as those described above. The instructions can also be stored by one or more storage devices, such as one or more computer- or machine -readable mediums (for example, the memory 1464, the expansion memory 1474, or memory on the processor 1452). In some implementations, the instructions can be received in a propagated signal, for example, over the transceiver 1468 or the external interface 1462.
[0111] The mobile computing device 1450 may communicate wirelessly through the communication interface 1466, which may include digital signal processing circuitry where necessary. The communication interface 1466 may provide for communications under various modes or protocols, such as GSM voice calls (Global System for Mobile communications), SMS (Short Message Service), EMS (Enhanced Messaging Service), or MMS messaging (Multimedia Messaging Service), CDMA (code division multiple access), TDMA (time division multiple access), PDC (Personal Digital Cellular), WCDMA (Wideband Code Division Multiple Access), CDMA2000, or GPRS (General Packet Radio Service), among others. Such communication may occur, for example, through the transceiver 1468 using a radio-frequency. In addition, short-range communication may occur, such as using a Bluetooth®, Wi-Fi™, or other such transceiver (not shown). In addition, a GPS (Global Positioning System) receiver module 1470 may provide additional navigation- and location-related wireless data to the mobile computing
device 1450, which may be used as appropriate by applications running on the mobile computing device 1450.
[0112] The mobile computing device 1450 may also communicate audibly using an audio codec 1460, which may receive spoken information from a user and convert it to usable digital information. The audio codec 1460 may likewise generate audible sound for a user, such as through a speaker, e.g., in a handset of the mobile computing device 1450. Such sound may include sound from voice telephone calls, may include recorded sound (e.g., voice messages, music files, etc.) and may also include sound generated by applications operating on the mobile computing device 1450.
[0113] The mobile computing device 1450 may be implemented in a number of different forms, as shown in the figure. For example, it may be implemented as a cellular telephone 1480. It may also be implemented as part of a smart-phone 1482, personal digital assistant, or other similar mobile device.
[0114] Various implementations of the systems and techniques described here can be realized in digital electronic circuitry, integrated circuitry, specially designed ASICs (application specific integrated circuits), computer hardware, firmware, software, and/or combinations thereof. These various implementations can include implementation in one or more computer programs that are executable and/or interpretable on a programmable system including at least one programmable processor, which may be special or general purpose, coupled to receive data and instructions from, and to transmit data and instructions to, a storage system, at least one input device, and at least one output device.
[0115] These computer programs (also known as programs, software, software applications or code) include machine instructions for a programmable processor, and can be implemented in a high-level procedural and/or object-oriented programming language, and/or in assembly/machine language. As used herein, the terms machine-readable medium and computer-readable medium refer to any computer program product, apparatus and/or device (e.g., magnetic discs, optical disks, memory, Programmable Logic Devices (PLDs)) used to provide machine instructions and/or data to a programmable processor, including a machine- readable medium that receives machine instructions as a machine -readable signal. The term
machine-readable signal refers to any signal used to provide machine instructions and/or data to a programmable processor.
[0116] To provide for interaction with a user, the systems and techniques described here can be implemented on a computer having a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user and a keyboard and a pointing device (e.g., a mouse or a trackball) by which the user can provide input to the computer. Other kinds of devices can be used to provide for interaction with a user as well; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form, including acoustic, speech, or tactile input.
[0117] The systems and techniques described here can be implemented in a computing system that includes a back end component (e.g., as a data server), or that includes a middleware component (e.g., an application server), or that includes a front end component (e.g., a client computer having a graphical user interface or a Web browser through which a user can interact with an implementation of the systems and techniques described here), or any combination of such back end, middleware, or front end components. The components of the system can be interconnected by any form or medium of digital data communication (e.g., a communication network). Examples of communication networks include a local area network (LAN), a wide area network (WAN), and the Internet.
[0118] The computing system can include clients and servers. A client and server are generally remote from each other and typically interact through a communication network. The relationship of client and server arises by virtue of computer programs running on the respective computers and having a client-server relationship to each other.
[0119] In view of the structure, functions and apparatus of the systems and methods described here, in some implementations, systems and methods for characterizing potentially cancerous and/or pre-cancerous tissue are provided. Having described certain implementations of methods and apparatus for provided methods, it will now become apparent to one of skill in the art that other implementations incorporating the concepts of the disclosure may be used.
Therefore, the disclosure should not be limited to certain implementations, but rather should be limited only by the spirit and scope of the following claims.
[0120] Throughout the description, where apparatus and systems are described as having, including, or comprising specific components, or where processes and methods are described as having, including, or comprising specific steps, it is contemplated that, additionally, there are apparatus, and systems of the disclosed technology that consist essentially of, or consist of, the recited components, and that there are processes and methods according to the disclosed technology that consist essentially of, or consist of, the recited processing steps.
[0121] It should be understood that the order of steps or order for performing certain action is immaterial so long as the disclosed technology remains operable. Moreover, two or more steps or actions may be conducted simultaneously.
EXAMPLES
Example 1 - Characterization of Breast Tumor Organization
[0122] This Example shows that the 2D WTMM method has the ability to characterize breast tumors and their microenvironment. Data was obtained from the Digital Database for Screening Mammography (DDSM), described above. Cases having a suspicious region, i.e. containing a benign or malignant tumor, were considered. From these cases, only cases having one set of microcalcifications were kept for further investigation and masses were disregarded for this analysis. Using the methods described above, we were able to fully characterize the lesion by the fractal dimension of the breast tumor and the roughness of the microenvironment given by the Hurst exponent.
[0123] In this Example, a total of 128 images, 78 of which are benign and 50 of which are malignant, were analyzed. Since each mammographic case contains two images, one corresponding to the mediolateral oblique (MLO) view, and the other to the cranio-caudal (CC) view with their projections shown in FIG. 1, we combined the information to reconstruct an estimate of the 3D structure of the tumor embedded into the breast tissue. The information was combined by creating a scatter plot of the minimum fractal dimension versus the maximum fractal dimension as shown in FIG. 2.
[0124] The Mclust function in R was applied to the two data sets, the benign and malignant scatter plots. The results from the function shown in FIG. 3, suggests that benign tumors are constructed from three subpopulations and malignant tumors are constructed from two subpopulations. Here, it is important to note the centers of the components identified by the function. Benign breast lesions' three components are centered at (.84,1.04), (.85,1.78) and (1.78,1.88), while malignant breast lesions' two components are centered at (1.07,1.62) and (1.49,1.52).
[0125] The center of the benign categories demonstrate how these tumors are noninvasive since the populations fractal dimensions lay in Euclidean space and are represented by lines, sheets and spheres. However, both malignant categories are shown by invasive fractal dimensions and represented by fractal structures in breast tissue.
Example 2 - Classification of Tumor Prone Tissue
[0126] This example illustrates how provided methods may be used to better understand the mechanisms that drive the differing organization of benign and malignant tumors, and specifically the microenvironment of a radiologist-identified suspicious region.
[0127] To gain a deeper understanding of the behavior and development of breast lesions, neighboring tissue was classified according to the two dimensional WTMM methods described above to survey the microenvironment of the tumors {also see Kestener et al, Wavelet-based multifractal formalism to assist in diagnosis in digitized mammograms, 2001, Image Anal Stereol, 20: 169-174 for an application of this methodology). Briefly, the 2D WTMM method was used to analyze mammary parenchyma. This analysis was used to generate a Hurst exponent (H) as a measure of the roughness of a sample. The Hurst score was used to classify breast tissue as fatty (H = 0.25 ± 0.05), dense (H = 0.65 ± 0.05).
[0128] To categorize the neighboring tissue, each image of a breast lesion was segmented into 9 subimages, corresponding to the same size and shape of the radiologist encircled tumor. Once the images were properly segmented, only the central part of the neighboring images were analyzed to counter any edge effects that may disrupt the analysis. The rough surface images were analyzed via the above methodology.
[0129] As shown in FIG. 4, the results allowed for the characterization of abnormal tissue corresponding to a benign tumor and malignant tumor. In FIG. 4, blue represents H = (0- 0.45), or fatty tissue, red represents H = (0.55 - 1), or dense tissue, and yellow represents H = (0.45; 0.55). It is important to note that in the cancer sample shown in FIG. 4 there are several yellow regions identified to have H = (0.45 - 0.55), where in benign and healthy samples, no such tissue classification was found to exist.
[0130] In order to statistically determine if the microenvironments of benign and malignant breast lesions are significantly different, we constructed a t-test to obtain a p-value. Our hypothesis is that the microenvironment of malignant breast tumors are categorized as H = 0.5, while benign tumors have a microenvironment categorized as H = 0.3 or H = 0.65, corresponding to healthy fatty or dense tissue, respectively. We obtained one value per case by averaging the calculating | H -0.5 | values of the neighboring tissue for both the CC and MLO views. From here, the H values were averaged for each case. By computing ΔΗ =<| H -0.5 |>, we can expect benign tissue to have a higher value characterizing the microenvironment compared to malignant tissue.
[0131] The boxplot in FIG. 5 shows the distribution of ΔΗ for both benign and malignant microenvironments. Here, the max and min values are shown on the plot as 0.33 and 0.07 for benign and 0.26 and 0.08 for cancer respectively. As shown in the boxplot, three outliers were identified for benign and one for cancer, meaning these values for ΔΗ are 1.5 times more than the maximum ΔΗ value.
[0132] Since the boxplot suggests there may be some non-normality of the data, a test for normality was performed in order to decide if a transformation of the data would be necessary before proceeding to the t-test. The Shapiro-Wilk test, first published by S. Shapiro and M. Wilk in 1965, tests the null hypothesis that a sample comes from a normally distributed population. The test statistic is given by:
where x is the sample mean, x(i) is the i order statistic, and
where m = (mi; im; .. . ; ηΐη)τ and mi is expected value of the order statistic and V is the covariance matrix.
[0133] The result from performing the Shapiro-Wilk Test on our data shows that the data was not normally distributed since p = 0.001 for benign and p = 0.033 for cancer. A log transformation was conducted and the results in Table 5 demonstrates normal data.
Table 5 - Results From Shapiro-Wilk Test
[0134] Once the data has been transformed in order for the data sets to follow a normal distribution, we checked to make sure the two data sets have equal variances by using the F-ratio
test. Here, the null hypothesis is that two normal populations have the same variance. The test statistic is:
F = (5.2)
s?. where X and Y are independent and normal, and S2x and S2y are the sample variances. Thus, we would reject Ho if F is too large or too small. The results from the F-ratio test in Table 6 tells us to reject H0.
Table 6 - Results from F-Ratio Test
[0135] Once we have checked the assumptions of the parametric test, we can perform the necessary t-test. We can define μβ as ΔΗ for the benign microenvironment and μC as ΔΗ for the cancer microenvironment. We will be testing the null hypothesis that μC > μΒ, making the alternative hypothesis μC < μΒ. The result, shown in Table 7 , was a p-value of 0.04. Thus, we can reject Ho and conclude there is a significant difference in the log(AH) values characterizing the microenvironments of benign and malignant breast lesions.
Table 7 - T-Test Results of Cancer and Benign ΔΗ
[0136] It is key to point out that H = 0.5 is not found on healthy mammograms in previous studies and less often on the benign sample studied here. Another key piece is that H = 0.5 corresponds to no correlation and a loss of structure. This suggests that the tissue composing the environment of malignant tumors has lost its structure and may even precede the onset or progression of the microcalcifications (MC). Thought to be an indicator of the early development of breast cancer. Thus, this Example shows that calculation of a Hurst exponent according to provided methods provides a powerful tool in the early detection and diagnosis of breast cancer.
Example 3 - Alternative Analysis of Prone Tissue
[0137] As indicated above in Example 2, obtaining one average ΔΗ value per subject by averaging the calculated | H -0.5 | values of neighboring tissue for both the MLO and CC views provided valuable insight that the loss of microenvironmental structure is indicative of the development of tumors, with a higher level of disorder being found in more serious (malignant) tumors.
[0138] In this example, rather than calculating a single average ΔΗ value per subject
(consisting of up to 16 neighbors per image, 8 in each of the MLO and CC views), the distribution of ΔΗ = | H -0.5 | for each of the neighboring regions was taken individually. This methodology provides a significantly larger sample size than the analysis in Example 2 because the sample size becomes the number of subjects x the number of neighboring images around each lesion, for a total of (up to) 16 ΔΗ = | H-0.5 | values per subject. When the data in Example 2 is analyzed using this method, the p-value between malignant and benign tissues is an even more significant p= 0.001.
[0139] The analysis in this Example is likely to be of significant utility in various embodiments because it is better able to account for variations within a subject and may be able to better define cancerous or pre-cancerous regions within a target tissue.
Example 4 - Characterization of Tissue Microenvironment
[0140] In this Example, data from analyzed breast tissue was characterized on a per neighboring region basis, as opposed to a per patient basis. For this Example, a total of 1,131 images of tissue classified as malignant, 1,131 images of tissue classified as benign, and 64 images of tissue classified as normal were analyzed. Specifically, the number of regions in analyzed microenvironments characterized as fatty, dense or disrupted were calculated from the samples examined in Example 1. Table 8 shows the results of this analysis:
Table 8 - Classification of Tissue Microenvironment
[0141] As shown in Table 8, the percentage of disrupted regions was highest in malignant tissue, followed by benign tissue, with normal tissue having the lowest levels of disruption across the tissue. Exemplary photographs of the differences in observed disruption may be found in FIG. 6. It is important to note that the percentage of disrupted regions for normal breast tissue is not 0%>. Without wishing to be held to a particular theory, it is possible these disrupted regions may be representative of a transitional state of the tissue that may be prone to the eventual development of a tumor. Longitudinal studies will help to determine the nature of disrupted regions in breast tissue classified as normal by current diagnostic
methodologies.
[0142] In order to more accurately calculate the proportion of disrupted regions in malignant, benign, and normal tissue (particularly given the low number of normal tissue images analyzed above), a second study of the same design as that described above was performed. In the second study a total of 1,081 malignant images, 1,122 benign images, and 740 images
classified as normal (more than ten times that number from the first study) were included. The results are shown below in Table 9:
Table 9 - Classification of Tissue Microenvironment
p = 9.4 x 10" as compared to benign tissue; p = 2.2 x 10" as compared to malignant tissue
** p = 0.0012 as compared to malignant tissue
[0143] As shown in Table 9, using a larger sample size of tissue classified as normal, the proportion of tissue exhibiting disrupted surrounding regions drops from 10.9% in Table 8 to 5.9% in Table 9. The lower observed proportion of images classified as normal exhibiting disrupted tissue further supports the use of roughness and disorder in surrounding tissue as an indicator of cancer or pre-cancer in a tissue, and possibly malignancy. Without wishing to be held to a particular theory, it is possible that tissue otherwise classified as normal that exhibits some degree of disrupted tissue in the surrounding microenvironment may be pre-cancerous or have an increased likelihood of developing cancer in the future as compared to tissue that does not exhibit disruption in the surrounding microenvironment.
Example 5 - Detection of Microcalcification (MC) Clusters and Calculation of their Fractal Dimension
Methods - Samples
[0144] In this Example, images that were analyzed were obtained from the Digital
Database for Screening Mammography (DDSM) at the University of South Florida. The databank contains over 2,600 studies made up of normal, benign, benign without call back and malignant mammograms all categorized by an expert radiologist. Each study has two images of each breast, consisting of a mediolateral-oblique (MLO) view and cranio-caudal (CC) view with any suspicious region circled by a radiologist. The suspicious region could contain a mass and/or microcalcifications (MC), but only the cases that were classified as having exactly one
tumor composed of only MC in the benign and malignant categories were looked at in this particular study.
[0145] In addition to only considering tumors consisting of MC, any mammographic images that contained artifacts inside the radiologist's encircled region were discarded due to the impact it has on the analysis. These artifacts could include scratches, hair, deodorant, patient movement, scanner artifacts (rollers slipped), pacemaker, breast implants, skin markers (for scars, moles, and nipples, as well as marked lumps of breast pain), metallic foreign bodies, and fingerprints. Some (but not all) of these effects were recorded under notes in the DDSM website. In this Example, a total of 59 cases were considered, corresponding to 118 images of size greater than 2562 pixels, 34 of which are benign (68 images) and 25 of which are malignant (50 images).
Methods - The 2D WTMM Method
[0146] In this Example, the two dimensional (2D) WTMM method is used to
characterize images of breast tissue, as described above. Most of the existing CAD methods, whether specifically designed for two dimensional (2D) mammograms or more recently, for three dimensional (3D) breast tomosynthesis, have been elaborated on the prerequisite that the background roughness fluctuations of normal breast texture are statistically homogeneous and uncorrected, which precludes their ability to adequately characterize background tissue. The majority of the fractal methods used to examine and classify mammographic breast lesions rely on the estimate of the Hurst exponent (or its various mathematical equivalents), which globally characterizes the self-similar properties of the landscape in question. However, the 2D WTMM method takes in account that the function defining the image may be multifractal, therefore requiring the use of the Holder exponent (M = Kah) to characterize the local regularity at a particular point.
[0147] The 2D WTMM method requires us to define a smoothing function, 0(x, y), in two dimensions that is a well-localized isotropic function around the origin. In this Example, we used the Gaussian function and define the wavelets as:
^ fe y) = ^ and ^2 fc ) = ^ .
The wavelet transform with respect to ψι and ψ2 is τψ[/Ί( ) )
Τψ [f] (b, a) (Κψ [/] ¾ α), ^[/] ¾ α))
[0148] The wavelet transform modulus maxima are defined as the locations b where M ■,ψ
[/] (b,a) is locally maximum in the direction of the argument Αψ [f] (b,a), at a given scale a > 0. The WTMM lie on connected chains and are thus called maxima chains (FIG. 7A-7F). One can then find the maxima along these WTMM chains. The WTMM maxima, or WTMMM are defined as the points along the maxima chains where the Μψ [/] (b,a) is locally maximum. The
WTMMM are linked through scales to form the space-scale skeleton (Fig. 7G). Hence, one can identify the singularities of the function as the loci x where the WTMMM lines of the WT skeleton (Fig. 7G) point to in the limit a→ 0+. Along these space-scale vertical lines the WTMMM behave as a power-law ~ah^ (M =kah) from which one can extract the local Holder exponent h(x). The multifractal formalism amounts to characterize the relative contributions of each Holder exponent value via the estimate of the so-called D(h) singularity spectrum defined as the fractal dimension of the set of points x where h(x) = h. To compute D(h) we therefore use wavelets to partition the surface by defining the partition function directly from the WTMMM in the skeleton:
where L(a) is the set of all vertical space-scale lines in the skeleton that exist at the given scale a > 0 and which contain maxima at any scale a' < a and q ε R. One can then define the exponent τ (q) from the power-law behavior of the partition function:
+
And the D(h) singularity spectrum of/ can be determined from the Legendre transform of the partition function scaling exponent
D (/i) = mm l q<7 {qh— r(q))
In practice, to avoid instabilities in the estimation of the singularity spectrum D(h) through the Legendre transform, we used h and D(h) as mean quantities defined in a canonical ensemble, i.e. with respect to their Boltzmann weights computed from the WTMMM:
Then one computes the expectation values: h(q, a) =∑£e£(a) ln |{x>t,0ei,fl,≤a^ [ ] (x. a') | W4,[f}(q.£, a),
D (tjf, a) =∑.£εχ(α) Wip [f] (q, L, a) In (w,p [f] (q, L, a) from which one derives i airr(Uqf).,
= imi
ddqq a- >Q+ h(q, a)/ \n -
D (q)— lima→0+ D (qf a) / In a and thus the singularity spectrum D(h) as a curve parameterized by q.
[0149] Homogeneous monofractal functions with singularities of unique Holder exponent
H are characterized by a linear τ (q) curve of slope H. A nonlinear τ (q) is the signature of nonhomogeneous multifractal functions, meaning that the Holder exponent is a fluctuating quantity that depends on x. Then the corresponding singularity spectrum has a characteristic single-humped shape. Note that for both mono- and multifractal functions
D (cj = 0) = - r(q = 0) = DF
Where DF (noted simply D throughout the text) is the fractal dimension of the support of singularities of f.
Methods - Statistical Tests
[0150] The Wilcoxon rank-sum test is a non-parametric statistical hypothesis test that is used as an alternative to Student's t-test when the population cannot be assumed to be normally distributed. It was used here to calculate the /^-values comparing the CC and MLO fractal dimensions and benign and malignant cases images since the benign data followed a bimodal distribution (FIG. 10). The calculations were done using the Wilcox test in R.
Methods - Bayesian Statistics
[0151] Bayes theorem states that
ρ(θ \ 0) = p(G \ 9) * ρ(θ) / p(G)
posterior likelihood prior evidence
where the prior, ρ(θ) , represents the strength of our belief in malignant lesions p(M) or benign lesions p(B) out of those that have been diagnosed by a radiologist. The posterior, p(9IG), represents the strength of our belief, having accounted for the geometrical evidence, G, where G represents the position of the lesion in the fractal dimension plot, either fractal (F) or Euclidean (E). The quotient of the likelihood over the evidence, p(9IG)/p(G), represents the support the evidence, G, provides for Θ. Since the prior reflects uncertainty in the parameter value θ, ρ(θ) was based on a Beta distribution with a specified mean and standard deviation. The Beta distribution is as follows:
Κθ,- α, β) = Γ(α+β) Θα-Ηΐ - θ →
[0152] The Probability Model: To estimate the mean of the Beta distribution for malignant cases, p(M), the prevalence of mammograms with a BI-RADS assessment score of 3, 4, and 5 were used as determined by the radiologists diagnostics multiplied by the historical probability of mammograms receiving those assessment scores result in malignant MC clusters respectively. One out of the 59 cases considered in this study received an assessment score of 3,
47 out of 59 received a 4, and 11 out of 59 received a 5. Based on historical data, the probability of malignant lesions given an assessment score of 3 is 2%, an assessment score of 4 is 26.5% (taken as the midpoint of the reported interval of [23%>-30%>], and an assessment score of 5 is 95%. Therefore the Beta distribution for p(M) was chosen with a mean of
1 * 0.02 + 47 * 0.265 + 11 * 0.95
= 0.3885
1 + 47 + 11
and the Beta distribution for p(B) with a mean of 1-0.3885 = 0.6115. However, since there isn't much certainty regarding these values, the Beta distribution for p(M) and p(B) were assigned a relatively large standard deviation of 0.25. This resulted in a Beta distribution for p(M) with parameters (α, β) = (1.09241,1.71223) and p(B) with parameters (α, β) = (1.71223,1.09241). The likelihood p(FlM), where F represents breast lesions characterized as being in the fractal 465 zone (and E likewise represents those in the Euclidean zone), is based on the 23 of the 25 malignant cases that were in the fractal zone; this likelihood is the probability that the data could be generated with parameter values Θ. Similarly, the likelihood p(ElB) is based on the 30 out of 34 benign cases that were in the Euclidean zone. To arrive at the posterior distributions p(MIF) and p(BIE), the R routine "BernBeta.R" was used.
[0153] The resulting posterior distribution for p(MIF) was a Beta distribution with parameters (a = 24.0924, β = 3.71223). Based on this distribution, the resulting 95% highest density interval was (0.742, 0.975. The posterior distribution for P(BIE) was a Beta distribution with parameters (a = 32.7122, β = 5.09241). Based on this distribution, the resulting 95% highest density interval was (0.757, 0.962). The highest density interval spans 95% of the posterior distribution such that every point inside the interval is deemed more credible. In other words, given the prior and the likelihood, observing the parameter value for the percentage of breast lesions characterized in the fractal zone that are malignant, there is a 95% probability that this parameter is between 0.742 and 0.975. Similarly, for the percentage of breast lesions characterized in the Euclidean zone that are benign, there is a 95% probability that this parameter is between 0.757 and 0.962.
[0154] As described in detail in the 2D WTMM Method section and in FIG. 7, the wavelet transform (WT) acts as a mathematical microscope to characterize spatial image
information over a continuous range of size scales. It is the gradient vector of the image smoothed by dilated versions of a Gaussian filter. At each size scale, the wavelet transform modulus maxima (WTMM) are defined by the positions where the modulus of the WT is locally maximal. These WTMM are automatically organized as maxima chains at the considered scale. Along each of these chains, further local maxima are found, i.e., the WTMM maxima
(WTMMM). This process is repeated for all size scales and the WTMMM from each scale are then linked to form the WT skeleton. As shown in FIGS. 8 and 9, the ability to consider (vertical) space-scale WTMMM lines in the WT skeleton individually is key, since it allows us to objectively discriminate between lines pointing to the tissue background from those pointing to the microcalcifications by considering how the WT modulus varies as a function of the scale parameter along each space-scale line. In FIGs. 8D and 9D, each space-scale line obtained from the WT skeleton is represented by plotting the evolution of the WT modulus, M (see the 2D WTMM section above), as a function of the scale parameter, a, in a log-log plot. This relationship between M and a is characterized by a power-law behavior via the equation
M = Kah where T is a pre-factor and h is the Holder exponent quantifying the strength of the singularity to which the space-scale line is pointing to. In log-log plots shown in FIGs. 8D and 9D, the slope of the curves therefore corresponds to h. By considering two types of information characterizing the behavior of a space-scale line, namely the strength of the modulus at the smallest scale, which is given by the log of the pre-factor, log(AT), as well as the slope (in a logarithmic representation) of the modulus variation across scales, h, a classification procedure is setup which results in two sets of space-scale lines that clearly segregate MC from background tissue. Isolated MC can be seen as Dirac-like singularities through the optics of the WT, for which h is theoretically known to be -1. However, while clustered MC may not appear as isolated
200 Dirac-like singularities, the edge that they form is still easily detectable through the space- scale lines, with a value of h ~ 0 (discontinuity). This means that for both isolated and clustered MC, we can expect the space-scale lines to behave as M = kah with h<0, which contrasts from the healthy background tissue, for which A-l/3 for fatty breast tissue and h~2/3 for dense breast tissue. However, relying only on h may not be sufficient, which is why the strength of the
WT modulus at the lowest scale, which quantifies the contrast between MC and background, is also needed. For the sample case presented in FIG. 8, the plot in Fig. 8D shows that neither log( T) nor h, taken individually, would have been sufficient to segregate between MC and background. However, for the sample case presented in FIG. 9, the plot in Fig. 9D shows that log( T) alone was sufficient. A more detailed discussion of both cases follows.
[0155] In FIG. 8 the background breast tissue is dense, which makes the contrast between background and MC weak (i.e. causing a low value for the WT modulus of red curves at the smallest scale in FIG. 8D). However, the roughness fluctuations of dense breast tissue are characterized by a relatively high smoothness level, which translates to blue curves with a large slope (i.e., a high h value, -2/3) for scales 10 < a < 200 pixels as compared to the red curves with negative slopes for scales a > 10 pixels that correspond to WTMMM lines that point to MC at small scale (a→ 0+) (FIG. 8D). In FIG. 9, the background breast tissue is fatty, which is characterized by a higher roughness level (i.e. a lower h value -1/3, although still positive), that reduces the discriminatory power of h. However, for MC embedded into fatty tissue, the contrast is high, which translates to a high value of log( T). Therefore, applying a threshold on both parameters, h and log(K), is key to segregating MC from their background regardless of the density (fatty or dense) of the composition of the breast tissue.
[0156] Once this segregation is done, the so-called singularity spectrum can then be calculated separately for each subset, which then allows us to consider the fractal dimension D of the lesion, characterizing its architecture.
[0157] We restricted the analysis of DDSM cases (see Methods sections above) to images having a radiologist encircled region that was larger than 2562 pixels for both views (CC and MLO) and also to make sure that the distribution of patient ages was comparable (i.e. 56.7
+/- 1 1.4 years old for the benign cases and 65.5 +/- 12.4 years old for the malignant cases). This resulted in an analyzed sample with a total of 59 cases (1 18 images), 34 of which are benign (68 images) and 25 of which are malignant (50 images). The histograms of fractal dimension values obtained are presented in FIG. 10. Note that blending the CC and MLO fractal dimensions together in these distributions would not guarantee an unbiased statistical analysis, which is why the fractal dimension values for the CC and MLO results were analyzed independently. FIG.10 demonstrates that benign MC clusters have a strong preference for Euclidean dimensions that
are either close to D=\ or to D=2 and that there is a very clear zone of avoidance in the fractal range, i.e., for 1<D<2, with an actual gap in the benign histograms for the bin centered at D=\ .5 for both views. For the malignant cases, it is the opposite: Euclidean zones are avoided and the data are very clearly centered in the fractal range for both views, with the peak of the
histograms at D=\ .5. A statistical comparison between benign vs. malignant MC clusters was performed using the Wilcoxon rank-sum test, which yielded /^-values of 0.009 for CC
comparisons and 0.014 for MLO comparisons for the benign vs. malignant fractal dimension distributions. These are statistically significant differences.
Example 6 - The CC-MLO Fractal Dimension Plot and Fractal Zone
[0158] In this Example, a fractal dimension plot between the MLO and CC views of a breast is calculated. Unless otherwise specified, all samples, sample handling, and analysis was as described above in Example 5. As shown in Example 5, the significance of the difference between benign and malignant is quite interesting. However, it is still only based on statistics of populations. The histograms in FIG. 10 show that, when each view is taken independently (CC or MLO), it is still possible, though unlikely, for a malignant lesion to have a Euclidean dimension, and vice-versa, for benign lesions to have a (non-integer) fractal dimension.
However, in order to work towards a potential CAD method that would be able to diagnose breast lesions individually instead of via the population statistics, we combined the information to indirectly infer the 3D structure of the tumor embedded into the breast tissue. This is presented in a novel plot called the "CC-MLO fractal dimension plot" shown in FIG. 11, where red dots represent malignant cases and green dots represent benign cases. The square centered at (1.5, 1.5) represents those cases for which both CC and MLO views have a fractal dimension that is within 12<D<\ .8. Note that only malignant cases are found in this internal square. However, having one of the two views with a score that is close to D=\ .5 may "compensate" for its other view being outside of the [1.2, 1.8] x [1.2, 1.8] square, i.e., as one view approaches D=1.5, the farther from 1.5 the other can be. Further justification is presented below and in FIG. 12. That is how the triangular regions that decay linearly as a function of distance from the internal square were defined. Therefore, the central square, combined with the four triangular regions extending
from it are what we define as the "fractal zone". Of the 59 cases considered in this study, 92% of malignant breast lesions studied (23 out of 25) were in the fractal zone while 88% of the benign lesions were in the Euclidean zones (30 out of 34).
Bayesian Analysis
[0159] The inferences from a Bayesian analysis are richer and more informative than null hypothesis significance testing. In particular, there is no reliance on /^-values. But also,
Bayesian models are designed to be appropriate to the data structure without having to make approximation assumptions typical in null hypothesis significance testing. The results reported above show that the vast majority of malignant breast lesions are fractal, and that the vast majority of benign breast lesions are Euclidean. However, the condition of interest is how breast lesions in the fractal zone can indicate malignancy, and how breast lesions in the Euclidean zone can indicate benignancy.
[0160] Bayesian inference derives the posterior probability as a consequence of two antecedents, a prior probability and a likelihood function derived from a probability model for the data to be observed. In this application, the model is based on historical radiology assessment scores using the BI-RADS system. A detailed description of this probability model as well as the mathematical model behind Bayes analysis is presented in the Methods section. Bayesian inference then becomes a computation of the posterior probability according to Bayes' rule. The interpretable output of this Bayesian analysis is the so-called 95% highest density interval (HDI), which is analogous to the 95% confidence interval in frequentist statistics. The 95% HDI from the resulting posterior distribution indicates that the percentage of breast lesions in the fractal zone that are malignant is between 74.2% and 97.5%. Alternatively, in terms of controlling for false positives, which is a major concern, as discussed in the Introduction, the percentage of breast lesions in the Euclidean zone that are benign is between 75.7% and 96.2%.
Interpretation of the 3D Geometrical Structure
[0161] Even though two different 2D views of a 3D object are insufficient to fully characterize its 3D geometry, it can nonetheless give a robust estimate, as shown in this
Example. The cases outside of the fractal zone can be categorized in two Euclidean subsets: 1) LINES, i.e. those that are approximately in the
area, which are therefore seen
as one-dimensional objects from both views (FIG. 12A, left panel); or 2) SHEETS, i.e. those that are either in the (Αχ=1, -DMLO=2) or (DCc=2, DMLO=1) areas, which are seen as a full two- dimensional object in one view, but as a one-dimensional object from the other view and also those that are in the (Dcc=2, DMLO=2) area, which are seen as a full two-dimensional object from both views (FIG. 12B, center panel). Although simplistic, these case models represent a good estimate of what the 3D Euclidean structure of a benign lesion may look-like.
[0162] For the cases that fall in the fractal zone, those malignant lesions that are in the
[1.2,1.8] x[1.2,1.8] square have a fractal signature that is seen from both views, whereas those that are in the triangular areas would represent fractal clusters that grow onto a 2D plane, i.e. seen as a fractal from one view, but seen either as a line (bottom or left triangular regions) or a plane (top or right triangular regions) from the other view (FIG. 12C, right panel). Interestingly, a diffusion-limited aggregate embedded in 3D space and for which 2 < D < 3 will have a 2D projection with D=2. Since no malignant lesions are found in the (Dcc=2, DMLO=2) area of the CC-MLO fractal dimension plot, we can safely hypothesize that all tumors are essentially limited to a 2-dimensional fractal structure (within the 3-dimensional breast tissue), for which 1 < D < 2. Without wishing to be held to a particular theory, this therefore leads us to conjecture that all breast tumors considered in this study, benign and malignant, fractal or Euclidean, would grow on 2-dimensional manifolds.
[0163] Examples 5 and 6 illustrate embodiments of provided methodologies which offer a way to accurately classify benign and malignant tumors based on their invasiveness as determined by the geometrical structure. By considering the organization of tumors via CAD systems, current mammographic practice may be improved by increasing accuracy, and potentially decreasing recall rates and costs. The inferred 3-dimensional geometry of the breast lesions based on the analysis of the mammographic images using the 2D WTMM methodology allows one to explore the invasiveness of the breast tumors and provide an interpretation of the severity of the lesion. By considering where each case falls on the CC-MLO fractal dimension plot, a score similar to the Breast Imaging-Reporting and Data System (BI-RADS) assessment score may be assigned to each case. Not only does this tool have the potential as a CAD, but it may also provide insight into the underlying mechanisms that drive overall growth and structure at the time of the screening mammogram.
[0164] Without wishing to be held to a particular theory, it is contemplated that, since the structure of the tumors are different, with benign lesions likely being smooth Euclidean objects and malignant lesions being branching objects (and possibly, for both cases, being restricted to growing along 2D manifolds within the 3D breast tissue environment), there may be a link to the cellular mechanisms at the lower levels in the system that drive the organization at the much larger scale of mammograms. Use of provided methods may be helpful in revealing these relationships and improving detection and treatment.
EQUIVALENTS AND SCOPE
[0165] Those skilled in the art will recognize, or be able to ascertain using no more than routine experimentation, many equivalents to the specific embodiments of the invention described herein. The scope of the present invention is not intended to be limited to the above Description, but rather is as set forth in the following claims:
Claims
1. A method comprising
providing a first region of tissue from a subject;
calculating a roughness exponent for the first region of tissue; and
comparing the roughness exponent of the first region of tissue to 0.5, wherein a difference of less than 0.2 between the roughness exponent of the first region of tissue and 0.5 indicates that tissue adjacent to the first region of tissue is cancerous or pre-cancerous.
2. The method of claim 1, further comprising treating the tissue adjacent to the first region of tissue if it is cancerous or pre-cancerous.
3. The method of claim 1, wherein the providing, calculating, and comparing steps are each performed a plurality of times, each on a different region of tissue from a subject.
4. The method of claim 3, wherein the plurality is at least 10.
5. The method of claim 3, wherein the plurality is at least 100.
6. The method of claim 3, wherein the plurality is at least 1,000.
7. The method of claim 1, wherein the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
8. The method of claim 1, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.15.
9. The method of claim 1, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.1.
10. The method of claim 1, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.05.
11. The method of claim 1 , wherein the roughness exponent is calculated from one or more images of the first region of tissue.
12. The method of claim 1, wherein the roughness exponent is calculated using one or more multi-scale analytical methods selected from a wavelet-transform modulus maxima, a wavelet leader, detrended fluctuation, and Fourier analysis.
13. A method comprising
providing a first view of a region of tissue;
providing a second view of the region of tissue;
calculating a first fractal dimension for the first view of the region of tissue; and calculating a second fractal dimension for the second view of the region of tissue;
wherein if the fractal dimension of at least one of the first fractal dimension and the second fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
14. The method of claim 13, further comprising treating the region of tissue if it is cancerous.
15. The method of claim 13, wherein the fractal zone is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
16. The method of claim 15, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8);
first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8);
second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2);
third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and
fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
17. The method of claim 15, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9);
first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9);
second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1);
third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and
fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
18. The method of claim 15, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7);
first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7);
second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3);
third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and
fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
19. The method of claim 13, further comprising
providing a third view of the region of tissue; and
calculating a third fractal dimension for the third view of the region of tissue;
wherein if the fractal dimension of at least one of the first fractal dimension, the second fractal dimension, and the third fractal dimension is in the fractal zone, the region of tissue is considered cancerous.
20. The method of claim 13, wherein the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
21. An apparatus comprising:
a memory for storing a set of instructions; and
a processor for executing the set of instructions, wherein the instructions, when executed, cause the processor to:
provide a representation of a first region of tissue from a subject; calculate a roughness exponent for the first region of tissue; and
compare the roughness exponent of the first region of tissue to 0.5.
22. The apparatus of claim 21, wherein the providing, calculating, and comparing steps are each performed a plurality of times, each on a different region of tissue from a subject.
23. The apparatus of claim 22, wherein the plurality is at least 10.
24. The apparatus of claim 22, wherein the plurality is at least 100.
25. The apparatus of claim 22, wherein the plurality is at least 1,000.
26. The apparatus of claim 21, wherein the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
27. The apparatus of claim 21, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.15.
28. The apparatus of claim 21, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.1.
29. The apparatus of claim 21, wherein the difference between the roughness exponent of the first region of tissue and 0.5 is less than or equal to 0.05.
30. The apparatus of claim 21, wherein the roughness exponent is calculated from one or more images of the first region of tissue.
31. The apparatus of claim 21 , wherein the roughness exponent is calculated using one or more multi-scale analytical methods selected from a wavelet-transform modulus maxima, a wavelet leader, detrended fluctuation, and Fourier analysis.
32. An apparatus comprising:
a memory for storing a set of instructions; and
a processor for executing the set of instructions, wherein the instructions, when executed, cause the processor to:
provide a first view of a region of tissue;
provide a second view of the region of tissue;
calculate a first fractal dimension for the first view of the region of tissue; and calculate a second fractal dimension for the second view of the region of tissue; and determine if one or more of the fractal dimension and the second fractal dimension is in the fractal zone.
33. The apparatus of claim 32, wherein the fractal zone is defined as a polygon consisting of a central square and a first, second, third, and fourth extending triangular region as plotted on a graph of the fractal dimension of the first view of the region of tissue by the fractal dimension of the second view of the region of tissue.
34. The apparatus of claim 33, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.2, 1.2), (1.2, 1.8), (1.8, 1.2)(1.8, 1.8);
first extending triangular region: (0.5, 1.5), (1.2, 1.2), (1.2, 1.8);
second extending triangular region: (1.5, 0.5), (1.2, 1.2), (1.8, 1.2);
third extending triangular region: (1.5, 2.3), (1.2, 1.8), (1.8, 1.8); and
fourth extending triangular region: (2.3, 1.5), (1.8, 1.2), (1.8, 1.8).
35. The apparatus of claim 33, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.1, 1.1), (1.1, 1.9), (1.9, 1.1)(1.9, 1.9);
first extending triangular region: (0.5, 1.5), (1.1, 1.1), (1.1, 1.9);
second extending triangular region: (1.5, 0.5), (1.1, 1.1), (1.9, 1.1);
third extending triangular region: (1.5, 2.3), (1.1, 1.9), (1.9, 1.9); and
fourth extending triangular region: (2.3, 1.5), (1.9, 1.1), (1.9, 1.9).
36. The apparatus of claim 33, wherein the central square and first, second, third and fourth extending triangular regions are defined by the following fractal dimensions, as plotted on a graph of the fractal dimension of the first view by the fractal dimension of the second view: central square: (1.3, 1.3), (1.3, 1.7), (1.7, 1.3)(1.7, 1.7);
first extending triangular region: (0.5, 1.5), (1.3, 1.3), (1.3, 1.7);
second extending triangular region: (1.5, 0.5), (1.3, 1.3), (1.7, 1.3);
third extending triangular region: (1.5, 2.3), (1.3, 1.7), (1.7, 1.7); and
fourth extending triangular region: (2.3, 1.5), (1.7, 1.3), (1.7, 1.7).
37. The apparatus of claim 32, wherein the processor further:
provides a third view of the region of tissue; and
calculates a third fractal dimension for the third view of the region of tissue.
38. The apparatus of claim 32, wherein the tissue is selected from breast tissue, brain tissue, lung tissue, kidney tissue, liver tissue, uterine tissue, dermal tissue, and pancreatic tissue.
Priority Applications (7)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
ES14788804T ES2930503T3 (en) | 2013-04-23 | 2014-04-23 | Improved Tissue Characterization Methods |
DK14788804.4T DK2988659T3 (en) | 2013-04-23 | 2014-04-23 | IMPROVED METHODS FOR TISSUE CHARACTERIZATION |
EP22199368.6A EP4151146A1 (en) | 2013-04-23 | 2014-04-23 | Improved methods of cancer detection |
EP14788804.4A EP2988659B1 (en) | 2013-04-23 | 2014-04-23 | Improved methods of tissue characterization |
PL14788804.4T PL2988659T3 (en) | 2013-04-23 | 2014-04-23 | Improved methods of tissue characterization |
US14/786,366 US10467755B2 (en) | 2013-04-23 | 2014-04-23 | Methods of cancer detection |
US16/577,238 US10769790B2 (en) | 2013-04-23 | 2019-09-20 | Methods of cancer detection |
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201361815209P | 2013-04-23 | 2013-04-23 | |
US61/815,209 | 2013-04-23 |
Related Child Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/786,366 A-371-Of-International US10467755B2 (en) | 2013-04-23 | 2014-04-23 | Methods of cancer detection |
US16/577,238 Division US10769790B2 (en) | 2013-04-23 | 2019-09-20 | Methods of cancer detection |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2014176340A1 true WO2014176340A1 (en) | 2014-10-30 |
Family
ID=51792363
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2014/035153 WO2014176340A1 (en) | 2013-04-23 | 2014-04-23 | Improved methods of cancer detection |
Country Status (6)
Country | Link |
---|---|
US (2) | US10467755B2 (en) |
EP (2) | EP2988659B1 (en) |
DK (1) | DK2988659T3 (en) |
ES (1) | ES2930503T3 (en) |
PL (1) | PL2988659T3 (en) |
WO (1) | WO2014176340A1 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112004464A (en) * | 2018-04-24 | 2020-11-27 | 皇家飞利浦有限公司 | Tumor margin assessment |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104303184B (en) * | 2012-03-21 | 2018-05-15 | 皇家飞利浦有限公司 | Integrate the Clinical workstations of imaging of medical and biopsy data and use its method |
EP2988659B1 (en) | 2013-04-23 | 2022-10-26 | University of Maine System Board of Trustees | Improved methods of tissue characterization |
CN107155300A (en) * | 2014-06-21 | 2017-09-12 | 迈克尔·坦蒂洛 | The method and apparatus selected for breast implant surgery and implant |
KR20170059629A (en) * | 2015-11-23 | 2017-05-31 | 성신여자대학교 산학협력단 | Apparatus for detecting hemorrhage in brain magnetic resonance imaging, method thereof and computer recordable medium storing the method |
JP6774813B2 (en) * | 2016-08-12 | 2020-10-28 | 日本電子株式会社 | Image processing equipment, image processing methods, and analyzers |
CN110364250B (en) * | 2019-04-26 | 2021-12-14 | 深圳大学 | Automatic labeling method, system and storage medium for mammary gland molybdenum target image |
CN110533636B (en) * | 2019-07-31 | 2023-08-04 | 平安科技(深圳)有限公司 | Image analysis device |
CN118229684B (en) * | 2024-05-23 | 2024-08-13 | 深圳大学总医院 | Method and system for identifying adrenal pheochromocytoma |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20040259144A1 (en) * | 2003-05-27 | 2004-12-23 | Nagabhushana Prabhu | Method for diagnosing and treating breast cancer |
US7058210B2 (en) * | 2001-11-20 | 2006-06-06 | General Electric Company | Method and system for lung disease detection |
US7761240B2 (en) * | 2004-08-11 | 2010-07-20 | Aureon Laboratories, Inc. | Systems and methods for automated diagnosis and grading of tissue images |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6081612A (en) | 1997-02-28 | 2000-06-27 | Electro Optical Sciences Inc. | Systems and methods for the multispectral imaging and characterization of skin tissue |
US20040151358A1 (en) * | 2003-01-31 | 2004-08-05 | Akiko Yanagita | Medical image processing system and method for processing medical image |
JP2010504129A (en) * | 2006-09-22 | 2010-02-12 | コーニンクレッカ フィリップス エレクトロニクス エヌ ヴィ | Advanced computer-aided diagnosis of pulmonary nodules |
ES2795036T3 (en) * | 2008-07-25 | 2020-11-20 | Fund D Anna Sommer Champalimaud E Dr Carlos Montez Champalimaud | Systems and methods for treating, diagnosing, and predicting the onset of a medical condition |
WO2011008262A2 (en) * | 2009-07-13 | 2011-01-20 | H. Lee Moffitt Cancer Center & Research Institute | Methods and apparatus for diagnosis and/or prognosis of cancer |
CN101727537A (en) * | 2009-11-16 | 2010-06-09 | 杭州电子科技大学 | Computer determining method of mammary gland CR image based on double visual angle information fusion |
US8923595B2 (en) * | 2010-01-26 | 2014-12-30 | Clarkson University | Method of identification of cancerous and normal cells |
US20130066199A1 (en) * | 2010-04-15 | 2013-03-14 | Cedars-Sinai Medical Center | Tumor margin detection method based on nuclear morphometry and tissue topology |
DE102010024732A1 (en) * | 2010-06-23 | 2011-12-29 | Siemens Aktiengesellschaft | Method and device for detecting tumor-bearing tissue in the gastrointestinal tract with the aid of an endocapsule |
WO2013049153A2 (en) * | 2011-09-27 | 2013-04-04 | Board Of Regents, University Of Texas System | Systems and methods for automated screening and prognosis of cancer from whole-slide biopsy images |
EP2988659B1 (en) | 2013-04-23 | 2022-10-26 | University of Maine System Board of Trustees | Improved methods of tissue characterization |
-
2014
- 2014-04-23 EP EP14788804.4A patent/EP2988659B1/en active Active
- 2014-04-23 EP EP22199368.6A patent/EP4151146A1/en active Pending
- 2014-04-23 DK DK14788804.4T patent/DK2988659T3/en active
- 2014-04-23 ES ES14788804T patent/ES2930503T3/en active Active
- 2014-04-23 WO PCT/US2014/035153 patent/WO2014176340A1/en active Application Filing
- 2014-04-23 US US14/786,366 patent/US10467755B2/en active Active
- 2014-04-23 PL PL14788804.4T patent/PL2988659T3/en unknown
-
2019
- 2019-09-20 US US16/577,238 patent/US10769790B2/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7058210B2 (en) * | 2001-11-20 | 2006-06-06 | General Electric Company | Method and system for lung disease detection |
US20040259144A1 (en) * | 2003-05-27 | 2004-12-23 | Nagabhushana Prabhu | Method for diagnosing and treating breast cancer |
US7761240B2 (en) * | 2004-08-11 | 2010-07-20 | Aureon Laboratories, Inc. | Systems and methods for automated diagnosis and grading of tissue images |
Non-Patent Citations (2)
Title |
---|
KESTENER ET AL.: "Wavelet-based multifractral formalism to assist in diagnosis in digitized mammograms", IMAGE ANAL STEROL, vol. 20, 24 August 2001 (2001-08-24), pages 169 - 179, XP055290510 * |
See also references of EP2988659A4 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112004464A (en) * | 2018-04-24 | 2020-11-27 | 皇家飞利浦有限公司 | Tumor margin assessment |
Also Published As
Publication number | Publication date |
---|---|
US20200082532A1 (en) | 2020-03-12 |
DK2988659T3 (en) | 2022-11-21 |
PL2988659T3 (en) | 2023-01-02 |
US20160055636A1 (en) | 2016-02-25 |
EP4151146A1 (en) | 2023-03-22 |
EP2988659A4 (en) | 2017-01-04 |
EP2988659A1 (en) | 2016-03-02 |
US10769790B2 (en) | 2020-09-08 |
EP2988659B1 (en) | 2022-10-26 |
US10467755B2 (en) | 2019-11-05 |
ES2930503T3 (en) | 2022-12-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US10769790B2 (en) | Methods of cancer detection | |
Singh et al. | An information fusion based method for liver classification using texture analysis of ultrasound images | |
Blanchet et al. | Discrimination between metastasis and glioblastoma multiforme based on morphometric analysis of MR images | |
US9710908B2 (en) | Method and system for assessing fibrosis in a tissue | |
Batchelder et al. | Wavelet-based 3D reconstruction of microcalcification clusters from two mammographic views: new evidence that fractal tumors are malignant and Euclidean tumors are benign | |
Huang et al. | Analysis of tumor vascularity using three-dimensional power Doppler ultrasound images | |
Lee et al. | Associating spatial diversity features of radiologically defined tumor habitats with epidermal growth factor receptor driver status and 12-month survival in glioblastoma: methods and preliminary investigation | |
Xv et al. | Machine learning-based CT radiomics approach for predicting WHO/ISUP nuclear grade of clear cell renal cell carcinoma: an exploratory and comparative study | |
Alvarenga et al. | Assessing the combined performance of texture and morphological parameters in distinguishing breast tumors in ultrasound images | |
Caresio et al. | Quantitative analysis of thyroid tumors vascularity: A comparison between 3‐D contrast‐enhanced ultrasound and 3‐D Power Doppler on benign and malignant thyroid nodules | |
Huang et al. | Multi-parametric magnetic resonance imaging-based radiomics analysis of cervical cancer for preoperative prediction of lymphovascular space invasion | |
Muhtadi | Breast tumor classification using intratumoral quantitative ultrasound descriptors | |
Di Re et al. | MRI radiomics in the prediction of therapeutic response to neoadjuvant therapy for locoregionally advanced rectal cancer: a systematic review | |
Breit et al. | Revisiting DCE-MRI: classification of prostate tissue using descriptive signal enhancement features derived from DCE-MRI acquisition with high spatiotemporal resolution | |
CN106570325A (en) | Partial-least-squares-based abnormal detection method of mammary gland cell | |
Panic et al. | Normalization strategies in multi-center radiomics abdominal MRI: systematic review and meta-analyses | |
Muhtadi et al. | Texture quantified from ultrasound Nakagami parametric images is diagnostically relevant for breast tumor characterization | |
Antila et al. | Automatic segmentation for detecting uterine fibroid regions treated with MR‐guided high intensity focused ultrasound (MR‐HIFU) | |
Yang et al. | Machine learning based on magnetic resonance imaging and clinical parameters helps predict mesenchymal-epithelial transition factor expression in oral tongue squamous cell carcinoma: a pilot study | |
RU2698546C1 (en) | Method for prediction of volume of intraoperative blood loss in organ-preserving operations in tumors of renal parenchyma | |
Zaylaa et al. | Advancing Breast Cancer Diagnosis through Breast Mass Images, Machine Learning, and Regression Models | |
Gangeh et al. | Quantitative ultrasound spectroscopy and a kernel-based metric in clinical cancer response monitoring | |
Rautela et al. | A comprehensive review on computational techniques for breast cancer: past, present, and future | |
Lapointe et al. | OC-0417: Functional imaging using dual energy computed tomography and its application in radiation oncology | |
Dılek et al. | The success of machine learning algorithms developed with radiomic features obtained from preoperative contrast-enhanced MRI in the prediction of short-term survival in patients with glioblastoma |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application |
Ref document number: 14788804 Country of ref document: EP Kind code of ref document: A1 |
|
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2014788804 Country of ref document: EP |