Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  548 1755 (N43E46 L= 63)  1487    267        3 ( 513/1713- 596/1796) -69.4 E51-E41 N39-N47
      2  702 1693 (N38E29 L= 80)  1211    172        3 ( 637/1654- 774/1723)  -9.5 E34-E22 N35-N40
      3  806 1637 (N34E19 L= 90)   344     48        1 ( 782/1614- 832/1661)  56.3 E20-E16 N32-N35
      4  885 1617 (N33E13 L= 96)   629    155        7 ( 834/1586- 955/1652) -12.5 E16-E07 N30-N35
      5  325 1580 (N30E57 L= 51)   596    232        1 ( 310/1529- 338/1642) -52.1 E62-E54 N25-N34
      6  763 1421 (N19E18 L= 90)   214     84        1 ( 726/1389- 784/1457) -51.3 E21-E16 N16-N21
      7  619 1339 (N13E27 L= 82)   304     83        0 ( 607/1297- 633/1385) -80.5 E28-E25 N10-N16
      8  980 1020 (S05E03 L=106)   176     66        0 ( 955/ 984-1014/1040)  45.0 E03-W00 S06-S03
      9  809  696 (S25E13 L= 96)   811    271        0 ( 700/ 605- 940/ 790)  34.5 E20-E04 S31-S18
     10  253  681 (S27E57 L= 52)   595    138        1 ( 233/ 634- 279/ 726)  50.2 E59-E53 S29-S23
     11 1700  665 (S22W50 L=159)   297     90        0 (1673/ 648-1732/ 683)  -9.5 W47-W53 S22-S20
     12 1115  548 (S33W09 L=118)   489    120        1 (1037/ 530-1173/ 566)   0.0 W03-W12 S34-S32
  

    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