Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  884 1703 (N47E18 L=252)   311     40        0 ( 860/1683- 914/1717)   0.0 E20-E15 N45-N48
      2 1204 1626 (N43W11 L=281)  2695    554        6 ( 942/1573-1423/1679)  -2.9 E11-W29 N39-N46
      3 1735 1574 (N41W69 L=339)   789     49        0 (1706/1558-1768/1586)   0.0 W63-W81 N40-N42
      4  340 1528 (N29E62 L=208)   431     58        1 ( 324/1503- 363/1555) -53.1 E64-E59 N27-N31
      5 1697 1508 (N36W56 L=326)  1039    103        0 (1672/1456-1740/1543)  21.8 W50-W65 N32-N38
      6 1531 1495 (N35W37 L=307)   268     44        0 (1511/1470-1554/1515)  31.0 W34-W39 N33-N36
      7 1812 1299 (N22W61 L=331)   808    154        1 (1784/1229-1835/1383)  78.7 W58-W63 N17-N27
      8 1758 1303 (N23W54 L=325)   706    117        0 (1726/1250-1783/1354) -63.4 W52-W56 N19-N25
      9 1130 1164 (N11W06 L=276)  1847    430        2 (1078/ 951-1157/1344) -85.0 W01-W07 S01-N22
     10  426 1145 (N06E40 L=230)  1628    336        1 ( 377/1000- 466/1312) -67.6 E46-E36 S02-N15
     11  687 1258 (N15E23 L=247)   219     87        0 ( 647/1240- 726/1276)  11.3 E25-E20 N13-N16
     12  595 1180 (N09E28 L=242)   562    162        0 ( 571/1117- 631/1250) -69.4 E30-E25 N05-N13
     13  669  771 (S15E21 L=249)  1073    321        1 ( 604/ 649- 707/ 922)  71.6 E26-E18 S23-S05
     14 1644  453 (S32W57 L=327)  1500    107        4 (1587/ 419-1699/ 481)  -6.7 W48-W65 S34-S30
     15 1449  372 (S39W41 L=311)   580     66        1 (1412/ 354-1489/ 395)   0.0 W37-W44 S40-S37
     16 1077  303 (S47W10 L=280)   302     48        0 (1051/ 281-1109/ 328)  45.0 W08-W12 S49-S44
     17  765  291 (S51E19 L=251)   348     52        0 ( 746/ 265- 785/ 320) -68.2 E20-E17 S53-S48
     18  581  277 (S55E44 L=226)  1235    131        1 ( 509/ 253- 632/ 311)  10.9 E62-E34 S58-S50
     19  671  288 (S52E30 L=240)   457     41        0 ( 653/ 268- 693/ 311)  39.8 E33-E27 S54-S49
  

    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