Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  334 1610 (N11E70 L= 22)  1007    224        1 ( 296/1563- 357/1646) -26.6 E75-E64 N08-N13
      2  415 1567 (N10E57 L= 34)   538    111        0 ( 381/1546- 452/1596)  36.9 E60-E54 N08-N12
      3 1070 1445 (N18E09 L= 82)   567    223        0 (1050/1350-1092/1549) -63.4 E12-E05 N12-N24
      4 1439 1456 (N31W13 L=104)   432     49        1 (1411/1417-1464/1490) -53.1 W10-W14 N28-N32
      5 1265 1307 (N16W06 L= 97)   459    159        0 (1241/1235-1284/1386) -58.0 W02-W08 N11-N20
      6  171 1250 (S13E68 L= 23)   643    101        1 ( 142/1219- 192/1289)  68.2 E70-E65 S15-S10
      7 1796 1094 (N21W47 L=138)   405    130        0 (1752/1078-1832/1112)  29.7 W43-W50 N18-N23
      8 1072 1022 (S05W03 L= 94)   224     67        0 (1050/ 986-1090/1063)  90.0 W01-W03 S07-S02
      9  404  857 (S32E36 L= 55)  2521    608        7 ( 218/ 746- 554/1012)  42.7 E53-E26 S43-S19
     10 1769  905 (N10W50 L=141)   392    129        1 (1731/ 887-1800/ 944)  45.0 W46-W52 N08-N12
     11 1804  756 (N03W60 L=151)   958    312        0 (1702/ 723-1879/ 811)  26.6 W50-W66 S01-N08
     12  772  700 (S32E06 L= 85)   772    248        2 ( 699/ 610- 862/ 808)  76.8 E08-E02 S39-S22
     13 1113  702 (S22W15 L=106)   539    215        0 (1001/ 655-1197/ 768)  56.3 W09-W17 S27-S15
     14 1389  668 (S15W32 L=123)   347    180        0 (1348/ 639-1429/ 720)  -8.1 W27-W34 S17-S13
     15 1769  646 (S03W62 L=153)   484     66        1 (1746/ 608-1792/ 674) -76.0 W60-W64 S05-S01
     16 1141  545 (S29W23 L=114)   246     93        0 (1094/ 534-1195/ 565)  45.0 W19-W24 S31-S26
     17 1647  501 (S15W59 L=150)   452    105        0 (1604/ 475-1697/ 545)  60.3 W55-W60 S17-S10
     18 1398  362 (S31W49 L=140)  1336    250        3 (1305/ 328-1496/ 396)  45.0 W42-W54 S35-S25
  

    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