Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  752 1811 (N42E40 L= 92)  1361    124        3 ( 706/1757- 809/1853) -17.4 E49-E33 N37-N45
      2 1017 1683 (N36E12 L=120)   337     56        2 ( 994/1650-1042/1713) -36.9 E14-E10 N33-N37
      3  484 1652 (N26E52 L= 80)   642    104        1 ( 439/1620- 527/1674)  -9.5 E58-E46 N24-N27
      4 1537 1617 (N42W30 L=162)   620     87        2 (1510/1582-1571/1649)  40.6 W26-W33 N38-N44
      5 1729 1228 (N18W44 L=176)   720    174        0 (1692/1155-1791/1286) -36.9 W40-W50 N15-N21
      6  133  996 (S17E70 L= 63)  4339    273        7 (  78/ 905- 207/1081)  29.4 E84-E59 S22-S11
      7  625 1059 (S09E25 L=108)   261     65        0 ( 596/1036- 663/1078)  26.6 E26-E22 S10-S08
      8 1564  995 (N01W34 L=166)   397    137        0 (1504/ 968-1624/1023)   7.1 W29-W37 N00-N03
      9 1868  940 (N05W61 L=194)   516     76        0 (1843/ 900-1886/ 978) -33.7 W57-W64 N02-N06
     10 1502  911 (S05W31 L=163)   264    107        0 (1444/ 898-1546/ 927)   7.1 W26-W34 S05-S03
     11  415  871 (S23E39 L= 93)   431     68        1 ( 391/ 834- 436/ 907) -68.2 E40-E37 S25-S20
     12  891  874 (S17E06 L=127)   297     93        0 ( 848/ 858- 922/ 903)  36.9 E08-E03 S18-S14
     13  171  783 (S30E68 L= 64)  1529    184        2 ( 150/ 699- 192/ 874)  79.7 E73-E65 S35-S24
     14 1526  677 (S18W39 L=171)   373     74        0 (1490/ 646-1556/ 711)  53.1 W36-W40 S19-S15
     15 1560  572 (S23W46 L=178)   481     81        1 (1526/ 550-1604/ 597)   0.0 W41-W49 S24-S22
     16  607  418 (S52E26 L=107)   842    133        0 ( 579/ 370- 651/ 472) -49.4 E28-E21 S55-S47
     17 1528  440 (S32W51 L=183)   385     65        1 (1502/ 410-1553/ 465) -45.0 W47-W53 S33-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