Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1258 1808 (N55E13 L=243)   639    113        1 (1204/1785-1308/1827)  16.7 E17-E07 N51-N57
      2  896 1785 (N38E37 L=220)   712    156        0 ( 815/1761- 958/1803)   8.7 E43-E30 N36-N40
      3  305 1523 (N07E68 L=189)  2272    173        4 ( 249/1457- 373/1577)  52.1 E74-E62 N02-N11
      4 1343 1341 (N22W10 L=267)   438    153        0 (1309/1300-1397/1395) -15.9 W07-W14 N20-N25
      5  233 1290 (S08E63 L=194)   620     81        1 ( 188/1270- 266/1313)  18.4 E68-E58 S09-S06
      6 1865  880 (N12W63 L=320)  1035    113        2 (1835/ 828-1889/ 936)  74.1 W60-W66 N08-N15
      7 1358  566 (S21W34 L=291)   366     98        0 (1343/ 512-1371/ 613) -50.2 W32-W37 S23-S17
      8  762  546 (S40E03 L=254)   364     93        0 ( 740/ 504- 797/ 593) -33.7 E05-W01 S41-S37
      9 1409  497 (S23W41 L=298)   308     57        0 (1384/ 481-1439/ 516)   0.0 W38-W43 S24-S21
  

    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