Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1633 1355 (N33W31 L=124)   493     68        2 (1609/1321-1655/1391) -36.9 W29-W33 N30-N35
      2 1813 1100 (N23W52 L=145)  2486    374        4 (1720/ 967-1909/1185)  46.6 W44-W67 N12-N30
      3  925 1106 (S02E08 L= 85)   237     47        0 ( 898/1080- 943/1132) -33.7 E09-E06 S03-S01
      4  407  770 (S35E36 L= 57)  2448    333        8 ( 319/ 665- 494/ 915)  52.7 E46-E30 S44-S23
      5 1762  649 (S04W62 L=155)   564     50        0 (1742/ 617-1778/ 677) -45.0 W59-W64 S05-S02
      6 1388  599 (S18W35 L=128)   668    179        0 (1317/ 565-1456/ 665)  60.9 W31-W37 S21-S12
      7 1077  546 (S30W18 L=111)  1192    385        3 ( 944/ 461-1222/ 623)  57.1 W12-W24 S39-S22
      8 1459  553 (S18W42 L=134)   622    116        1 (1408/ 521-1510/ 579)  56.3 W39-W44 S21-S15
      9 1649  507 (S15W60 L=153)   430     91        0 (1612/ 483-1687/ 542)  63.4 W57-W62 S17-S11
  

    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