Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  860 1770 (N36E31 L= 90)   834    143        5 ( 800/1720- 921/1802) -17.1 E37-E24 N33-N37
      2 1041 1604 (N28E12 L=109)   345    105        2 (1007/1569-1095/1631)  51.3 E13-E09 N25-N30
      3 1389 1574 (N35W12 L=133)   929    205        4 (1295/1517-1455/1605)  -4.8 W04-W16 N32-N36
      4 1097 1349 (N13E03 L=119)   226    103        0 (1068/1300-1130/1393) -45.0 E05-W00 N10-N15
      5  478 1254 (S04E38 L= 83)   376    143        0 ( 440/1198- 496/1304)  81.9 E40-E36 S07-N00
      6 1655 1178 (N16W36 L=157)   579    131        0 (1600/1135-1700/1223) -14.0 W31-W39 N14-N17
      7  189 1159 (S13E62 L= 59)   574     68        3 ( 153/1132- 226/1184)   0.0 E67-E57 S14-S11
      8  722 1070 (S10E18 L=103)   521    231        0 ( 626/1013- 800/1131)  41.6 E23-E14 S14-S05
      9 1886 1071 (N17W60 L=181)   516    105        2 (1870/1025-1904/1128)  74.1 W57-W62 N13-N20
     10 1714  949 (N05W45 L=166)   781    184        1 (1692/ 860-1730/1047) -71.6 W41-W47 N00-N09
     11  194  998 (S22E60 L= 61)   821    117        0 ( 148/ 953- 242/1038)  22.6 E66-E54 S24-S19
     12 1180  894 (S11W12 L=133)   247     85        0 (1153/ 853-1199/ 928) -53.1 W09-W12 S13-S09
     13 1217  792 (S16W16 L=137)   237     75        1 (1192/ 756-1238/ 824)  78.7 W14-W17 S18-S13
     14  919  767 (S24E01 L=120)   304    109        0 ( 875/ 748- 970/ 802)   0.0 E04-W02 S24-S22
     15 1700  683 (S10W52 L=173)  1245    167        6 (1639/ 591-1749/ 763) -59.0 W46-W56 S15-S04
     16 1756  641 (S10W60 L=181)   742    128        1 (1717/ 597-1795/ 695)  74.1 W56-W62 S13-S06
     17 1832  642 (S08W69 L=190)  1107     88        3 (1804/ 597-1871/ 688) -68.2 W65-W75 S10-S05
     18 1638  475 (S23W59 L=180)  1381    221        0 (1528/ 450-1728/ 514)  33.7 W50-W65 S27-S17
     19 1472  310 (S37W59 L=180)   751    117        2 (1400/ 286-1517/ 345)  29.1 W53-W63 S39-S34
  

    OPENR, u, binfile, /COMP, /GET_LUN
    head = FLTARR(6) ;[xsize,ysize,radius,hh,mm,ss]
    READU, u, head
    data = INTARR(FIX(head[0]),FIX(head[1]))
    READU, u, data
    FREE_LUN, u
    

• each filament has the value of Fil_ID
• all values are in pixel

  c = WHERE(data EQ Fil_ID, area)
  f = '(2(A,I6))'
  PRINT, FORM=f, "Area of Fil_ID", Fil_ID," =", area 

  data = ROTATE(data,7)
  Solar_y = 2048-Solar_y