Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1681 1647 (N26W69 L=341)   530     58        0 (1649/1633-1711/1659) -14.0 W64-W74 N24-N26
      2  463 1605 (N49E36 L=235)   345     32        0 ( 440/1590- 491/1620)   0.0 E39-E33 N47-N50
      3 1144 1481 (N26W16 L=288)   652    129        1 (1101/1427-1205/1523)  24.0 W12-W21 N23-N27
      4  450 1428 (N36E34 L=238)   652    154        3 ( 372/1400- 518/1458) -14.9 E42-E27 N33-N38
      5 1709  997 (S13W44 L=315)   350     89        1 (1679/ 968-1745/1036)  26.6 W41-W47 S14-S10
      6 1841  833 (S25W57 L=329)   780    164        0 (1806/ 759-1876/ 908)  48.8 W53-W61 S29-S21
      7  340  512 (S17E61 L=210)  3580    371       10 ( 246/ 411- 402/ 639)  61.9 E69-E54 S24-S09
      8 1610  569 (S38W35 L=306)   462     76        1 (1585/ 534-1631/ 609) -59.0 W32-W37 S40-S35
  

    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