Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1136 1822 (N56E17 L= 36)  1915    250        2 (1033/1772-1225/1853)  51.3 E25-E07 N48-N63
      2 1709 1116 (N18W44 L= 97)  3464    570       15 (1568/ 990-1825/1204)   4.6 W31-W56 N10-N25
      3  212 1134 (S11E61 L=352)   830    134        2 ( 175/1069- 260/1183)  45.0 E65-E55 S15-S07
      4 1556 1101 (N14W32 L= 85)   347     70        1 (1535/1069-1570/1139)  90.0 W30-W33 N11-N16
      5 1889 1019 (N17W66 L=119)   566     75        1 (1872/ 991-1909/1058) -26.6 W63-W70 N14-N18
      6  205  726 (S37E65 L=348)  1247    117        0 ( 157/ 691- 258/ 763)  12.1 E83-E54 S39-S33
      7 1198  631 (S21W20 L= 73)   511    145        0 (1131/ 593-1255/ 677)  54.5 W16-W22 S24-S17
      8 1702  650 (S09W56 L=109)   397     51        2 (1678/ 630-1732/ 672)  36.9 W53-W58 S10-S07
      9 1716  521 (S16W66 L=119)  2020    220        3 (1648/ 446-1801/ 575)  50.7 W60-W77 S22-S11
     10 1374  424 (S30W41 L= 94)  1728    443        1 (1184/ 370-1547/ 487)  32.5 W27-W51 S37-S22
     11 1144  241 (S50W38 L= 91)   510     60        1 (1112/ 218-1174/ 266) -11.3 W32-W42 S51-S48
     12 1313  193 (S49W62 L=115)  1081     92        4 (1267/ 157-1344/ 220)  -4.2 W51-W79 S50-S46
  

    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