Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  786 1357 (N27E07 L=336)   258     78        0 ( 752/1334- 828/1381) -51.3 E08-E04 N24-N29
      2 1841 1212 (S06W60 L= 42)  2632    344        7 (1799/1058-1881/1340)  51.8 W54-W68 S14-N01
      3 1699 1177 (S04W45 L= 28)  2075    416        0 (1591/ 985-1760/1306) -76.6 W40-W48 S16-N05
      4  777 1249 (N21E09 L=333)   447    208        1 ( 678/1227- 863/1269) -14.0 E15-E03 N18-N22
      5 1130 1213 (N11W10 L=352)   241     94        0 (1088/1194-1180/1231) -36.9 W08-W12 N09-N12
      6 1844  877 (S25W54 L= 37)   926    131        3 (1820/ 800-1870/ 941)  82.9 W51-W57 S29-S21
      7 1308  874 (S12W13 L=355)   393    100        1 (1287/ 833-1339/ 917)  53.1 W11-W14 S14-S09
      8 1526  813 (S20W26 L=  8)   232     80        1 (1503/ 775-1554/ 846) -63.4 W24-W27 S23-S17
      9 1042  540 (S25E10 L=332)  4432   1150        4 ( 730/ 174-1428/ 827)  19.4 E55-W19 S44-S17
     10  885  732 (S11E14 L=328)   299    178        0 ( 805/ 690- 932/ 785)  18.4 E19-E10 S11-S08
     11 1819  703 (S36W53 L= 35)   448     61        1 (1798/ 672-1842/ 739) -40.6 W49-W57 S38-S32
     12  457  460 (S18E53 L=290)   703    136        3 ( 419/ 411- 510/ 512) -82.9 E56-E50 S21-S13
  

    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