Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1602 1536 (N21W54 L=133)  2014    359        2 (1460/1498-1718/1591) -33.7 W44-W64 N15-N28
      2  321 1416 (N46E48 L= 31)   573    114        2 ( 289/1387- 351/1438) -29.7 E51-E44 N43-N47
      3 1289 1284 (N16W22 L=101)   387    118        1 (1231/1249-1322/1328)   8.1 W17-W24 N14-N17
      4  316 1128 (N27E47 L= 32)   666    108        0 ( 276/1096- 359/1166) -39.8 E50-E42 N24-N29
      5  895  614 (S14E17 L= 62)   966    219        2 ( 844/ 518- 944/ 703)  45.0 E22-E11 S19-S09
      6  304  558 (S08E64 L= 15)  2446    319        6 ( 245/ 517- 393/ 625) -18.4 E71-E55 S11-S03
      7 1700  594 (S38W39 L=118)   639    130        1 (1665/ 562-1736/ 621) -45.0 W35-W43 S41-S34
      8 1702  533 (S43W40 L=119)   993    151        2 (1661/ 500-1738/ 571) -37.9 W35-W44 S46-S38
      9  643  310 (S30E47 L= 32)   553    161        1 ( 610/ 279- 663/ 363)  45.0 E50-E43 S32-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