Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1047 1620 (N39W22 L=161)   795    196        0 (1003/1534-1096/1688)  33.7 W16-W28 N33-N41
      2 1466 1586 (N23W46 L=185)  1057    190        3 (1395/1541-1550/1618) -37.9 W40-W50 N18-N26
      3  498 1238 (N31E28 L=112)   890    346        0 ( 369/1207- 650/1274) -22.2 E38-E16 N26-N36
      4 1513 1149 (S01W31 L=171)   893    217        3 (1421/1115-1600/1182) -28.6 W24-W36 S05-N01
      5 1666  967 (S16W37 L=176)   680    104        3 (1618/ 932-1709/1001)   0.0 W33-W41 S17-S14
      6 1401  797 (S17W15 L=155)   347    117        2 (1386/ 756-1429/ 855)  81.9 W14-W17 S20-S13
      7 1071  678 (S15E07 L=133)   644    163        0 ( 979/ 643-1145/ 705) -41.2 E11-E03 S19-S11
      8 1169  567 (S24E04 L=135)   307     83        2 (1155/ 524-1184/ 613)  71.6 E05-E02 S27-S21
      9  692  461 (S18E37 L=103)   669    278        4 ( 613/ 402- 798/ 509)   3.6 E44-E28 S20-S15
     10  767  398 (S24E35 L=104)   400     72        2 ( 740/ 365- 797/ 435) -90.0 E36-E33 S27-S21
     11 1022  214 (S46E31 L=108)  1006    191        3 ( 984/ 167-1053/ 258)  32.5 E36-E25 S49-S41
  

    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