Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1391 1597 (N43W23 L=117)   367     39        1 (1363/1576-1420/1613)  18.4 W20-W26 N41-N44
      2 1206 1369 (N25W08 L=102)  1137    478        3 (1063/1276-1354/1485) -24.2 E02-W18 N18-N30
      3  308 1445 (N19E62 L= 32)   645    143        1 ( 269/1395- 370/1475) -13.2 E69-E52 N16-N20
      4 1641 1416 (N33W44 L=138)   516    116        0 (1593/1370-1681/1460) -32.0 W39-W47 N29-N35
      5 1117 1243 (N16W04 L= 97)   219    113        0 (1104/1189-1130/1295) -80.5 W02-W04 N12-N18
      6 1033  719 (S17W04 L= 98)   454    137        2 ( 985/ 685-1102/ 755)  20.6 W00-W08 S19-S15
      7 1492  616 (S19W37 L=131)   488    108        1 (1440/ 595-1548/ 632)   7.1 W32-W41 S20-S18
      8  337  555 (S38E54 L= 39)  1418    177        2 ( 270/ 513- 419/ 596)  11.7 E72-E43 S42-S33
      9 1563  383 (S35W57 L=151)   735    122        1 (1511/ 348-1618/ 424)  45.0 W52-W60 S38-S31
     10 1047  202 (S59W19 L=113)   705    127        1 (1023/ 155-1062/ 260) -61.7 W14-W25 S65-S52
  

    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