Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  794 1798 (N62W15 L=229)  1110    210        2 ( 726/1736- 863/1848)   9.9 W04-W27 N58-N64
      2 1562 1715 (N27W66 L=280)  1414    169        4 (1519/1655-1620/1777)  56.3 W61-W71 N22-N31
      3 1476 1659 (N27W52 L=267)  1821    355        2 (1324/1604-1648/1693) -33.7 W41-W63 N18-N32
      4 1537 1419 (N13W42 L=256)   326     68        1 (1518/1377-1561/1450)  31.0 W39-W45 N10-N14
      5 1415 1369 (N14W32 L=246)   412    197        0 (1305/1347-1487/1385) -20.0 W25-W36 N10-N16
      6 1763 1133 (S09W49 L=263)   474     83        0 (1720/1105-1805/1161) -39.8 W45-W52 S11-S06
      7  192 1080 (N28E59 L=155)   687     83        0 ( 157/1043- 227/1116) -29.1 E64-E55 N25-N30
      8 1606  742 (S27W28 L=242)   833    204        0 (1546/ 683-1694/ 813)  12.1 W22-W36 S28-S25
      9  164  707 (N06E74 L=140)  3095    230        3 ( 122/ 587- 212/ 790) -90.0 E80-E70 S01-N11
     10  801  748 (S04E19 L=195)   293     83        0 ( 769/ 718- 835/ 783)  31.0 E21-E16 S05-S02
     11  409  674 (N00E47 L=167)  1703    299        3 ( 324/ 621- 515/ 771) -51.3 E52-E39 S03-N07
     12  613  672 (S04E33 L=181)   226     43        0 ( 579/ 656- 639/ 688) -14.0 E35-E31 S04-S03
     13 1070  628 (S18E08 L=206)   479    142        0 (1035/ 596-1138/ 661) -33.7 E09-E03 S19-S14
     14 1629  510 (S44W28 L=243)  1268    165        1 (1548/ 492-1708/ 530) -15.5 W20-W38 S47-S39
     15  881  429 (S25E26 L=188)  2219    435        2 ( 734/ 303-1064/ 516)   6.1 E41-E13 S30-S19
     16  460  405 (S17E57 L=157)   728    102        3 ( 433/ 361- 487/ 460)  60.3 E61-E54 S19-S12
     17  591  285 (S27E56 L=158)   723     68        1 ( 560/ 259- 628/ 310)   5.7 E62-E52 S29-S25
  

    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