Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1673 1403 (N08W52 L=312)  1054    282        4 (1609/1339-1719/1483)  48.8 W46-W58 N03-N11
      2  509  936 (N10E33 L=228)   251     63        0 ( 479/ 910- 531/ 954) -45.0 E34-E31 N08-N11
      3  667  822 (N00E25 L=235)   534    132        0 ( 598/ 800- 729/ 845) -23.2 E29-E21 S01-N03
      4  280  773 (N05E55 L=205)  1136    128        3 ( 241/ 739- 334/ 816) -35.5 E59-E50 N03-N08
      5 1893  746 (S37W58 L=319)   713    171        0 (1870/ 730-1925/ 772) -32.0 W55-W63 S39-S34
      6  474  468 (S16E51 L=210)  1075    187        3 ( 425/ 426- 535/ 526) -69.4 E55-E47 S19-S11
      7  400  437 (S17E59 L=201)  2126    277        5 ( 347/ 370- 469/ 509) -84.3 E64-E55 S21-S11
      8  728  321 (S31E41 L=220)   293     73        0 ( 712/ 296- 743/ 349)  31.0 E42-E37 S33-S29
  

    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