Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1557 1625 (N32W54 L= 34)   593    194        0 (1514/1592-1589/1673)  16.7 W48-W58 N30-N34
      2 1285 1495 (N30W27 L=  7)   760    173        3 (1223/1448-1366/1534)  -5.2 W22-W33 N27-N31
      3  375 1477 (N44E47 L=294)   353     91        0 ( 346/1458- 414/1490)   8.1 E49-E42 N43-N45
      4 1744 1441 (N17W61 L= 41)  1753    242        2 (1664/1404-1819/1487) -26.6 W55-W70 N12-N21
      5 1289 1207 (N13W19 L=359)   213     43        1 (1265/1187-1310/1230)  33.7 W17-W20 N11-N14
      6  690 1134 (N18E19 L=321)   668    216        0 ( 582/1084- 773/1195) -36.9 E26-E14 N13-N22
      7 1376  957 (S03W20 L=  0)  1960    537        6 (1320/ 761-1415/1157)  75.4 W14-W23 S14-N09
      8 1608  997 (S06W37 L= 17)   296     72        0 (1583/ 961-1636/1021) -90.0 W34-W38 S07-S03
      9  143  935 (N10E70 L=270)  1035    263        1 ( 103/ 894- 174/ 996)  23.6 E81-E64 N07-N14
     10 1524  876 (S11W29 L=  9)   276     70        0 (1508/ 840-1541/ 916)  68.2 W27-W30 S13-S08
     11 1570  702 (S23W31 L= 11)  1884    262        1 (1533/ 592-1596/ 838)  72.6 W27-W33 S30-S14
     12 1122  369 (S38E05 L=335)  1387    315        2 (1012/ 287-1190/ 479)  39.1 E16-W00 S43-S30
  

    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