Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  614 1745 (N27E52 L=284)   875    109        1 ( 586/1701- 649/1789) -35.0 E57-E46 N23-N30
      2  424 1712 (N21E70 L=266)  1849    113        4 ( 381/1678- 462/1763) -66.8 E77-E64 N18-N26
      3  985 1462 (N17E14 L=323)   315     55        0 ( 966/1428-1000/1493) -53.1 E15-E11 N14-N18
      4 1687 1222 (N23W37 L= 13)  1224    319        2 (1633/1086-1740/1374) -47.1 W30-W43 N16-N31
      5  874  981 (S13E07 L=329)   223     37        0 ( 853/ 960- 903/ 998)  33.7 E08-E05 S14-S11
      6  857  764 (S26E03 L=333)   271    111        1 ( 831/ 718- 896/ 811)  74.1 E03-E01 S29-S22
      7  181  732 (S43E69 L=267)   865     99        0 ( 153/ 684- 216/ 768)  13.4 E82-E60 S45-S40
      8  934  617 (S33W06 L=342)   284     56        0 ( 911/ 597- 966/ 646) -14.0 W04-W08 S33-S30
      9  533  403 (S61E25 L=312)   442     38        1 ( 511/ 382- 555/ 433) -63.4 E28-E21 S62-S57
     10  693  334 (S61E03 L=334)   471     80        0 ( 660/ 305- 729/ 362)  -5.7 E07-W03 S62-S58
  

    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