Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  689 1729 (N26E45 L= 31)   640    117        0 ( 671/1672- 709/1777) -37.9 E50-E41 N22-N29
      2 1100 1615 (N30E13 L= 63)   755    272        0 (1015/1496-1181/1696) -25.2 E21-E04 N24-N32
      3  759 1605 (N20E33 L= 43)  1030    277        5 ( 681/1526- 841/1661) -17.4 E40-E24 N16-N22
      4 1324 1626 (N38W01 L= 77)   334    107        1 (1288/1575-1352/1657) -38.7 E01-W04 N35-N40
      5  156 1020 (S27E67 L=  9)  1596    149        6 ( 109/ 974- 223/1068)  -3.0 E76-E57 S29-S23
      6 1768  664 (S04W61 L=137)   375     21        0 (1752/ 642-1785/ 685) -26.6 W58-W62 S04-S02
      7 1500  590 (S17W43 L=119)   282     98        0 (1452/ 575-1546/ 616)  45.0 W39-W45 S18-S13
      8 1502  441 (S24W52 L=128)   489    107        0 (1459/ 404-1538/ 488)  90.0 W50-W53 S27-S20
      9 1644  455 (S18W63 L=139)   375     58        0 (1616/ 437-1672/ 470)  31.0 W59-W65 S19-S16
     10 1627  364 (S23W72 L=148)  2013    101        5 (1588/ 320-1668/ 407)  80.5 W67-W81 S25-S19
     11 1374  388 (S31W47 L=123)   347    114        0 (1310/ 377-1421/ 399)  29.7 W42-W49 S33-S29
  

    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