Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  502 1717 (N70E38 L=114)   931     91        3 ( 470/1679- 544/1744)   4.1 E49-E27 N66-N72
      2 1510 1602 (N25W53 L=205)  2181    339        6 (1352/1570-1667/1637) -24.8 W40-W66 N18-N31
      3 1618 1410 (N11W49 L=201)   322     64        0 (1587/1384-1647/1435) -63.4 W47-W51 N09-N13
      4  514 1274 (N34E28 L=124)   360    129        0 ( 451/1251- 569/1295) -26.6 E33-E23 N31-N36
      5  465 1147 (N27E33 L=118)   966    297        1 ( 419/1061- 554/1282) -66.4 E36-E27 N20-N36
      6  644 1137 (N23E20 L=131)   234     49        1 ( 619/1112- 670/1161) -45.0 E21-E18 N21-N24
      7  248  839 (N13E59 L= 93)   729    127        0 ( 205/ 815- 307/ 874) -15.3 E63-E52 N10-N15
      8  418  773 (N06E45 L=107)   564    150        1 ( 385/ 704- 443/ 826)  81.9 E46-E42 N02-N09
      9  688  716 (S03E27 L=124)   335     73        1 ( 661/ 683- 714/ 742) -63.4 E28-E26 S04-N00
     10  554  547 (S10E42 L=110)   342     60        0 ( 522/ 526- 587/ 569)   9.5 E45-E39 S10-S08
     11  412  468 (S12E59 L= 93)   758     99        0 ( 393/ 424- 435/ 511)  36.9 E62-E53 S15-S09
     12 1184  423 (S34E07 L=145)   365     80        1 (1148/ 395-1213/ 444) -26.6 E09-E04 S35-S31
  

    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