Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  638 1661 (N21E45 L=248)   543    177        0 ( 596/1613- 682/1701) -21.8 E49-E39 N18-N22
      2  280 1480 (N02E67 L=226)   985    173        2 ( 267/1423- 307/1539) -41.2 E71-E61 S01-N06
      3  510 1464 (N06E44 L=249)   286    134        0 ( 474/1422- 531/1498) -23.2 E48-E41 N03-N07
      4  154 1148 (S19E66 L=227)  2312    336        2 ( 114/ 996- 202/1268) -83.3 E69-E63 S28-S11
      5  193 1145 (S19E61 L=232)  1034    132        1 ( 159/1092- 221/1204)  81.9 E63-E57 S22-S15
      6  646 1062 (S14E23 L=270)   504    133        2 ( 583/1016- 693/1116)  60.3 E26-E20 S16-S09
      7  310 1022 (S24E46 L=247)   579    184        3 ( 245/1002- 368/1048)   5.2 E52-E41 S25-S22
      8  570  929 (S23E25 L=268)   560    229        0 ( 475/ 896- 652/ 972)  30.3 E31-E19 S26-S18
      9 1112  803 (S15W11 L=304)   300    109        0 (1095/ 753-1126/ 866) -66.8 W09-W12 S18-S11
     10 1066  676 (S24W13 L=306)   252     94        1 (1019/ 654-1119/ 703)   0.0 W09-W16 S24-S22
     11  367  603 (S50E39 L=254)   528     91        2 ( 346/ 569- 396/ 633)  68.2 E41-E35 S52-S47
     12  482  566 (S49E26 L=267)   974    226        1 ( 401/ 541- 570/ 599)  25.2 E35-E17 S52-S43
  

    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