Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1015 1588 (N34W04 L= 11)  4145    785        9 ( 713/1375-1302/1741) -36.9 E20-W20 N17-N49
      2 1498 1554 (N28W39 L= 46)  2164    467        0 (1346/1426-1671/1714) -55.7 W31-W49 N18-N41
      3 1564 1670 (N36W52 L= 59)   596    126        0 (1496/1647-1615/1697) -21.8 W46-W57 N34-N38
      4  433 1528 (N35E42 L=324)  1344    334        1 ( 296/1486- 571/1569) -11.6 E63-E29 N30-N39
      5  751 1502 (N30E15 L=352)   540    171        2 ( 686/1449- 821/1538) -33.7 E20-E10 N26-N33
      6  779 1353 (N20E13 L=354)   415    179        0 ( 733/1326- 830/1407) -38.7 E16-E09 N18-N23
      7  257 1255 (N18E54 L=313)  2592    449        8 ( 150/1092- 401/1398) -28.8 E80-E39 N07-N29
      8  512 1312 (N20E32 L=335)   267     81        1 ( 473/1282- 546/1337) -31.0 E34-E29 N17-N21
      9  741 1168 (N09E16 L=351)   348    204        0 ( 656/1151- 828/1208) -15.3 E21-E10 N07-N11
     10  175 1149 (N13E63 L=304)   463     80        1 ( 154/1105- 186/1186) -71.6 E65-E60 N09-N15
     11  936  876 (S10E06 L=  0)  1336    382        5 ( 878/ 700-1031/1017)  58.6 E11-W00 S20-S01
     12  368  728 (S14E48 L=319)   911    180        0 ( 335/ 643- 409/ 815)  42.3 E53-E42 S19-S09
     13 1125  742 (S19W04 L= 11)   438    113        0 (1094/ 690-1169/ 791)  50.2 W02-W07 S22-S16
     14 1764  739 (S24W53 L= 60)   318     72        0 (1744/ 713-1788/ 772) -33.7 W50-W56 S25-S21
     15  262  569 (S23E67 L=300)  1941    170        5 ( 229/ 485- 297/ 647)  45.0 E74-E60 S28-S18
     16  571  426 (S36E42 L=325)   557    112        1 ( 530/ 380- 608/ 475)  28.6 E48-E36 S38-S32
     17  434  335 (S41E66 L=300)   909     79        0 ( 389/ 318- 468/ 359)   0.0 E76-E58 S41-S39
  

    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