Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  842 1856 (N40E48 L=206)   616    102        0 ( 829/1806- 856/1904) -32.7 E56-E42 N35-N44
      2  692 1805 (N32E54 L=200)   486     91        1 ( 674/1767- 709/1845) -37.9 E58-E49 N28-N35
      3 1596 1672 (N56W25 L=278)   800     84        0 (1573/1637-1616/1717)  65.6 W21-W28 N51-N62
      4  660 1555 (N14E38 L=216)   660    203        1 ( 605/1452- 703/1618) -36.9 E44-E31 N08-N17
      5  370 1169 (S15E45 L=209)  2096    562        1 ( 187/ 960- 551/1295)  46.1 E61-E33 S30-S03
      6 1585 1220 (N21W28 L=282)   775    219        1 (1496/1138-1638/1275) -14.0 W20-W33 N17-N23
      7  171 1106 (S22E66 L=188)   548     68        1 ( 146/1071- 195/1144)  45.0 E69-E62 S24-S19
      8  407  986 (S25E38 L=216)   351    111        1 ( 366/ 950- 449/1022)  45.0 E41-E35 S27-S21
      9  597  817 (S30E20 L=233)  1305    368        2 ( 452/ 753- 771/ 863)   9.1 E32-E07 S31-S27
     10 1325  477 (S27W38 L=292)   382     93        1 (1284/ 440-1368/ 508)  68.2 W36-W39 S30-S24
     11 1545  332 (S27W67 L=321)   514     67        0 (1511/ 319-1582/ 344)  23.2 W62-W70 S27-S24
  

    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