Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1 1597 1709 (N26W65 L= 48)   610    122        0 (1569/1666-1618/1749)  71.6 W63-W66 N23-N29
      2 1686 1421 (N09W55 L= 38)   783    162        2 (1650/1353-1721/1490)  77.5 W52-W57 N04-N13
      3 1245 1347 (N19W22 L=  5)   461    152        0 (1172/1312-1306/1397)   0.0 W16-W27 N17-N20
      4 1772  909 (S21W45 L= 28)   302     24        0 (1751/ 892-1792/ 924) -45.0 W43-W46 S22-S19
      5 1204  631 (S20E00 L=343)   212     72        0 (1191/ 592-1219/ 670)  68.2 E01-W01 S22-S17
      6 1164  530 (S25E05 L=338)   367     87        0 (1120/ 510-1207/ 563)   0.0 E07-E01 S26-S23
      7  933  310 (S34E29 L=314)   929    268        1 ( 879/ 230- 989/ 402)  32.9 E37-E20 S39-S28
      8 1133  160 (S56E29 L=314)  1147    208        0 (1045/ 123-1235/ 192) -23.6 E39-E17 S58-S51
  

    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