Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  861 1401 (N24E06 L= 22)   418    155        0 ( 798/1354- 926/1428) -32.0 E10-E02 N20-N26
      2 1776  806 (S22W52 L= 80)   398     95        0 (1755/ 754-1790/ 853)  71.6 W50-W53 S24-S18
      3 1713  651 (S31W48 L= 77)   314     41        1 (1693/ 628-1734/ 672) -36.9 W46-W50 S32-S29
      4  278  503 (S24E71 L=317)   716    128        2 ( 233/ 456- 305/ 560)  66.8 E76-E67 S26-S19
      5 1579  386 (S50W44 L= 73)   992    114        2 (1523/ 344-1634/ 420)   0.0 W36-W51 S53-S46
      6 1367  328 (S52W21 L= 50)   919    104        1 (1311/ 290-1412/ 378) -28.3 W15-W28 S56-S47
      7 1454  330 (S53W31 L= 60)   456     34        1 (1431/ 310-1482/ 352) -20.6 W27-W35 S54-S50
  

    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