Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1410 1740 (N57W09 L=218)   825     88        2 (1376/1711-1445/1774) -24.0 W04-W13 N52-N59
      2 1188 1592 (N36E05 L=204)  1043    302        1 (1119/1456-1275/1704) -37.4 E13-W04 N28-N42
      3 1693 1485 (N44W41 L=250)   560     56        3 (1667/1455-1723/1513)   0.0 W38-W45 N41-N46
      4 1512  767 (S06W36 L=245)   300     40        1 (1493/ 739-1535/ 793) -56.3 W34-W37 S07-S04
      5  331  712 (S38E45 L=163)   336     72        1 ( 309/ 688- 373/ 728)  20.6 E48-E40 S40-S36
      6  396  556 (S48E40 L=169)   770     83        0 ( 369/ 515- 423/ 600)  48.8 E45-E36 S52-S44
      7  428  483 (S53E38 L=170)   555     52        1 ( 400/ 455- 461/ 511)  26.6 E43-E32 S55-S50
      8 1528  390 (S27W58 L=266)  2327    255        3 (1403/ 350-1644/ 443)  32.9 W48-W66 S32-S20
      9  516  398 (S58E28 L=181)   469     33        1 ( 494/ 373- 537/ 422)  51.3 E31-E24 S60-S55
  

    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