Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  678 1741 (N30E47 L=244)   419    110        1 ( 647/1713- 704/1772) -23.2 E50-E43 N28-N31
      2  238 1259 (S09E61 L=231)  1505    267        4 ( 205/1147- 266/1379) -66.8 E66-E56 S16-S02
      3  315  973 (S25E46 L=246)  2002    523        6 ( 171/ 931- 462/1024)  17.1 E59-E33 S30-S21
      4  793  932 (S15E11 L=281)  1201    338        0 ( 636/ 903- 956/ 973)  16.7 E21-E01 S18-S11
      5  126  891 (S33E64 L=227)   847    201        2 (  99/ 867- 155/ 915) -18.4 E67-E60 S34-S31
      6 1702  575 (S10W60 L=351)  1150    166        1 (1633/ 517-1774/ 642)  57.5 W55-W63 S14-S03
      7 1371  389 (S30W44 L=336)   303     59        0 (1339/ 378-1397/ 411)  45.0 W42-W46 S31-S28
      8 1113  384 (S39W28 L=319)   490     67        0 (1063/ 366-1150/ 402)  38.7 W24-W30 S41-S37
  

    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