Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  615 1485 (N31E24 L=193)   437     90        2 ( 577/1463- 660/1510) -23.2 E27-E20 N28-N32
      2 1600 1487 (N22W44 L=261)   334     82        0 (1563/1469-1632/1509)  18.4 W41-W47 N21-N23
      3 1253 1374 (N18W17 L=234)   306     97        1 (1201/1350-1309/1400) -26.6 W14-W20 N15-N19
      4  462 1339 (N23E35 L=182)   513    130        2 ( 417/1295- 496/1385) -39.8 E38-E32 N19-N25
      5 1481 1169 (N04W29 L=246)   404    101        1 (1465/1108-1499/1217) -90.0 W28-W30 N00-N06
      6 1740 1114 (S01W48 L=265)   964    185        1 (1698/1078-1808/1182) -29.1 W45-W54 S03-N02
      7  517 1084 (N06E31 L=186)   356     39        1 ( 488/1055- 543/1110)  26.6 E33-E29 N04-N07
      8 1672 1027 (S06W42 L=259)   233     55        1 (1654/1004-1687/1053) -90.0 W40-W43 S07-S04
      9 1763  965 (S10W50 L=266)   705    134        3 (1739/ 894-1797/1034)  63.4 W47-W53 S14-S06
     10 1594  902 (S13W36 L=252)   494     99        1 (1561/ 850-1620/ 940) -80.5 W33-W37 S16-S10
     11  316  756 (S11E51 L=166)  3079    617        8 ( 211/ 571- 423/ 895) -70.7 E60-E43 S22-S02
     12  674  847 (S09E23 L=194)   207     56        0 ( 659/ 812- 689/ 881)  78.7 E23-E22 S11-S06
     13 1504  799 (S19W29 L=246)   332    171        0 (1460/ 730-1542/ 845)  50.2 W25-W31 S22-S15
     14 1548  724 (S24W32 L=249)   954    246        1 (1450/ 603-1639/ 832)  47.3 W25-W39 S30-S17
     15  656  562 (S26E29 L=187)   428     44        1 ( 624/ 538- 684/ 586) -26.6 E31-E27 S27-S24
     16 1579  418 (S46W42 L=258)   915    109        2 (1553/ 379-1606/ 467)  90.0 W38-W45 S48-S41
     17 1483  402 (S46W32 L=249)   332     63        1 (1448/ 386-1516/ 416)  -9.5 W29-W35 S47-S44
  

    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