Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1054 1878 (N57W49 L=315)   785     75        3 (1017/1848-1089/1907)  78.7 W44-W56 N53-N59
      2  755 1696 (N55W05 L=271)   438     99        1 ( 709/1664- 803/1717)   0.0 E01-W09 N53-N56
      3 1607  787 (S24W29 L=296)   674    182        0 (1523/ 739-1675/ 830)   0.0 W22-W35 S25-S23
      4 1318  689 (S21W08 L=274)   410    178        0 (1286/ 627-1376/ 726)  39.8 W04-W11 S24-S19
      5  437  683 (N00E45 L=221)   297     64        0 ( 413/ 664- 470/ 707) -33.7 E46-E43 N00-N02
      6 1374  624 (S27W10 L=276)   355     94        0 (1351/ 579-1396/ 674)  45.0 W07-W12 S29-S24
      7 1446  565 (S33W13 L=280)   422    101        0 (1395/ 548-1495/ 590) -14.0 W09-W17 S34-S31
      8 1705  512 (S48W39 L=305)   889    133        2 (1671/ 474-1733/ 571) -55.0 W34-W44 S52-S42
      9 1094  497 (S27E11 L=255)   473    120        1 (1066/ 446-1129/ 553)  35.5 E14-E07 S29-S23
     10  505  348 (S21E59 L=207)   660     53        2 ( 477/ 318- 533/ 380)  59.0 E61-E56 S23-S18
  

    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