Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  908 1722 (N55W05 L=105)   391    127        1 ( 846/1689- 961/1750)  13.0 E01-W12 N53-N57
      2 1431 1350 (N19W32 L=132)   665    239        0 (1367/1291-1502/1385) -45.0 W27-W35 N14-N22
      3 1150  656 (S18W02 L=102)   945    262        2 (1067/ 556-1227/ 743)  39.8 E04-W08 S23-S13
      4 1709  641 (S31W46 L=146)   578    103        1 (1662/ 610-1744/ 672)  14.0 W41-W49 S32-S29
      5 1480  636 (S25W25 L=125)   324     59        0 (1451/ 611-1502/ 664)  26.6 W22-W26 S26-S23
      6  791  505 (S23E24 L= 76)   713    241        0 ( 708/ 444- 877/ 583)  28.1 E31-E16 S26-S18
      7  387  546 (S16E55 L= 45)   458     69        0 ( 371/ 511- 409/ 578)  33.7 E57-E51 S18-S14
      8 1316  505 (S32W12 L=112)   503    179        0 (1252/ 444-1366/ 570)  37.9 W06-W15 S35-S27
      9 1393  449 (S38W18 L=117)   531    128        0 (1341/ 418-1449/ 485)  12.5 W13-W22 S39-S35
     10  894  422 (S30E19 L= 81)  1022    272        1 ( 786/ 371-1004/ 482)  14.7 E29-E10 S33-S27
     11  488  408 (S27E53 L= 47)   863     72        2 ( 458/ 371- 522/ 445)  63.4 E55-E49 S29-S23
  

    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