Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1094 1573 (N33W18 L= 75)  1661    361        3 ( 950/1535-1267/1612) -24.2 W07-W28 N27-N38
      2 1591 1532 (N19W51 L=108)  1368    240        0 (1486/1492-1671/1602)   2.6 W41-W63 N17-N22
      3  463 1421 (N38E32 L= 24)  1405    205        5 ( 385/1344- 534/1487) -40.6 E41-E25 N31-N44
      4  474 1191 (N23E32 L= 25)   366     68        1 ( 445/1162- 512/1209)   0.0 E34-E28 N20-N24
      5  157 1013 (N18E65 L=351)  6924    391       16 (  95/ 862- 224/1174)  78.1 E79-E55 N08-N27
      6  531 1090 (N16E29 L= 28)   661    124        2 ( 484/1035- 575/1132) -54.5 E31-E26 N11-N18
      7 1217  724 (S19W05 L= 62)   698    242        1 (1129/ 672-1334/ 790) -47.7 W01-W12 S24-S13
      8 1762  657 (S37W49 L=106)   922     85        3 (1725/ 614-1791/ 693) -78.7 W44-W53 S39-S33
      9  767  460 (S26E30 L= 27)   446    112        1 ( 747/ 409- 785/ 517)  54.5 E33-E27 S29-S22
     10  396  371 (S25E69 L=347)   631     48        0 ( 374/ 344- 414/ 399)  45.0 E72-E66 S26-S22
  

    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