Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  967 1766 (N37E26 L=271)   759     90        0 ( 932/1739-1023/1795)   8.1 E29-E22 N34-N39
      2  360 1607 (N14E64 L=233)  1835    265        1 ( 331/1543- 381/1680) -42.3 E69-E57 N09-N19
      3  452 1462 (N06E47 L=250)   769    182        2 ( 391/1421- 515/1508) -14.0 E54-E42 N04-N08
      4  875 1340 (N07E16 L=281)  1175    273        2 ( 843/1215- 905/1472) -60.3 E19-E11 N00-N15
      5 1678 1124 (N16W38 L=335)   934    140        0 (1639/1046-1713/1197) -45.0 W33-W41 N11-N19
      6  463 1095 (S15E36 L=261)   497    183        0 ( 387/1076- 539/1134)   0.0 E42-E30 S15-S13
      7  694 1071 (S12E20 L=277)   260     65        0 ( 656/1047- 735/1094) -11.3 E22-E17 S12-S11
      8  227 1025 (S23E55 L=242)  1592    335        0 ( 144/ 950- 323/1085)  31.0 E62-E46 S28-S18
      9 1637  864 (S00W41 L=338)   354     83        1 (1595/ 841-1673/ 881)  31.0 W38-W43 S02-N01
     10 1719  547 (S14W62 L=359)  1173    166        2 (1657/ 518-1761/ 590)  40.6 W56-W66 S16-S10
  

    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