Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  553 1734 (N36E59 L=247)   767    107        2 ( 526/1690- 582/1783) -18.4 E69-E51 N32-N39
      2  713 1609 (N30E35 L=271)   321     99        1 ( 675/1576- 747/1639) -18.4 E39-E30 N28-N31
      3 1145 1539 (N34E02 L=304)   255     45        1 (1129/1510-1168/1563) -76.0 E03-W00 N31-N35
      4 1232 1474 (N31W05 L=312)   479    155        0 (1175/1411-1287/1517) -26.6 W00-W10 N27-N33
      5  888 1417 (N21E16 L=290)   324    110        0 ( 843/1382- 926/1457) -23.2 E19-E12 N19-N22
      6 1649  685 (S10W49 L=355)   705    138        1 (1591/ 657-1720/ 712)  29.1 W44-W55 S12-S07
      7  760  620 (S30E10 L=296)   781    264        1 ( 650/ 584- 885/ 650)  22.4 E18-E01 S34-S26
      8 1227  467 (S32W26 L=332)  1657    274        3 (1169/ 368-1299/ 567) -32.9 W18-W36 S37-S25
      9 1437  382 (S33W47 L=353)   908    102        3 (1387/ 354-1483/ 413)  -5.2 W41-W52 S34-S31
  

    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