Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  555 1644 (N17E49 L= 83)   856    249        2 ( 522/1547- 610/1737) -66.8 E52-E44 N10-N24
      2  908 1683 (N29E28 L=104)   319    103        0 ( 883/1649- 932/1723) -39.8 E31-E25 N26-N31
      3  626 1624 (N17E43 L= 89)   524    224        1 ( 601/1561- 668/1695) -66.0 E46-E37 N13-N22
      4  687 1526 (N13E34 L= 98)   608    211        1 ( 652/1436- 713/1614) -61.4 E37-E31 N07-N18
      5  465 1318 (S04E42 L= 91)   919    240        2 ( 423/1208- 493/1409) -51.3 E48-E38 S10-N00
      6 1223 1345 (N17W02 L=135)   242    100        0 (1196/1288-1247/1382)  90.0 W01-W03 N12-N19
      7 1535 1183 (N17W25 L=157)   401    122        3 (1515/1135-1557/1252) -56.3 W22-W27 N13-N20
      8  293 1104 (S20E50 L= 83)   313     84        1 ( 272/1063- 309/1150)  80.5 E51-E48 S22-S16
      9 1360 1046 (N04W18 L=150)   533    234        2 (1258/1006-1477/1102)  45.0 W12-W23 N00-N10
     10  773 1035 (S13E14 L=118)   189     83        0 ( 741/ 993- 798/1078) -26.6 E17-E11 S14-S10
     11  129 1036 (S26E65 L= 67)   854    163        0 ( 107/ 993- 162/1073)  51.3 E68-E62 S28-S23
     12  577  994 (S20E26 L=107)   430    161        1 ( 510/ 933- 641/1049)  52.1 E29-E22 S24-S15
     13 1567  784 (S04W39 L=171)   264     76        0 (1527/ 770-1599/ 796)  36.9 W36-W41 S05-S02
     14 1027  673 (S26W11 L=143)   422    113        2 ( 981/ 637-1097/ 704)  56.3 W07-W14 S28-S22
     15  879  622 (S33W03 L=136)   288     46        2 ( 859/ 599- 907/ 649)  90.0 W01-W05 S35-S31
     16 1652  483 (S16W60 L=192)  1899    364        7 (1545/ 420-1764/ 536)  45.0 W54-W69 S22-S08
     17 1329  284 (S38W50 L=182)  1694    424        3 (1245/ 203-1414/ 373)  -5.0 W39-W62 S39-S35
  

    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