Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  502 1241 (N30E28 L=148)   330     54        1 ( 476/1212- 528/1274) -63.4 E29-E25 N28-N32
      2  501 1132 (N24E30 L=146)   358     74        2 ( 482/1090- 526/1173)  59.0 E31-E27 N21-N26
      3  180  988 (N22E63 L=113)  1518    145        2 ( 152/ 921- 202/1063)  69.4 E67-E59 N18-N26
      4 1647  777 (S26W33 L=209)   272     87        0 (1619/ 734-1661/ 824)  51.3 W30-W34 S28-S23
      5 1399  756 (S20W15 L=191)   472    232        0 (1324/ 674-1484/ 818)  17.1 W08-W21 S23-S17
      6 1713  640 (S38W38 L=214)   403     32        1 (1693/ 620-1738/ 663) -38.7 W36-W41 S40-S35
      7 1311  527 (S32W04 L=180)   439    171        1 (1226/ 489-1373/ 562)   0.0 E02-W08 S33-S30
      8 1109  348 (S38E16 L=160)   352     76        0 (1068/ 330-1145/ 365)   0.0 E19-E12 S39-S37
  

    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