Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1645 1423 (N38W32 L=167)  1689    266        1 (1545/1328-1727/1485) -15.5 W22-W40 N33-N41
      2 1454 1373 (N28W17 L=152)   763    133        2 (1401/1322-1495/1414) -29.7 W12-W20 N25-N29
      3 1313 1361 (N23W08 L=143)   356    102        0 (1284/1314-1336/1410)  90.0 W06-W08 N18-N26
      4  365 1347 (S02E52 L= 83)  1492    202        3 ( 301/1257- 408/1398) -30.3 E59-E45 S07-N00
      5 1193 1266 (N13W03 L=138)   328    116        1 (1134/1244-1245/1290)  50.2 W00-W05 N10-N16
      6 1811 1178 (N28W51 L=186)   392    117        0 (1780/1121-1829/1222) -45.0 W47-W53 N24-N29
      7 1728 1096 (N20W43 L=178)   590     98        1 (1698/1051-1752/1133) -53.1 W40-W46 N17-N21
      8 1492  908 (N02W30 L=165)   752    212        3 (1466/ 814-1517/1025) -70.0 W27-W32 S03-N08
      9 1228  907 (S06W14 L=150)  1319    328        3 (1105/ 837-1307/ 965)  18.4 W08-W20 S10-S01
     10 1644  915 (N07W41 L=176)   279     72        0 (1606/ 900-1678/ 927)  18.4 W37-W43 N05-N08
     11  883  737 (S25E00 L=135)   392    163        0 ( 853/ 676- 903/ 813) -58.0 E04-W02 S28-S20
     12  532  788 (S31E25 L=110)   246     27        0 ( 510/ 771- 555/ 805)  33.7 E26-E23 S32-S29
     13 1730  578 (S08W63 L=198)   448     76        0 (1714/ 540-1748/ 612) -68.2 W61-W65 S10-S05
     14 1688  507 (S13W64 L=200)   833     92        1 (1640/ 474-1725/ 534)  35.5 W60-W68 S16-S11
     15 1478  332 (S30W60 L=195)  1847    180        3 (1415/ 287-1557/ 379)   2.4 W50-W74 S31-S27
     16 1142  194 (S51W52 L=187)   465     60        1 (1131/ 164-1159/ 229) -22.6 W46-W59 S53-S48
  

    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