Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  587 1642 (N26E48 L=253)   549    126        1 ( 555/1606- 628/1673) -68.2 E50-E43 N23-N28
      2  842 1509 (N23E23 L=278)   247    102        1 ( 817/1464- 861/1557) -39.8 E26-E20 N20-N25
      3 1199 1449 (N27W01 L=302)   936    303        1 (1074/1356-1311/1507) -10.0 E07-W10 N23-N30
      4 1373 1430 (N30W14 L=314)   264    106        0 (1339/1385-1408/1475) -33.7 W10-W16 N28-N32
      5 1873 1185 (N27W64 L=  4)  2224    228        1 (1832/1085-1897/1290)  81.3 W57-W68 N20-N34
      6  441 1066 (S11E38 L=263)   477     85        2 ( 400/1039- 474/1090)  38.7 E40-E35 S13-S09
      7 1758  990 (N12W50 L=351)   612     94        0 (1734/ 941-1780/1035)  80.5 W48-W52 N08-N14
      8 1438  549 (S21W39 L=339)  2082    423        3 (1362/ 432-1553/ 646)  82.9 W34-W49 S28-S12
      9  466  556 (S44E33 L=267)   343     50        1 ( 451/ 527- 485/ 582)  78.7 E35-E31 S46-S41
     10  914  535 (S34W04 L=305)   671    222        0 ( 820/ 505-1016/ 562)  26.6 E02-W10 S37-S30
  

    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