Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  537 1741 (N26E58 L=330)   776    113        2 ( 506/1689- 559/1787) -45.0 E63-E54 N21-N29
      2  424 1686 (N20E66 L=323)  1189     86        4 ( 382/1642- 454/1715) -56.3 E71-E61 N16-N22
      3 1143 1635 (N32E09 L= 19)   472    130        0 (1126/1572-1161/1701) -56.3 E12-E06 N28-N37
      4 1612 1423 (N33W28 L= 57)  1710    286        2 (1523/1299-1695/1527) -31.0 W20-W36 N26-N38
      5 1706 1216 (N22W39 L= 67)   376     55        0 (1671/1193-1735/1238)   9.5 W35-W41 N21-N23
      6  482  952 (S23E32 L=356)  1015    222        2 ( 399/ 899- 589/ 994)  34.7 E38-E25 S27-S18
      7  887  709 (S28E00 L= 28)   314    162        0 ( 823/ 681- 970/ 731)  40.6 E03-W04 S31-S25
      8  569  673 (S39E22 L=  6)   451    109        1 ( 531/ 628- 612/ 719)  58.0 E25-E19 S42-S34
      9  388  428 (S62E52 L=336)   525     88        0 ( 362/ 404- 435/ 458) -10.3 E57-E43 S63-S59
  

    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