Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1346 1715 (N49E01 L=336)   801    107        2 (1310/1671-1388/1755) -32.0 E05-W03 N46-N52
      2  713 1636 (N22E39 L=299)   695    107        2 ( 682/1595- 756/1671)  59.0 E41-E35 N19-N25
      3 1081 1466 (N22E09 L=328)  1263    471        3 (1041/1277-1134/1669) -71.6 E14-E04 N09-N36
      4  904 1624 (N26E26 L=311)   308     82        1 ( 878/1581- 926/1668) -29.7 E30-E23 N24-N29
      5  557 1569 (N14E46 L=291)  1157    203        3 ( 493/1497- 616/1663) -80.5 E49-E44 N08-N20
      6 1598 1509 (N42W26 L=  3)  2126    302        4 (1485/1411-1708/1586)  -9.9 W14-W37 N39-N45
      7 1395 1248 (N18W15 L=353)   753    352        0 (1287/1108-1435/1351) -32.7 W06-W20 N10-N21
      8  198 1161 (S17E63 L=274)  1552    209        3 ( 160/1081- 225/1258)  84.8 E67-E60 S22-S11
      9 1554 1165 (N18W28 L=  5)   671    114        0 (1531/1104-1575/1223) -66.8 W25-W29 N14-N21
     10 1759  905 (N10W50 L= 28)   308     67        2 (1740/ 869-1786/ 939) -56.3 W48-W53 N08-N11
     11  943  721 (S24W04 L=342)  1593    422        2 ( 876/ 617-1084/ 883) -32.0 E02-W15 S28-S15
     12 1179  794 (S13W15 L=353)   245     69        1 (1166/ 749-1196/ 828) -53.1 W13-W16 S15-S11
     13  557  716 (S36E22 L=316)   873    275        1 ( 435/ 693- 696/ 740)  23.2 E32-E11 S39-S30
     14 1138  654 (S22W18 L=355)   281     81        0 (1094/ 635-1182/ 673)  45.0 W15-W19 S24-S20
     15  662  564 (S43E10 L=327)  1383    299        3 ( 537/ 533- 780/ 613)  30.1 E21-E02 S47-S35
     16 1617  461 (S18W61 L= 38)   533     95        0 (1570/ 443-1661/ 478)  33.7 W57-W64 S19-S15
     17 1483  450 (S23W50 L= 28)   298     54        0 (1461/ 432-1513/ 471)   0.0 W47-W53 S23-S21
  

    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