Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  298 1181 (S13E51 L=357)  1556    192        2 ( 246/1095- 363/1244)  60.9 E56-E46 S18-S09
      2  654 1035 (S15E22 L= 26)   546    187        0 ( 601/ 951- 702/1118) -37.9 E26-E17 S18-S11
      3 1721 1062 (N14W42 L= 90)   319     68        1 (1698/1026-1734/1102) -63.4 W40-W43 N11-N16
      4 1341  971 (S02W19 L= 67)   969    235        0 (1237/ 934-1447/1048)   0.0 W11-W25 S04-N00
      5 1860  999 (N15W57 L=105)   526    116        1 (1829/ 977-1885/1027)   8.1 W52-W59 N13-N17
      6 1893  937 (N13W62 L=111)   532    126        3 (1867/ 899-1913/ 965) -18.4 W59-W66 N11-N14
      7 1766  834 (N03W52 L=101)   940    162        1 (1735/ 756-1792/ 908) -54.5 W49-W56 N00-N07
      8  398  837 (S32E37 L= 11)  1082    325        0 ( 262/ 811- 508/ 864)  10.8 E49-E28 S35-S29
      9 1167  720 (S20W16 L= 65)   278    107        0 (1114/ 704-1222/ 735)  26.6 W13-W19 S22-S19
     10 1419  673 (S16W33 L= 81)   652    245        0 (1319/ 648-1526/ 721)  30.3 W26-W38 S18-S10
     11 1690  540 (S15W59 L=108)  3438    550       10 (1611/ 471-1774/ 646) -74.7 W50-W72 S19-S08
     12 1129  557 (S31W20 L= 68)   477    133        0 (1057/ 536-1196/ 581)  36.9 W15-W23 S33-S27
  

    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