Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  840 1381 (N13E16 L=112)   224     63        2 ( 822/1349- 864/1416)  90.0 E16-E14 N11-N15
      2 1456  787 (S14W30 L=159)   423    144        0 (1419/ 715-1498/ 832)  69.4 W29-W32 S18-S10
      3 1715  644 (S17W55 L=183)   435     62        2 (1691/ 612-1741/ 675) -31.0 W52-W58 S18-S15
      4 1417  527 (S30W36 L=164)   262    106        0 (1389/ 494-1459/ 570)  71.6 W34-W37 S32-S26
      5 1545  512 (S28W47 L=175)   498    127        0 (1493/ 479-1600/ 565)  60.3 W43-W49 S31-S24
      6 1651  526 (S25W56 L=184)   500    103        0 (1607/ 492-1691/ 554)  45.0 W52-W58 S28-S22
      7 1424  466 (S34W39 L=167)   362     93        0 (1375/ 448-1472/ 489)   0.0 W34-W43 S34-S32
      8 1642  424 (S32W63 L=191)   561     96        0 (1584/ 409-1683/ 443)  15.3 W57-W68 S33-S30
      9 1320  173 (S57W62 L=190)   601     57        1 (1299/ 148-1343/ 204)  -9.5 W55-W73 S58-S54
  

    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