Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1333 1553 (N38W10 L=153)   655    182        1 (1254/1502-1411/1602)  -4.4 W03-W16 N35-N40
      2  914 1499 (N25E17 L=126)   381    124        0 ( 861/1454- 966/1546) -24.0 E21-E12 N23-N27
      3 1678 1379 (N34W40 L=183)  2117    372        2 (1546/1311-1803/1458)  -7.4 W28-W52 N30-N37
      4 1405 1171 (N15W21 L=164)  1046    334        2 (1360/1043-1482/1254) -72.9 W16-W25 N07-N20
      5 1154  977 (S02W08 L=151)   422    141        0 (1071/ 950-1208/1014)  39.8 W04-W11 S04-N01
      6  507  703 (S31E29 L=114)   306     59        1 ( 485/ 672- 531/ 733)  68.2 E30-E27 S33-S28
      7 1160  692 (S18W15 L=158)   225     34        0 (1139/ 668-1181/ 713) -33.7 W13-W16 S19-S17
      8 1624  657 (S11W48 L=191)   731    178        2 (1566/ 630-1681/ 691)  -5.2 W42-W53 S12-S09
      9 1250  251 (S46W42 L=185)   360     43        1 (1226/ 235-1269/ 273)  90.0 W39-W44 S47-S44
  

    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