Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  500 1676 (N23E53 L= 41)   666    129        0 ( 475/1627- 521/1726) -54.5 E55-E50 N19-N26
      2  900 1512 (N19E18 L= 75)   507    214        1 ( 853/1419- 932/1597) -48.0 E23-E14 N13-N23
      3  251 1314 (S03E60 L= 33)  3124    488        6 ( 166/1096- 329/1505) -85.4 E70-E55 S16-N09
      4 1259 1291 (N13W08 L=101)   425    190        1 (1210/1226-1341/1374) -37.9 W04-W13 N10-N17
      5 1723 1297 (N26W40 L=133)   426     92        1 (1694/1247-1744/1337) -45.0 W37-W42 N23-N28
      6  764 1168 (S04E18 L= 75)   317    138        0 ( 711/1109- 809/1224)  69.4 E19-E16 S08-N00
      7  428 1153 (S10E40 L= 53)   261    103        0 ( 409/1129- 461/1203) -45.0 E41-E36 S11-S07
      8 1818 1001 (N11W53 L=146)   963    213        1 (1788/ 907-1841/1106) -57.5 W48-W57 N05-N16
      9  294 1055 (S18E49 L= 44)   319     80        0 ( 282/1017- 305/1099) -76.0 E50-E47 S19-S15
     10  740  837 (S23E13 L= 80)   227     92        0 ( 695/ 824- 793/ 851)  18.4 E16-E10 S24-S21
     11 1545  363 (S31W57 L=150)   626    169        1 (1485/ 343-1603/ 383)  29.1 W51-W60 S33-S28
  

    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