Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  576 1683 (N31E45 L=293)   806    154        4 ( 552/1617- 600/1746) -51.3 E50-E39 N26-N36
      2  778 1686 (N33E28 L=309)   404    141        1 ( 719/1647- 820/1733) -26.6 E35-E23 N30-N36
      3  644 1468 (N17E30 L=307)   341    118        2 ( 621/1411- 664/1522) -60.3 E32-E27 N13-N20
      4 1387 1437 (N24W18 L=355)   250    115        0 (1360/1389-1405/1487)  81.9 W16-W19 N20-N27
      5 1255 1329 (N15W10 L=347)   407    141        0 (1208/1270-1305/1373) -35.5 W06-W13 N12-N17
      6 1718 1322 (N23W44 L= 21)   401    125        0 (1699/1274-1749/1371) -63.4 W42-W46 N20-N26
      7 1624 1308 (N20W36 L= 13)   330     77        2 (1604/1280-1643/1351) -53.1 W34-W37 N18-N22
      8  591  919 (S17E26 L=311)   321     74        0 ( 569/ 875- 619/ 950)  63.4 E27-E24 S19-S14
      9  268  909 (S20E53 L=285)   341     44        1 ( 244/ 892- 295/ 931)   9.5 E55-E49 S21-S18
     10  199  847 (S24E62 L=276)   545     91        1 ( 168/ 809- 221/ 887) -51.3 E65-E58 S26-S21
     11  354  840 (S24E45 L=293)   460    119        1 ( 304/ 819- 401/ 873)  18.4 E49-E40 S25-S21
     12  824  764 (S24E09 L=328)   485    196        1 ( 771/ 687- 857/ 857)  74.7 E12-E07 S29-S18
     13 1712  652 (S17W55 L= 32)   405     40        1 (1688/ 629-1734/ 677)  56.3 W53-W57 S18-S14
  

    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