Color  xc   yc  (LatLon LCarr)  Area Length Branches (  x0/y0     x1/y1  )  Tilt Location
      1  897 1518 (N28E17 L= 75)  1021    417        0 ( 806/1355-1009/1672) -38.0 E29-E06 N19-N37
      2 1543 1227 (N20W31 L=123)   411    116        1 (1506/1190-1585/1274) -29.7 W27-W34 N18-N22
      3 1889 1222 (N26W70 L=162)  1339     94        2 (1847/1187-1935/1245)  -2.5 W62-W86 N23-N27
      4 1769  912 (N05W54 L=146)   464     79        0 (1735/ 881-1800/ 935)  -8.1 W50-W57 N03-N05
      5 1567  706 (S11W41 L=133)   391    144        0 (1548/ 632-1590/ 768) -83.7 W40-W42 S15-S06
      6 1356  594 (S21W30 L=122)   232    108        1 (1348/ 536-1374/ 643) -66.8 W28-W31 S24-S17
      7 1096  458 (S35W16 L=108)   672    129        1 (1078/ 396-1130/ 524) -68.2 W13-W19 S39-S29
  

    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