Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1341 1716 (N49E02 L=336)   745    101        1 (1306/1673-1384/1750) -32.0 E05-W03 N46-N51
      2 1078 1489 (N23E10 L=328)  1201    406        2 (1039/1320-1131/1677) -70.2 E15-E06 N12-N37
      3  708 1637 (N22E39 L=299)   652    102        2 ( 679/1598- 748/1671)  90.0 E42-E36 N19-N25
      4  554 1576 (N14E47 L=291)  1269    200        1 ( 490/1499- 612/1666) -81.3 E50-E44 N08-N21
      5 1594 1511 (N43W25 L=  3)  2127    290        2 (1480/1416-1704/1589) -10.3 W14-W37 N39-N45
      6 1355 1324 (N21W11 L=349)   284    120        0 (1298/1287-1396/1355)   0.0 W06-W14 N19-N21
      7  195 1178 (S16E64 L=274)  2413    221        3 ( 146/1082- 238/1267)  85.2 E73-E61 S22-S10
      8 1412 1202 (N16W17 L=356)   304    137        0 (1394/1120-1426/1254) -60.3 W15-W19 N11-N18
      9 1547 1166 (N18W27 L=  5)   640    106        0 (1524/1104-1565/1223) -66.8 W24-W29 N14-N21
     10 1756  904 (N10W50 L= 28)   457     67        2 (1735/ 865-1786/ 944) -53.1 W47-W53 N07-N12
     11  934  728 (S24W03 L=342)  1578    416        1 ( 870/ 626-1076/ 894) -36.3 E02-W14 S28-S14
     12 1171  796 (S13W15 L=353)   230     76        1 (1159/ 752-1185/ 832) -53.1 W13-W16 S15-S11
     13  543  718 (S36E23 L=315)   830    260        0 ( 434/ 695- 678/ 738)  25.3 E32-E13 S39-S30
     14 1128  658 (S22W17 L=355)   282     84        0 (1085/ 641-1173/ 676)  38.7 W14-W19 S24-S20
     15  657  572 (S42E11 L=327)  1371    298        3 ( 532/ 539- 797/ 629)  34.4 E21-E02 S47-S33
     16 1667  496 (S14W63 L= 41)  1024    249        0 (1574/ 446-1776/ 560)  42.7 W57-W71 S19-S07
     17 1477  454 (S23W49 L= 28)   290     60        0 (1455/ 434-1506/ 474)   0.0 W46-W52 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