Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  957 1720 (N47W12 L=112)   281     64        1 ( 926/1692- 980/1742) -36.9 W09-W13 N44-N49
      2 1630 1563 (N22W56 L=156)   794    116        3 (1576/1519-1678/1607) -45.0 W52-W59 N18-N25
      3  243 1458 (N43E62 L= 38)   544     52        1 ( 214/1435- 267/1475)   5.7 E68-E58 N40-N43
      4  571 1388 (N31E24 L= 76)   465    176        0 ( 519/1311- 631/1433)  36.9 E28-E18 N27-N33
      5 1821 1045 (S12W55 L=155)   379     92        0 (1787/1005-1847/1072) -45.0 W51-W57 S14-S09
      6 1549 1013 (S09W32 L=132)   354     81        1 (1528/ 966-1565/1055)  59.0 W30-W33 S11-S06
      7  321  870 (N04E49 L= 51)   877    146        1 ( 267/ 827- 372/ 932)  20.0 E54-E43 N02-N07
      8 1097  847 (S10W01 L=101)   670    181        1 (1048/ 763-1151/ 913)  41.2 E03-W05 S14-S06
      9  509  758 (S05E37 L= 63)   449     76        1 ( 487/ 719- 547/ 788)  18.4 E39-E33 S07-S03
     10 1830  708 (S33W58 L=158)   933    106        0 (1800/ 658-1858/ 751) -33.7 W53-W65 S37-S29
     11 1707  652 (S34W44 L=144)   789    205        2 (1635/ 578-1763/ 727)  32.5 W37-W49 S38-S30
     12  518  610 (S14E41 L= 59)   511     74        0 ( 492/ 586- 554/ 644) -76.0 E43-E37 S15-S11
  

    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