Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1449 1674 (N50W09 L=334)   822    120        3 (1406/1632-1494/1710) -26.6 W04-W14 N46-N53
      2 1249 1402 (N23W02 L=328)   965    446        4 (1220/1250-1294/1614) -73.7 E03-W06 N13-N37
      3  699 1509 (N14E34 L=292)   825    221        2 ( 640/1436- 742/1612) -81.3 E36-E32 N08-N21
      4 1064 1547 (N26E13 L=312)   244    126        0 (1022/1494-1101/1596) -26.6 E17-E09 N23-N28
      5 1709 1449 (N43W38 L=  4)  2035    266        4 (1596/1360-1807/1526)  -4.4 W25-W51 N39-N45
      6 1545 1205 (N20W26 L=351)   464    216        0 (1473/1112-1582/1261) -31.0 W19-W30 N15-N22
      7 1696 1061 (N17W41 L=  6)   649    109        2 (1678/1005-1718/1124) -66.8 W38-W43 N13-N20
      8 1112  633 (S24W17 L=342)  1241    413        2 (1051/ 540-1237/ 803) -34.5 W10-W28 S28-S14
      9  706  653 (S35E09 L=316)   741    300        0 ( 584/ 639- 864/ 669)  29.1 E18-W00 S39-S29
     10  789  503 (S42W01 L=327)  1192    310        2 ( 650/ 477- 924/ 536)  31.4 E09-W09 S46-S35
     11 1593  420 (S21W62 L= 27)   861    151        1 (1548/ 367-1651/ 465)  82.9 W59-W65 S24-S16
  

    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